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Abstract: In the lead-up to the Square Kilometre Array (SKA) project, several next-generation 
radio telescopes and upgrades are already being built around the world. These include APERTIF 
(The Netherlands), ASKAP (Australia), eMERLIN (UK), VLA (USA), e-EVN (based in Europe), 
LOFAR (The Netherlands), Meerkat (South Africa), and the Murchison Widefleld Array (MWA). 
Each of these new instruments has different strengths, and coordination of surveys between them 
can help maximise the science from each of them. A radio continuum survey is being planned on 
each of them with the primary science objective of understanding the formation and evolution of 
galaxies over cosmic time, and the cosmological parameters and large-scale structures which drive 
it. In pursuit of this objective, the different teams are developing a variety of new techniques, 
and refining existing ones. To achieve these exciting scientific goals, many technical challenges 
must be addressed by the survey instruments. Given the Umited resources of the global radio- 
astronomical community, it is essential that we pool our skills and knowledge. We do not have 
sufficient resources to enjoy the luxury of re-inventing wheels. We face significant challenges in 
calibration, imaging, source extraction and measurement, classification and cross-identification, 
redshift determination, stacking, and data-intensive research. As these instruments extend the 
observational parameters, we will face further unexpected challenges in calibration, imaging, and 
interpretation. If we are to realise the full scientific potential of these expensive instruments, it 
is essential that we devote enough resources and careful study to understanding the instrumental 
effects and how they will affect the data. We have established an SKA Radio Continuum Survey 
working group, whose prime role is to maximise science from these instruments by ensuring we 
share resources and expertise across the projects. Here we describe these projects, their science 
goals, and the technical challenges which are being addressed to maximise the science return. 

Keywords: telescopes — surveys — stars: activity — galaxies: evolution — galaxies: formation — 
cosmology: observations — radio continuum: general 



1 Introduction & Background 

The Square Kilometre Array (SKA) is a proposed ma- 
jor internationally-funded radio telescope ( |Dewdney| 
et al.||200"9 1 which is expected to be completed in the 



next decade. It will be many times more sensitive than 
any existing radio telescope, covering centimetre to me- 
tre wavelengths, and will answer fundamental ques- 
tions about the Universe by surveying the radio sky and 
studying individual objects in detail ( [Carilli & Rawl-| 
ings||2004 1. It wiU consist of many antennas, consti- 



tuting an effective collecting area of about one square 
kilometre, deployed over two sites in Australia and South 
Africa. 

In the lead-up to the SKA, several next-generation 
radio telescopes and upgrades are being constructed 
around the world, including APERTIF (The Nether- 
lands), ASKAP (Australia), eMERLIN (UK), e-EVN 
(based in Europe), LOFAR (The Netherlands), Meer- 
kat (South Africa), MWA (Australia) and the VLA (Karl 
G. Jansky Very Large Array, in the USA). Large con- 
tinuum surveys are being planned for many of these 
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telescopes. 

The predicted sensitivities and areas for these sur- 
veys are shown in Figure[T| alongside existing 1 .4 GHz 
continuum radio surveys. The largest published radio 
survey, shown in the top right, is the wide but shal- 



low NRAO VLA Sky Survey (NVSS) (Condon etal 



|1998[ ). The most sensitive published radio survey is the 



deep but narrow Lockman Hole observation ( Owen & 



[Morison 2008 ) in the lower left, which has recently 
been overtaken by an even deeper observation of the 
same field ( Condon et al.||20l"2 l. All current surveys 
are bounded by a diagonal line that roughly marks the 
Hmit of available telescope time of current-generation 
radio telescopes. The region to the left of this line 
is currently unexplored, and this area of observational 
phase space presumably contains as many potential new 
discoveries as the region to the right. 

Most of these projects have multiple science goals, 
but they have one goal in common, which is to survey 
the radio continuum emission from galaxies, in order 
to understand the formation and evolution of galax- 
ies over cosmic time, and the cosmological parame- 
ters and large-scale structures that drive it. In pursuit 
of this goal, the different teams are developing tech- 
niques such as multi-scale deconvolution, source ex- 
traction and classification, and multiwavelength cross- 
identification. Furthermore, these projects share spe- 
cific scientific goals, some of which require further def- 
inition before a well-planned survey can be executed. 
Finally, each of these new instruments has different 
strengths, and coordination of surveys between them 
can help maximise the scientific return from each. 

The radio continuum surveys planned for these in- 
struments typically reach orders of magnitude deeper 
than traditional surveys. Consequently, the resulting 
surveys will be far more than just a deep version of 
existing surveys, but will be qualitatively different, in 
most cases being dominated by star-forming galaxies. 
For example, only a total of about 10 deg^ of the sky 
has been surveyed at 1 .4 GHz to the planned 10 /iJy/beam 
rms of many of these surveys, in fields such as the Hub- 
ble, Chandra, ATLAS, COSMOS and Phoenix deep 
fields (Hopkins et al. 2003; Huynh et al. 2005; Biggs 



& Ivison,|2006, ,Norri s et al. 2006; Middelberg e t al. 
2008al|Milleret al.|2008 ; Schinnerer et al. 2007 , Mot " 



rison et al.'2010'; TTales et al. |2012b| |Banfield et al 
2012; Guglielmino et al. 2012| ). 



Surveys at this depth extend beyond the traditional 
domains of radio astronomy, where sources are pre- 
dominantly radio-loud galaxies and quasars, into the 
regime of star-forming galaxies. At this depth, even the 
most common active galactic nuclei (AGN) are radio- 
quiet AGNs, which make up most of the X-ray extra- 
galactic sources. As a result, the role of radio astron- 
omy is changing. Previous radio-astronomical surveys 
had their greatest impact in the niche area of radio- 
loud AGNs. For example, the numbers of sources in 
NVSS ( [Condon et aL|1998^ and FIRST ( [Becker et al.| 



[1995[ ) are overwhelmingly dominated by AGNs, and 
only a small fraction of their sources are star-forming 
galaxies. In contrast, most sources detected by the 
next-generation surveys are star-forming galaxies, so 
that these surveys are dominated by the same galax- 
ies as are studied by optical and infrared (IR) surveys. 
As a result, these next-generation radio surveys are an 
increasingly important component of multiwavelength 
studies of galactic evolution. 

These surveys, however, will also encounter a new 
set of technical and scientific challenges. For example, 
the new continuum surveys such as EMU (see 
andWODAN(see 
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should cover the whole sky at 
sensitivities previously attained only in small fields 
10 /iJy/beam rms) by instantaneously covering wide 
fields (- 30 deg2) with large bandwidths (- 300 MHz) 
and long integration times (^12 hours), and not by 
virtue of larger collecting areas or lower system noise 
temperatures. To reach the sensitivities predicted by 
the radiometer equation, these surveys must have ex- 
ceptionally high dynamic ranges (e.g., about 42 db for 
EMU) and low systematic errors. This requires very 
high quality hardware, calibration techniques, and imag- 
ing algorithms. In particular, the primary beams must 
have < 1% uncertainties in reconstructed position, size, 
and shape. The ASKAP dishes, for example, appear to 
have good surface and pointing accuracies, and their 
3-axis mounts maintain a constant footprint of the pri- 
mary beam on the sky. The electronically formed pri- 
mary beams of phased-array feeds also have to be very 
stable, and the survey astronomers need to work closely 
with the engineers to ensure that amplitude and phase 
calibration on strong sources in the array field-of-view 
can be transferred from one primary beam to another. 
Finally, the success of a survey in the astronomical 
community depends as much on the survey results be- 
ing made available to everyone in an easily-accessible 
fashion (such as through a web portal, or virtual ob- 
servatory tools, and with useful and intuitive query en- 
gines) as on whether the survey actually meets every 
performance goal, such as sensitivity. 

Moreover, some of these projects differ from ear- 
lier surveys in that one goal is to cross-identify the de- 
tected radio sources with major surveys at other wave- 
lengths, and produce public-domain Virtual Observa- 
tory (VO) accessible catalogues as "value-added" data 
products. This is facilitated by the growth in the num- 
ber of major surveys spanning all wavelengths, dis- 
cussed below in 



4.6 



To address these challenges, a new group, the SKA 
PAthfinder Radio Continuum Survey (SPARCS) Work- 
ing Group has been established, with the following goals: 

• To coordinate developments of techniques, to avoid 
duphcation of effort and ensure that each project 
has access to best practice. 

• To hold cross-project discussions of the specific 
scientific goals, to ensure cross-fertilisation of 
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ideas and optimum survey strategies. 

• To coordinate the surveys in their choice of area, 
depth, location on the sky, and other survey pa- 
rameters, to maximise the scientific return from 
the surveys. 

• To distil the SKA pathfinder experiences in order 
to provide the most relevant and up to date input 
for SKA planning. 

SPARCS was initially proposed during the SKA 
meeting in Manchester in March 2010, was formed 
in May 2010, and held its first meeting at the Lorentz 
Center in Leiden, The Netherlands, in February 2011. 
This review paper documents the outcomes of that work- 
shop, together with related developments since the work- 
shop. The authors are the speakers and facilitators of 
the workshop, together with other domain experts who 
were unable to attend the workshop in person. 

This paper describes the SKA pathfinder contin- 
uum surveys in § 2, setting out their scientific goals 
in § 3, and identifying the challenges to achieve these 
goals in § 4. § 5 concludes by summarising how these 
challenges are being addressed, and how SPARCS is 
helping to coordinate and facilitate this process. We re- 
fer to the collection of all surveys discussed in this pa- 
per as the SPARCS surveys, and the sources detected in 
those surveys as the SPARCS sources. We assume the 
convention for spectral index a that S (x i/" through- 
out. 



2 Pathfinders and their Continuum 
Surveys 

2.1 APERTIF-WODAN 

APERTIF, the new Phased Array Feed (PAF) receiver 
system for the Westerbork Synthesis Radio Telescope 
(WSRT) will dramatically enlarge the instantaneous field 
of view (FOV) of the WSRT ( |Oosterloo et al.|2009| by 
replacing the current single front-end feeds by PAFs. 
Each of the PAFs consists of 121 Vivaldi elements and 
will detect the radiation field (in dual polarisation) in 
the focal plane of each dish over an area of about one 
square metre at an observing frequency of 1.4 GHz. 
Because of this, many beams can be formed simulta- 
neously for each dish making it possible to image an 
area of about 8 square degrees on the sky, which is an 
increase of about a factor of 30 compared to the cur- 
rent WSRT. Its large 300 MHz bandwidth will not only 
cater for sensitive continuum imaging, but is also cru- 
cial for efficient HI and OH emission surveys and for 
studies of polarised emission from large areas. 

The WODAN (Westerbork Observations of the Deep 
APERTIF Northern-Sky) survey ( [Rottgering et al.|2010b| l 
will use APERTIF to survey the northern 25% of the 
sky (i.e. North of dechnation +30°) that is inaccessible 
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Figure 1 Comparison of existing and planned deep 
1.4 GHz radio surveys. The horizontal axis shows the 
5-0- sensitivity, and the vertical axis shows the sky cov- 
erage. The diagonal dashed line shows the approx- 
imate envelope of existing surveys, which is largely 
determined by the availability of telescope time. For 
example, to extend NVSS to the sensitivity of EMU 
would have required over 600 years of (pre-upgrade) 
VLA time, so in practice this would not have been pos- 
sible, and this line therefore represents a hard limit to 
the sensitivity of traditional surveys. The squares in the 
top-left represent the new radio surveys discussed in 
this paper. Surveys at other frequencies are not shown 
in this diagram, as their relative sensitivity depends on 
the assumed spectral index of the sources, although we 
do include SUMMS at 843 MHz, without making any 
correction for spectral index. A similar comparison of 



low-frequency surveys can be found in Tingay et al. 
( |20T2l l. 



to ASKAP, to a target rms sensitivity of 10 /xJy/beam 
at a spatial resolution of 15 arcsec, although confusion 
noise with a 15- arcsecbeam may increase the observed 
rms noise level to about 20 /iJy/beam. 



2.2 APERTIF-BEOWULF and APERTIF- 
FRIGG 

BEOWULF and FRIGG are both polarisation projects 
proposed for APERTIF. 

The BEOWULF (B -field Estimation and Observa- 
tional Wide-field Understanding of Large-scale Faraday- 
structure) survey ( Scaife et al.|2010 1 will use APER- 
TIF to survey an area of 64 deg'^ within the Perseus- 
Pisces super-cluster of galaxies in order to construct 
a finely spaced rotation measure (RM) grid tuned for 
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cosmological and large-scale structure studies of mag- 
netism within both clusters and the filamentary inter- 
cluster medium. This survey will complement other 
wide area polarisation surveys: by measuring RMs to 
5,000 background sources it will give an RM-grid sam- 
pling of 3.5 arcmin. The BEOWULF survey will pro- 
vide vital insights into whether magnetic fields were 
primordial in clusters of galaxies, or whether the galax- 
ies themselves later injected the field. The survey is 
ideally suited to APERTIF whose wide FOV enables 
observations of large scale structures, and whose fre- 
quency coverage is well matched to the comparatively 
high RMs expected from such a region. 

FRIGG (Faraday Rotation Investigation of Galax- 



ies and Groups) survey ( Beck et al. 2010 1 will observe 



large areas around 4 large northern nearby spiral galax- 
ies, the bridge between M31 and M33, and 5 galaxy 
groups. Its goal is to measure grids of Faraday rota- 
tion measures of background sources and faint diffuse 
emission, in order to unravel the extent and dynamical 
importance of cosmic magnetic fields in galaxies and 
in intergalactic space in galaxy groups. 

2.3 ASKAP-EMU 



The Australian SKA Pathfinder (A SKAP: [Johnston et| 
|aD200T 2008 ; Deboer et al. 2009 ) is a new radio tele- 
scope being built on the Australian SKA site in West- 
ern Australia, at the Murchison Radio-astronomy Ob- 
servatory, with a planned completion date of 2013-14. 
It will consist of 36 12-metre antennas distributed over 
a region 6 km in diameter. Each antenna is equipped 
with a PAF ( [Bunton & HaypOl O") of 96 dual-polarisation 
pixels operating in a frequency band of 700-1 800 MHz. 
As a result, ASKAP will have a field of view up to 30 
deg^. To ensure good calibration, the antennas are a 
novel 3-axis design, with the feed and reflector rotating 
to mimic the effect of an equatorial mount, ensuring a 
constant position angle of the PAF and sidelobes on the 
sky. 

In continuum mode, ASKAP will observe a 300 
MHz band, split into 1 MHz channels, with full Stokes 
parameters measured in each channel. The data will 
be processed in a multi-frequency synthesis mode, in 
which data from each channel are correctly gridded in 
the uv plane. As well as producing images and source 
catalogues, the processing pipeline will also measure 
spectral index, spectral curvature, and aU polarisation 
products across the band. All data processing steps, 
from the output of the correlator to science -quality im- 
ages, spectra, and catalogues, are performed in auto- 
mated pipelines ( Cornwell et al.|[20TT I running on a 
highly distributed parallel processing computer. These 
steps include flagging bad data, calibration, imaging, 
source-finding, and archiving. 



Southern Sky, extending as far North as +30°. EMU 
will cover roughly the same fraction (75%) of the sky 
as the benchmark NVSS survey ( |Condon et al.|1998| l, 
but will be 45 times more sensitive, and will have an 
angular resolution (10 arcsec) 4.5 times better. Be- 
cause of the excellent short-spacing uv coverage of AS- 
KAP, EMU will also have higher sensitivity to extended 
structures. Like most radio surveys, EMU will adopt 
a 5-(T cutoff, leading to a source detection threshold 
of 50 /iJy/beam. EMU is expected to generate a cata- 
logue of about 70 million galaxies, and all radio data 
from the EMU survey will be placed in the public do- 
main as soon as the data quality has been assured. 



Together, EMU and WODAN (see SjZTji will pro- 
vide full-sky 1.3 GHz imaging at ^ 10-15 arcsec res- 
olution to an rms noise level of 10 /xJy/beam, pro- 
viding an unprecedented sensitive all-sky radio survey 
as a legacy for astronomers at all wavelengths. The 
EMU and WODAN surveys will overlap by a few de- 
grees of declination to provide a comparison and cross- 
validation, to ensure consistent calibration, and to check 
on completeness and potential sources of bias between 
the surveys. 

2.4 ASKAP-POSSUM 

POSSUM (Polarisation Sky Survey of the Universe's 
Magnetism) is an all-sky ASKAP survey of linear po- 
larisation (Gaensler et al. 2010| l. It is expected that 
POSSUM will be commensal wifli EMU, and that the 
two surveys will overlap considerably in their analysis 
pipelines and source catalogues. POSSUM will pro- 
vide a catalogue of polarised fluxes and Faraday ro- 
tation measures for approximately 3 million compact 
extragalactic sources. These data will be used to de- 
termine the large-scale magnetic field geometry of the 
Milky Way, to study the turbulent properties of the in- 
terstellar medium, and to constrain the evolution of 
intergalactic magnetic fields as a function of cosmic 
time. POSSUM will also be a valuable counterpart to 
EMU, in that it will provide polarisation properties or 
upper limits to polarisation for all sources detected by 
EMU. 

2.5 e-MERLIN 



The e-MERLIN aiTay (Garrington et al. 2004 1, op 



The Evolutionary Map of the Universe (EMU: Nor 



|ris et al.||2011b| l will use ASKAP to make a deep (10 
yuJy/beam rms) radio continuum survey of the entire 



erated by the University of Manchester, is a signifi- 
cant upgrade to the existing telescopes which form the 
MERLIN array in the UK. Consisting of seven tele- 
scopes, spread across the UK with a maximum base- 
line of 2 1 7 km, e-MERLIN provides high angular reso- 
lution (10-150 mas) imaging and spectroscopy in three 
broad cm-wavebands (1.3-1.8 GHz, 4-8 GHz and 22- 
24 GHz). The e-MERLIN project includes a major up- 
grade to the existing telescope hardware, including the 
installation of new receivers, analogue and digital elec- 
tronics, optical- fibre links to each telescope and a pow- 



www.publish.csiro.au/joumals/pasa 



5 



erful new digital correlator at Jodrell Bank Observa- 
tory (JBO). This upgrade increases the usable band- 
width by more than two orders of magnitude, com- 
pared to the old MERLIN system, so that the contin- 
uum sensitivity is increased by a factor of between 10 
and 30. In addition, the increase in bandwidth will en- 
able multi-frequency synthesis, dramatically improv- 
ing the uv coverage for continuum observations and 
enabling simultaneous spectral-index imaging. 

The e-MERLIN upgrade is well underway with new 
receivers in service, the dedicated optical fibre network 
and the digital transmission equipment operational and 
returning 30 Gbit/s of data from each telescope back to 
a new correlator. Currently e-MERLIN is undergoing 
the final stages of instrument commissioning and will 
be fully operational during 2012. 

These upgrades to the system allow e-MERLIN to 
typically provide /iJy sensitivity with up to 4 GHz of 
instantaneous bandwidth. This sensitive continuum imag 
ing can be simultaneously combined with powerful spec- 
tral line and polarisation capabilities provided by the 
new highly flexible WIDAR correlator which is a twin 
of the new VLA correlator and has been developed by 
DRAG in Canada. 

The e-MERLIN array is an important SKA pathfi- 
nder instrument covering a unique scientific niche due 
to its baseline lengths, which are intermediate between 
those provided by the VLA and VLBI arrays. Techni- 
cally this provides a test-bed for the development of 
fibre-based data transport and time distribution over 
distances spanning many hundreds of kilometres. Sci- 
entifically, results stemming from its high angular res- 
olution will provide significant input into the eventual 
locations of dish receptors as part of SKA Phase 1 . 

e-MERLIN is an open-user instrument with twice 
yearly proposal deadlines. As with many other new ra- 
dio facilities, however, it has pre-allocated a large frac- 
tion of observing time to large key-science projects. 
Eleven large e-MERLIN 'legacy' project^ have been 
allocated a total of 5000 hours (approximately 50% of 
the available observing time) over the first 5 semesters 
of operation. These key projects span a wide range of 
astrophysics from studies of planet formation, pulsar 
proper motions, stars, nearby normal and star-forming 
galaxies, classical radio galaxies, strong gravitational 
lenses, galaxy clusters, and the physics of high redshift 
star-formation and AGN dominated galaxies. In each 
of these programmes the strengths of e-MERLIN, in 
terms of its sensitivity and critically high angular reso- 
lution at centimetre wavelengths, are 

being exploited to provide a spatially resolved view 
of of these sources, complementing studies planned or 
underway with other SKA pathfinder instruments, and 
new facilities at other wavelengths. 



2.6 VLA 



The Karl G. Jansky Very Large Array (VLA: |Napier| 
2006 ) represents a major upgrade to the VLA. Although 
the upgraded VLA was not designed to be a survey 
telescope or an SKA pathfinder, its high sensitivity, 
continuous frequency coverage, flexible wideband cor- 
relator, and scaled array configurations make the VLA 
a useful complement to the dedicated pathfinders, and 
to the all-sky continuum surveys that are being planned. 
For example, the VLA could observe "reference fields" 
covering a few small areas of sky with higher angular 
resolution (e.g. 5 arcsec) and sensitivity (e.g. 5 /iJy/beam 
rms) than the L4GHz EMU and WODAN surveys in 
order to join those surveys seamlessly and provide a 
reference for resolution corrections, flux densities, po- 
sitions, etc. of the faintest survey sources. 

The VLA is a general-purpose user instrument, open 
to researchers around the world. The VLA upgrade 
is approaching completion, and algorithms for making 
20 cm images with large fractional bandwidth and high 
dynamic range are now being tested on real data. The 
ultimate continuum sensitivity limits (source confusion 
or dynamic range) of the full SKA depend on the sur- 
face density on the sky of sub-/iJy sources. If dynamic 
range permits, by imaging a single field at 20 cm to the 
confusion limit with 5 arcsec resolution, the VLA in its 
B configuration will place a tight statistical constraint 
on the source density of objects fainter than 1 /iJy. 

2.7 eVLBI 

The SKA will be a real-time instrument with prefer- 
ably 25% of its collecting area forming long baselines 
that extend up to to 3000 km. Data caching and trans- 
portation, distribution of the clock signal, and opera- 
tions of this array will present a great technical chal- 
lenge. The European VLBI Network (EVN) as an SKA 
pathfinder is addressing these issues within the frame- 
work of the NEXPReS project ( htt p://www.nexpres.eu/) l. 
The primary goal of NEXPReS is to support and de- 
velop real-time electronic-VLBI (e-VLBI) operations 
in the EVN and on global scales. 

Science observations with the e-EVN have routinely 
been carried out since 2006. The possibility of operat- 
ing a real-time VLBI array with baselines exceeding 



12,000 km was recently demonstrated by iGiroletti et 



al. (201 1 1. Relaxing the data storage limitations at the 
telescopes creates an obvious advantage for the EVN. 
This allows for more flexible operations to carry out 
transient science (see §3.7i, it may provide seamless 



data transport at > 1 Gbps data rates for superior sensi- 
tivity, and makes it possible to conduct automated ob- 
servations for efficient VLBI surveys, or to respond to 
external triggers. 

The importance of long baselines (and even the pos- 
sibility of baselines to space) for the SKA have been 



'http://www.e-merlin.ac.uk/legacy/ 



discussed by |Fomalont & Reid| ( |2004[ | and |Gurvits 
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( |2004| l. Present day observations of faint, compact ra- 
dio sources are beginning to shed light on the science 
that could be done with one or two orders of magnitude 
better sensitivity. The mJy sky is thought to be domi- 
nated by AGN, compact on milliarcsec scales, while in 
the sub-mJy regime star-forming galaxies become in- 
creasingly important, and a number of objects are now 
known which show both strong AGN activity and vig- 
orous star-formation activity (e.g. Afonso et al.|2001 



|Perez-Torres"eraLl[20T()| |Norris et al.||2012| l. VLBI. 



and the long baselines of SKA, will play a key role in 
distinguishing between AGN and SF galaxies. 

The long baselines of the SKA will be particularly 
important for studying the details of feedback processes 
between AGN and star-formation ( |Croton et al.|2006l l, 
the role of binary AGN mergers in galaxy evolution 



and feedback (e.g. Dotti, Sesana & Decarli^2012j , the 
first supermassive black holes at the highest redshifts 
( |Falcke et al.][2004l [Prey et al.|[20T0l [MTT) and sev- 
eral other applications (see Godfrey et al. (2012 1 for a 
review). 

In traditional VLBI observations, the field of view 
was severely limited by time- and bandwidth smear- 
ing. Since the sub-mJy sky is densely populated with 
radio sources, VLBI observations are moving from the 
traditional mode of observing single, bright, widely- 
separated objects to a mode which allows observations 
of wide fields. This transition poses a significant tech- 
nological challenge, because the combination of high 
resolution and wide fields requires high spectral and 
temporal resolution in correlation, to prevent loss of 
coherence arising from averaging effects (bandwidth 
and time smearing). Early wide-field VLBI experi- 
ments (e.g. [Garrett, Wrobel & Morganti|2005] ), were 
limited by the hardware VLBI correlators, and by lim- 
ited computing resources. Recently, software correla- 
tors running on cheap commodity computers have pro- 
vided the power and flexibility needed for VLBI obser- 
vations of wide fields ( [Deller et al.|201 1) . 

Since sources detectable with VLBI, even in the 
sub-mJy regime, are relatively sparsely distributed, it 
has proven to be more efficient to generate multiple 
phase centres with small fields of view directed to- 
wards known locations of radio emission, rather than 
to image huge blank portions of the sky. Alternative, 
equivalent approaches have also been proven to be suc- 
cessful ( [Morgan et al.|2011| l. The multi-phase centre 
technique has only recently been brought to produc- 
tion (Middelberg et a l.|201 la I but has already demon- 
strated that galaxies can contain a radio-active AGN 
even when an AGN is not evident in comprehensive 



radio, optical, and X-ray observations (Middelberg et 
|al.|2012| l. 

Wide-field VLBI observations require new calibra- 
tion steps to reach the full sensitivity of the data, such 
as primary beam corrections and multi-phase centre 
self-calibration, which help to make the technique use- 
ful for a wide range of applications. For example, the 



multi-phase centre technique can be used to image many 
sources around a faint target to aid in self-calibration. 
While no single source may be sufficiently strong for 
self-calibration, the ensemble of sources can provide a 
robust solution. The technique may also be used to pro- 
vide an astrometric reference, or to image many back- 
ground sources to probe a foreground object. 

2.8 LOFAR 

LOFAR, the Low Frequency Radio Array, is a pan- 
European radio phased-array telescope that is currently 
being commissioned. The two types of antenna, one 
optimised for 30 - 80 MHz and the other for 110 - 240 
MHz, are grouped together in stations the size of soc- 
cer fields. 40 stations are distributed over an area of 
diameter of 100 km in The Netherlands, and a further 
eight stations are located in Germany, UK, Sweden, 
and France. The signals from the antennas are digi- 
tised to form many beams on the sky, making LOFAR 
an extremely efficient survey instrument. LOFAR has 
already generated images that are the deepest ever at 
these low frequencies. 

A key motivation of LOFAR is to provide the entire 
international astronomical community with surveys of 
the radio sky that have a long-lasting legacy value for 
a broad range of astrophysical research. The LOFAR 
continuum survey ( Rottgering et al.|2010a| l will cover 
the northern half of the sky (i.e. North of declination 
0°). LOFAR will be especially complementary to WOD- 
AN and EMU in surveying the sky at high sensitivity 
and resolution but at a much lower frequency. 

The three fundamental areas of astrophysics that 
have driven the design of the planned LOFAR surveys 
are: (i) the formation of massive galaxies at the epoch 
of reionisation, (ii) magnetic fields and shocked hot gas 
associated with the first bound clusters of galaxies, and 
(iii) star formation processes in distant galaxies. The 
areas, depths and frequencies of the surveys have been 
chosen to contain: (i) ^100 powerful radio galaxies 
close to or at the epoch of reionisation, (ii) ^100 ra- 
dio halos at the epoch when the first massive bound 
galaxy clusters appeared, and (iii) ^100 protoclusters. 
The resulting survey parameters are based on estimates 
of luminosity functions for powerful radio galaxies by 



Wilman et al. (2008), for radio halos by EnBlin & 



Rottgering (2002) and Cassano et al.| (i2010b) , and for 
protoclusters by Venemans et al. ( 2007 1. 

To achieve the goals of the LOFAR surveys, a three- 



tiered approach has been adopted (Rottgering et al. 
|2010a| l. Tier 1 is a 27r steradian survey reaching an 
rms of 0.07mJy at 15-65 and 120-180 MHz, and is 
designed to detect ^ 100 cluster halos at z > 0.6 and 
^ 100 z > 6 radio galaxies. Tier 2 is a deep survey 
over 500 deg^ at 30, 60, and 150 MHz. At 150 MHz, 
this will require 55 pointings, to be centred on the fol- 
lowing science targets: 

• 25 weU-studied extragalactic-fields with superb 
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existing multiwavelength data; 

• 15 fields centred on clusters or super-clusters; 

• 15 fields centred on nearby galaxies. 

Tier 3 is an "ultra-deep" survey at 150 MHz cov- 
ering a single pointing of 100 deg^ reaching the con- 
fusion level of 7 /xJy/beam rms. Figure |2] shows the 
resulting depth versus frequency. 



2.9 Meerkat/Mightee 
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Figure 2 Flux limits (5 a) of the proposed LOFAR sur- 
veys compared to other radio surveys. The triangles 
represent existing surveys. The lines represent differ- 
ent power-laws (with a — —1.6 and —0.8) to illustrate 
how, depending on the spectral indices of the sources, 
the LOFAR surveys will compare to others. 

WODAN and EMU have an enormous synergy with 
the LOFAR surveys: virtually all the 5 x 10^ radio 
sources from the LOFAR all sky surveys will have their 
flux density at 1400 MHz measured (i.e. Figure[T]). To- 
gether they will yield data for all radio loud AGN in 
the northern hemisphere, and most luminous starbursts 
up to z = 2. The resulting densely populated radio 
colour-colour diagrams will be a powerful tool to spec- 
trally discriminate between very rare radio sources with 
extreme radio spectra, such as diffuse emission from 
clusters and very distant radio galaxies. Nearby re- 
solved sources will have spectral index and spectral 
curvature maps, a very rich source of information to 
constrain many physical parameters. As the combined 
surveys will cover the entire sky, measurements of the 
Integrated Sachs-Wolfe effect, galaxy auto-correlation 
functions and cosmic magnification will substantially 
tighten cosmological model parameters ( [Raccanelli et| 
lafpOl l). 

LOFAR will also support a Key Science Project on 
Cosmic Magnetism 

(MKSif]), 

which is likely to make 
significant contributions to our knowledge of the mag- 
netic sky (see §|3.5[). 



MeerKAT ( Jonas 2009 1 is the South African SKA path- 



finder telescope. Meerkat Phase 1 will consist of 64 
dishes, each 13.5 m in diameter, equipped with single- 
pixel receivers. An offset Gregorian dish configuration 
has been chosen because the unblocked aperture pro- 
vides good optical performance, sensitivity and imag- 
ing quality. It also facilitates the installation of multi- 
ple receiver systems in the primary and secondary fo- 
cal areas. MeerKAT's 64 dish array layout will be dis- 
tributed over two components. A dense inner compo- 
nent will contain 70% of the dishes with a Gaussian uv- 
distribution with a dispersion of 300m. The outer com- 
ponent will contain the remaining 30% of the dishes, 
having a Gaussian Mv-distribution with a dispersion of 
2500 m and a longest baseline of 8 km. A potential 
future extension (Meerkat Phase 2) could see 7 addi- 
tional dishes being added to extend the longest base- 
lines to about 20 km. MeerKAT will support a wide 
range of observing modes, including deep continuum, 
polarisation and spectral line imaging, pulsar timing, 
and transient searches. The plan is to provide standard 
data products, including an imaging pipeline. A num- 
ber of "data spigots" will also be available to support 
user-provided instrumentation. 

70% of observing time on MeerKAT, for the first 
few years, is allocated for large survey projects of 1000 
hours or more, while the remaining 30% is reserved for 
smaller PI driven proposals (of which 5% will be di- 
rector's discretionary time). Proposals for key projects 
were solicited in 2010, and 10 projects were selected. 
Two (a pulsar timing and a deep HI survey) were cho- 
sen as top priority, and a further eight as second prior- 
ity, one of which is MIGHTEE. 



The MIGHTEE survey (van der Heyden & Jarvis 



http://www.mpifr-bonn.mpg.de/staff/rbeck/MKSP/mksp.html 



|2010[ ) aims to probe to much fainter flux densities (0. 1- 
1 ^Jy/beam rms) than the EMUAV ODAN surveys, over 
smaller areas (^ 35 deg^) at higher angular resolution, 
at an observing frequency of 1.4 GHz. The higher sen- 
sitivity and resolution will enable exploration of AGN, 
star-forming galaxies and galaxy clusters from the epoch 
of reionisation through to the present day. The MIGH- 
TEE survey strategy will follow a 3 -tiered approach 
both in survey area and in sensitivity, using the longer 
baselines to be able to probe below the confusion limit 
of the shorter baseline surveys such as EMU and WODAN. 

Tier 1 will cover 1000 deg^ using Meerkat Phase 1 
at 1 .4 GHz to a target rms flux-density of 5 /iJy/beamThis 
will allow detailed studies of the distant Universe and 
the evolution of the lower-luminosity radio source pop- 
ulations into the epoch of reionisation. Tier 2 will ex- 
ploit the longer baselines of Meerkat Phase 2 and in- 
creased bandwidth to observe a single 35 deg^ point- 
ing down to 1 /zJy/beam rms. This will enable stud- 
ies of the evolution of Milky Way type galaxies up to 
z ~ 4 and pioneering weak lensing analyses at radio 
wavelengths. These observations wiU also be used to 
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test techniques to increase dynamic range in prepara- 
tion for the full SKA. Tier 3 will cover a much smaller 
area (1 deg^) to 0.1 /iJy/beam rms. This will allow 
us to isolate the level of AGN activity in star-forming 
galaxies and to investigate the morphological proper- 
ties of AGN from z ~ 0.5 to z ~ 6. 



2.10 Murchison Widefield Array 



The Murchison Widefield Array (MWA; [Lonsdale et 
|al.|2d09 ; Tingay et al. 2012[ l is a low-frequency synthe- 
sis telescope under construction at Boolardy, Western 
Australia, adjacent to ASKAP. It is due for comple- 
tion in November 2012, making it the first of the SKA 
precursor telescopes to commence full operation. It 
consists of 2048 dual-polarisation dipole antennas, ar- 
ranged as 128 "tiles", each consisting of a 4 x 4 array 
of dipoles designed to operate in the 80-300 MHz fre- 
quency range. Each tile performs an analogue beam- 
forming operation, narrowing the field of view to an 
electronically steerable 25 degrees at 150 MHz. The 
majority of the tiles (112) are scattered across aroughly 
1.5 km core region, forming an array with very high 
imaging quality, and a field of view of several hundred 
square degrees at a resolution of several arcmin. The 
remaining 16 tiles are located outside the core, yielding 
baselines up to 3 km to allow higher angular resolution. 
The MWA calibration and imaging is undertaken on 
24 IBM iDataPlex dual Xeon servers, each housing 2 
NVIDIA Tesla M2070 Graphics Processing Units. The 
correlated data are calibrated in real time using novel 
position-dependent self-calibration algorithms. 

MWA is focused on four key science projects ( Bow| 
|man et aT]|2012| l: (a) the detection and characterisa- 
tion of 3-dimensional brightness temperature fluctua- 
tions in the 21cm line of neutral hydrogen during the 
Epoch of Reionisation (EoR) at redshifts from 6 to 10, 
(b) solar imaging and remote sensing of the inner helio- 
sphere via propagation effects on signals from distant 
background sources, (c) high-sensitivity exploration of 
the variable radio sky, and (d) wide-field galactic and 
extragalactic astrophysics including a deep all-sky sur- 
vey over the MWA frequency range with full polarime- 
try and spectral resolution. An initial survey of 2400 
square degrees on the prototype '32T' system has re- 



cently been completed ( Bemardi et al. 2012 1 and a full 
sky survey south of +30° is being planned for the full 
array. 

2.11 Other survey telescopes 

We acknowledge that many other radio telescopes, such 
as the Giant Metre Wave Radio Telescope (GMRT: 
Ananthakrishnan & Pramesh Rao||2001 1 and the Very 



Long BaseUne Array (VLBA: [Napier et al.|1994[ ) are 

also making very significant contributions both to the 
continuum survey science and to the survey techniques 



necessary for the SKA, but are outside the scope of this 
paper 

3 Science Goals 
3.1 The Radio Sky 

The goal of most SKA pathfinder radio continuum sur- 
veys is to make a deep radio continuum survey of a 
significant fraction of the sky. For example, EMU and 
WODAN together wiU cover the entire sky at 1 .4 GHz 
to a depth of about 10-15 /iJy/beam rms with an an- 
gular resolution of 10-15 arcsec. LOFAR will com- 
plement this with a low-frequency survey of compara- 
ble intrinsic sensitivity, and VLA and Meerkat will be 
able to probe individual objects, or survey smaller ar- 
eas, to a much greater depth. The largest comparable 
survey currently existing is the benchmark NVSS sur- 
vey ( Condon et al.[[T998) l, compared to which EMU 
plus WODAN will be about 40 times more sensitive, 
with an angular resolution about four times better, and 
with higher sensitivity to extended structures. EMU 
and WODAN together are expected to detect and cata- 
logue about 100 million sources. 



< Windhorst 1993 compilation 
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Figure 3 The Euclidean normalised differential radio 
source counts at 1 .4 GHz, based on and updated from 



the distribution shown in Hopkins et al. 


(2003|l. The 


solid curve is the polynomial fit from 


lopkins et al. 



(j2003|, and the dashed curve is an updated polyno 
mial fit. The horizontal dot-dashed line represents a 
non-evolving population in a Euclidean universe. The 
shaded region shows the prediction based on fluctua- 
tions due to weak confusing sources ( a "P(D) anal- 
ysis") from [Condon[ ( [1974[ l and [Mitchell & Condon[ 
( fT985] i. 



At the sensitivity level of these surveys, the radio 
sky has already been studied by a number of deep sur- 
veys, and the flux-density-dependent surface-density 
statistics have been extensively measured (Figure |3]i. 
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However, these surveys cover a total of only a few 
square degrees jHopkins et al. 2003 ; Huynh et al. 2 005 [ 



Biggs & Ivison 2006„Norris et al. 2006; Middelbe rget] 
al.||2008at [Miller et al.||2008r[ Schinnerer et al. 2^7) 
Mon-ison et al."2010; 'Hales et al. 2012b; BanfielJetj 
al.(2012, Guglielmino et al. ,2012j . We know broadly 
what types of galaxy we will detect, and their approx- 
imate redshift distribution (see Figure [4]i but a number 
of significant questions remain: 

• What types of object dominate the source counts 
at low fluxes? 

• Are there undiscovered populations? 

• What causes the scatter between different sur- 
veys at low fluxes? 

• Do we need better simulations to guide research 
in this area? 

• How do environment and large-scale structure 
affect the evolution of radio sources? 
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Figure 4 Expected redshift distribution of sources with 
SiA > 10 /zJy/beam, based on the SKADS simu- 
lations ( IWilman et aLl|2008l [20T0l l. The five lines 
show the distributions for star-forming galaxies (SFG), 
starburst galaxies(SB), radio-quiet quasars (RQQ), and 
radio-loud galaxies of Fanaroff-Riley types 1 and II 
(FRI & FR2; [Fanaroff & Riieyl[T97?l l. The vertical 
scale shows the total number of sources expected to be 
detected. 



At high flux densities, the source counts (Figure [3]l 
are dominated by AGN, following a smooth power law 
distribution from Jansky levels down to about 1 mJy. 
Below 1 mJy, the normalised source counts flatten, sug- 
gesting an additional population of radio sources. The 
source counts below 1 mJy are less well-constrained, 
with a scatter between different surveys which can reach 
a factor of ^ 2 — 3. 

While some fraction of this scatter may be due to 
cosmic variance (e.g. |Prandoni et al.|20T2] l, that can- 
not be the only cause. For comparison, in the X-ray 



band, the scatter between the GOODS Nort h and South 
fields is only ~ 25% (~ 3a). Furthermore, Condon et 
al. (2012 1 have shown that some of the largest varia- 
tions can be accounted for by errors in the many nec- 
essary corrections for effects such as resolution bias, 
clean bias, Eddington bias, and completeness correc- 
tions, together with imaging errors such as excessive 
deconvolution, bandwidth smearing, and insufficient 
beam sampling in the image plane. There is also an in- 
consistency in that some authors measure the numbers 
of radio components, while other measure the numbers 
of sources, each of which may consist of several com- 
ponents, the numbers of which in turn may vary as a 
function of flux density. 

Resolving this discrepancy is critical, since this fac- 
tor of 2 — 3 scatter introduces large uncertainties in 
the comparison of observed number counts with de- 
tailed, model-based predictions. Obviously, radio sur- 
veys need to move beyond the small-area surveys that 
have previously dominated the literature. It is also es- 
sential that, before the large radio surveys start, we ob- 
tain consensus on how to image and analyse radio sur- 
vey data to produce reliable and consistent source flux 
densities and source counts. 



3.1.1 The radio sky at ^Jy levels: AGN or star for- 
mation? 

After years of intense debate, it is now clear that while 
star-forming (SF) galaxies contribute significantly to 
the sub-mJy counts, they do not dominate it above ^ 
200 iiJy. Recent re sults (e.g. jS eymour et al. 2008| 
[Smolcid et aL][2008l [Padovarii et"al.||2009| ), which are 
summarised in Figure |5] show that both SF galaxies 
and AGN contribute significantly to the source counts. 
About half of these AGN are radio-quiet, characterised 
by relatively low radio-to-optical flux density ratios and 
radio powers, as compared to radio quasars (see also 
Prandoni 20101. Above 1 mJy, the contribution of radio- 



quiet objects to the source counts is insignificant, and 
they are dominated by radio-loud AGN. 

To probe the radio-quiet AGN component, it is there- 
fore necessary to analyse the deepest radio fields (S < 
100 — 200 fiiy), where most sources are typically as- 
sociated with star-forming galaxies. Radio-quiet AGN 
share many properties with star-forming galaxies: they 
have similar radio luminosities (10^^^^* W/Hz), steep 
radio spectra (a < —0.5), and similar infrared/radio 
flux ratios (e.g. |Roy et al.|1998[ ). 

In addition radio-quiet AGN are often associated 
with faint optical galaxies characterised by Seyfert-2- 
like spectra, containing both star-forming and AGN 
components. Even when available, if optical spectroscopy 
is not of sufficient quality, it can be difficult to derive 
reliable emission line ratios, and distinguish between 



starburst and Seyfert 2 spectra, through BPT (Baldwin 
|et al.|198T] ) diagrams, (see e.g. Prandoni et al.|2009] r 
This makes it difficult to distinguish radio-quiet AGN 
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from star- forming galaxies, even with multi wavelength 
information, although the availability of X-ray data can 
help (| Padovani et al.|201 1 



Spitzer IRAC colours, when available, can be very 
effective in separating radio quiet AGN from star-forming 
galaxies as shown in Fig ure [6} where the well-known 
IRAC colour-colour plot ( Lacy et al. 2004| i is exploited 
for the First Look Survey (FLS) deep radio field. A 
fraction of sources with typical AGN IRAC colours ap- 
pear to be genuine radio-quiet AGN, and account for 
^ 45% of the overall AGN component in the FLS. 




13"XCDF - Seymour et al, SOOB. 
HDFN - Uullow et &1. 3005 
SXDF - Simpson ei al. 200e 
COSMOS - Smoicic et al SOOfl . 
Lockman - Ibar et al. 2009 
ECDFS - Padovani et al. 2Q09 . 



0.1 1 
flux density (mJy) 

Figure 5 Differential fraction of star-forming galaxies 
as a function of 1.4 GHz flux density, taken from |Nor^ 
|ris et al.| ( |201 Ih) . Shaded boxes, and the two lines for 
Padovani et al., show the range of uncertainty in the 
survey results. Arrows indicate constraints from other 
surveys. These results show that the fraction of star- 
forming galaxies increases rapidly below 1 mJy and, at 
the 50 ^Jy survey limit of EMUAVODAN, about 75% 
of sources will be star-forming galaxies. 



A number of other methods can also be used to dis- 
tinguish AGN from SF galaxies in radio surveys (e.g. 
Norris et al.|2008| l, including 



Radio morphology (e.g. Big gs & Ivison|[2008 



Biggs, Younger & lvison^20l0 1 



• Radio spectral index (e.g. Ibar et al. 2009 20101, 

• Radio-far- infrared ratio (e.g. |Norris et al.|2006t 



Middelberg et al.|2008al l. 



Radio-near-infrared ratio (e.g. |Willott et al.|2003 1 
Radio polarisation (e.g. Hales et al.|2012b 
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Figure 6 IRAC colour-colour plot of the FLS radio 
sources with (724 = log{S24^/ Siaghz > 0.6 (fifled 
symbols), from Prandoni et al.| ( |2009| l. Colours refer 



to optical spectral classification: star-forming galax- 
ies (blue); narrow/broad line AGN (green); early-type 
galaxies (red); galaxies with narrow emission line, 
which do not have a secure optical classification (ma- 
genta); sources with no optical spectroscopy available 
(black). Arrows indicate upper/lower limits. The lines 
indicate the expected IRAC colours as a function of 
redshift for different source types (see legend). The 
expected location for AGN is highlighted in pink. For 
reference we also show IRAC colours of: (a) all FLS 
IRAC-identified radio sources (no optical identifica- 
tion selection applied, cyan dots); (b) the entire FLS 
IR-selected star/galaxy population (no radio selection 
applied, black dots); and (c) a sample of high redshift 
obscured (type-2) quasars (Martinez-Sansigre et al 
j2006j), green crosses. 



• X-ray hardness ratio (e.g. Rosati et al.|2002 1 

• X-ray power (e.g. [Padovani et al.|201 l| l, 

• Radio source brightness measured by VLBI (e.g. 



Kewley et aL]|2000| [Parra et al.||2()T()| [Alexan- 
droff et al.|2012l L 



None of these techniques is foolproof and univer- 
sal, and a combination of techniques is necessary to 
provide unambiguous classification. In preparation for 
the daunting task of classifying tens of millions of ra- 
dio sources, we need to define reliable methods for 
classifying radio sources. 



Radio variability (e.g. [Norris et al.|2008tpiat^ 3.1.2 Composite Galaxies 



terjee et al.|2"0T0[ [Murphy et al.|2012[ i. 



• Optical and IR colours, including the use of spec- 
tral energy distribution (SED) templates (e.g. 
Lacy et al.l2004) . 



It is well-estabhshed (e.g. [Roy et al.|1998[ ) tiiat the ra- 
dio emission from Seyfert galaxies contains significant 
contributions from both AGN and star forming activity, 
and there is growing recognition ([Norris et al.[[2008[ 



Optical Une ratios (e.g. Baldwin et al.|198l] l, 



[Lutz et al.|20T0l[Shao et al. [20 10[ [Seymour et al.|201 l| l 
that this may also be true of high-luminosity galaxies, 
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particularly at high redshift. Many of these galaxies 
are not simply "star-forming" or "AGN" but include a 
significant contribution from both. In extreme cases 
( Norris et al.||2012 i a galaxy may appear to be "pure 
AGN" at one wavelength and "pure SF" at another. 

Particularly at high redshifts, such composite AGN/SF 
systems constitute a significant fraction of radio sources, 
and so a simple classification into AGN or SF galaxy 
is inadequate. Instead the contribution to the galaxy's 
luminosity from both AGN and SF activity must be 
assessed, to measure the relative contributions from 
underlying physical properties such as black hole and 
galaxy mass, star- formation rate, environment, etc. How- 
ever, the relative contribution will vary depending on 
the observing band, and so comprehensive multiwave- 
length data are required. However, such detailed infor- 
mation will be hard to obtain for very large samples, 
perhaps necessitating a simpler classification scheme, 
which may lead to oversimplifications or even to incor- 
rect interpretation of the data. 

It is therefore essential to develop methods which 
fit the SED of a given radio source with both a star- 
burst and AGN component (e.g. Afonso et al.|200l"| l. 
This approach is clearly only possible for fields with 
good multiwavelength data, but will inform analysis of 
fields with lower quaUty ancillary data. Observations 
from Spitzer and Herschel cover the bulk of the energy 
output of most galaxies and can be used to distinguish 
the starburst and AGN components. For example, the 
far-IR is directly related to the star formation rate, and, 
if the FlR-radio correlation remains true at all redshifts 
(e.g. |Mao et al.|2(3lla| and references therein), can be 
used to determine what fraction of the radio emission 
is due to star formation. Hence, the remaining radio 
emission would be due to AGN activity. The AGN 
bolometric luminosity can be estimated from the mid- 
IR component which is not necessarily proportional 
to the AGN emission at radio wavelengths. Once the 
starburst and AGN fractional contribution to individual 
sources are determined, it is possible to construct de- 
rived results such as luminosity functions and relative 
contributions by type to the source counts. 

3.1.3 Are there undiscovered populations of radio 
sources? 

The cosmic radio background (CRB) has been well- 
studied at high frequencies by instruments such as COBE 
and WMAP, but is less well-studied at low frequencies. 
Fixsen et al. ( |201 l| l measured the background sky 



temperature at five frequencies between 3 and 90 GHz 
using an in situ calibrator on board the balloon-borne 
ARCADE2 experiment. These authors found an ex- 
cess radio sky temperature above that due to the Cos- 
mic Microwave Background (CMB) below 10 GHz. The 
CRB was measured to have an excess of 50 ± 7 mK at 
3.3 GHz above a CMB temperature of 2.730±0.004 K, 
which is a factor of five higher than expected from 



known source populations. 

The ARCADE2 result is either due to an instru- 
mental or calibration problem (e.g. incorrectly sub- 
tracting the Galactic foreground emission), or it is a 
startling result which will necessitate a drastic revision 
of our models of extragalactic radio sources. Given 



that Fixsen et al. (2011 1 appear to have taken careful 
steps to avoid errors, there is a prima facie case that 
this result is correct, and the radio-astronomical com- 
munity has been galvanised to search for this putative 
new population. Since it is inconsistent with known 
radio source populations (e.g. |De Zotti et al. |2010 



Padovani 2011 1, it must, if confirmed, be caused by 



another population. For example, the radio emission 



could be caused by dark matter annihilation ( Fornengo 



|et al.| 20lT), in which case the emission would trace the 
dark matter distribution of cluster galaxies, resulting in 
a scale size of ^ arcmin. Other mechanisms are also 
possible, such as diffuse emission from clusters or ha- 
los, or very low surface brightness emission from ex- 
tended AGN radio lobes, or from a population of dwarf 
galaxies. 



Condon et al. 



(2012 1 have conducted a deep sur- 
vey with the VLA which shows that the ARCADE2 
result cannot be caused by a population of compact 
objects spatially associated with galaxies. It is there- 
fore either due to a diffuse population of objects which 
cannot be detected with the shortest baselines of the 
VLA observations, or they are due to an instrumental 
or calibration error. A number of observational pro- 
grams continue with the VLA and ATCA to determine 
the cause of this result. 

3.1.4 Simulations 

There is widespread agreement that the SKADS con- 



tinuum sky ( Wilman et al. 2008 1 represents the current 
state of the art in simulations of the radio sky. Such 
simulations have two distinct functions. 

• They represent useful approximations to the real 
sky which can be used to test algorithms and de- 
sign surveys. 

• They can be viewed as an attempt to encapsulate 
the current "best knowledge" of the sky, and can 
therefore be used to test models of radio source 
evolution. 

These two functions have different implications for 
future updates of the SKADS continuum sky. For ex- 
ample, the SKADS continuum sky was constructed from 
the best information then available of the radio lumi- 
nosity functions of different populations of radio ob- 
jects, and deeper and wider surveys have since become 
available. The SKADS simulation also assumes that 
sources are members of distinct populations, omitting 
the composite sources in which both star formation and 
AGN activity contribute to their radio luminosity. Such 
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small differences from the real sky are unlikely to be 
important for the design of a survey, but will be impor- 
tant if SKADS is to be used to test theoretical models. 

A useful, and relatively simple, modification to the 
SKADS continuum sky would be to change the as- 
sumed luminosity function for one or more of the pop- 
ulations in order to assess how that might influence 
the results of various continuum surveys. This change 
might be implemented in a parametric way, allowing 
users to specify parameters of a luminosity function. 

More difficult would be an attempt to incorporate 
"best knowledge," with the aim of producing "realis- 
tic" sky images. Experience from the Herschel mission 
bolsters the concern that identifying and incorporating 
information can be extremely difficult. Similarly, in- 
corporating polarisation information would be useful, 
so that one could calculate Faraday rotation measures 
(RMs), but would be quite difficult in practice because 
of the need to trace the geodesic of a photon in or- 
der to determine the RM. Nevertheless, there are ef- 
forts underway to insert realistic source shapes into the 
SKADS sky, and the SKADS team is investigating how 
to add the signature of the cosmic dipole. 

It would be very useful to include models of more 
complex structures in the SKADS sky, such as those of 
large radio sources, diffuse cluster sources, and struc- 
tures near the Galactic plane, which could potentially 
limit the lowest Galactic latitude to which an "all-sky" 
survey could probe. However, as described in Section 



4.3 imaging a field containing many compact sources 



may have similar issues as to imaging a field contain- 
ing large scale structure. 

3.2 The evolution of star formation 

The cosmic star formation history has been studied thor- 
oughly over the past decade and a half, building up a 
remarkably consistent picture. As shown in Figure |5] 
only a small fraction of radio sources above 1 mJy are 
star-forming galaxies. Below 1 mJy, a much larger vol- 
ume of the star forming galaxy population becomes de- 
tectable, and the population becomes significant at flux 
densities below about 0.2 mJy. Star-forming galaxies 
will therefore dominate the deep and wide area contin- 
uum surveys proposed with the ASKAP, APERTIF and 
MeerKAT telescopes, probing the entire sky to depths 
of about 50 /^Jy/beam. 

The star formation history is now well established 
up to almost the epoch of reionisation. The space den- 
sity of star formation rate declines by an order of mag- 
nitude between a redshift of unity and zero, and is al- 
most flat at higher redshifts, with the suggestion of a 
decline above a redshift of z ^ 5 (Hopkins & Bea com| 
|2006| ). At higher redshifts there is still some discrep- 
ancy, with evidence from Lyman dropout sources sug- 
gesting a dramatic decline above z=6 (|Bouwens et al 



|2008[ ), and evidence from gamma-ray bursts suggest- 
ing a much slower decline ( Kistler et al.|,2009^ . The 



mass-dependence of the star formation rate density has 
also begun to be constrained, with the contribution of 
more massive galaxies to the star formation rate den- 



sity increasing at higher redshift (e.g. Juneau et al. 



[20051 iFeuiiier et al.|200 5 ; Mobas her et al.|2009| l. En 
vironmental effects as well have been shown to change 
with redshift, with star formation in high density, clus- 
ter like environments being enhanced at early times 



(Elbaz et al. 2007 1 but suppressed at later times (jLewis 
let al.|2002HGomez et al.|2003] l. 



Together these lines of evidence point to a picture 
where star formation at the earliest times is dominated 
by massive systems in the densest regions, which rapidly 
exhaust their fuel, ultimately becoming the progeni- 
tors of today's brightest cluster galaxies and cD galax- 
ies. At later epochs the star formation is dominated by 
somewhat lower-mass systems in somewhat less dense 
environments, similar to cluster outskirts, which in turn 
exhaust their fuel, until at the lowest redshifts the star 
formation is dominated by low-mass galaxies living in 
the least dense environments (Figure |7]i. This picture 
is still simplistic, and omits significant details such as 
how the evolving luminosity function of star forming 
galaxies changes with galaxy mass and environment, 
and whether and how the stellar initial mass function 
varies (e.g. Gunawardhana et al.|201 1 1. 



3.2.1 Measuring star formation rates from radio 
observations 

At redshifts between 1 ^ z 3, when galaxies were 
undergoing rapid evolution, the star formation activ- 
ity appears to have been dominated by dusty, heavily 



obscured, starbursting galaxies (e.g. Chary & Elbaz 



Tool ; Caputi et al. 2007; M agnelli et al. 2010; Mur- 
,phy et al.. 2011j . Optical and near-infrared observa- 
tions are therefore seriously hampered by dust extinc- 
tion. Thus, deep radio continuum surveys provide an 
important tool for measuring the cosmic star formation 
history of the Universe. 

A key result from the next generation of wide and 
deep radio surveys will therefore be the measurement 
of the radio luminosity due to star formation in a wide 
range of galaxies. Converting this radio luminosity to 
a star formation rate (SFR) depends on a conversion 
factor which is principally based on the infrared/radio 
correlation ( |Yun et al. |[200T] |Beir||2003| |Seymour| 
et al.|2008) l. However, this correlation has not yet been 



well determined at the high end of the luminosity func- 
tion, since the number of ultra-luminous IR galaxies 
(ULlRGs) in these papers is very small, or at high red- 
shifts, although evidence so far suggests the correlation 
remains constant ( |Mao et al.|201 la| l. 

At high redshifts, however, we know that the co- 
moving star formation rate density is increasingly dom- 
inated by more luminous galaxies (with naturally higher 
SFRs). In Figure[8]we show the comoving infrared lu- 
minosity density, a proxy for the comoving SFR den- 
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z>3-6 



, ' 1<Z<3; 



z<0.5 



Figure 7 An illustration of the density-dependence of 
the downsizing of galaxy star formation rates, using 
optical images from a variety of sources. At the highest 
redshifts (z>3), star formation is occurring predom- 
inantly in massive galaxies (those that become local 
massive ellipticals) that live in the most overdense re- 
gions (that evolve into today's massive clusters). At 
lower redshifts (1 < z < 3), where EMUAVODAN 
will be sensitive to the most extreme star forming 
systems, star formation is dominated by lower mass 
systems, in less dense environments. By the current 
epoch, star formation is limited primarily to low-mass 
galaxies in the outskirts of clusters, and in the lowest- 
density environments. 



sity, as a function of redshift. We also separate the con- 
tribution by infrared luminosity and find that ULIRGs 
represent an increasing fraction of the total star for- 
mation budget above a redshift of unity. Hence, if we 
are to use deep radio surveys to probe the star forma- 
tion history of the Universe at high redshift, we must 
obtain a considerably more accurate conversion of the 
radio luminosity to SFR for galaxies with high SFRs. 
This conversion factor could also depend upon other 
parameters such as stellar mass, environment, metal- 
licity etc. 

3.2.2 Measuring SFRs at 10 GHz 

Radio continuum emission from galaxies typically arises 
from two processes that are both tied to the massive 
star formation rate (SFR). At low frequencies (e.g., 
^ 2 GHz), the radio continuum is dominated by non- 
thermal synchrotron emission arising from cosmic-ray 
(CR) electrons that have been accelerated by shocks 
from supernova remnants and are propagating through 
a galaxy's magnetised interstellar medium (ISM). This 
physical hnk to massive star formation provides the 
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Figure 8 The comoving infrared luminosity density of 
star forming galaxies as a function of redshift, sepa- 
rated into four infrared luminosity ranges. The shaded 
regions are from jLe Floc'h et al. ( |2005) and the 
bold points are derived from the analysis of faint ra- 
dio sources in the 13H field (Seymour e t al.|[2008| l 
where the star-forming radio luminosities have been 
converted to infrared luminosities using the relation of 
|Bell I ( |2003[). The faint, grey poi nts are from the com- 
pilation of Hopkins & Beacom] ( |2006) l converted to 
IR luminosity density. The figure shows that ULIRGs 
make an increasing contribution to the total star forma- 
tion budget above redshifts of unity. 



foundation for the far-infrared (FIR)-radio correlation 
(e.g. Helou et al.|1985[|Condon|1992[|Yun et al. |2001 



Murphy et al. 2006|. However, this link is not at all 



direct given that there are a large number of physical 
processes affecting the propagation of CR electrons 
and the heating of dust (e.g. CR diffusion, magnetic 
field strength/structure, dust grain sizes and composi- 
tion) that must conspire together to keep this relation 
intact ( |Bell |2003| l. 

High frequency (^ 10 — 100 GHz) radio emission, 
on the other hand, offers a relatively clean way to quan- 
tify the current star formation activity in galaxies. At 
these frequencies, emission is generally optically thin 
and dominated by free-free radiation, which is directly 
proportional to the ionising photon rate of young, mas- 
sive stars. While this picture could be complicated by 



the presence of anomalous dust emission (e.g. 


Kogut 


etal. 


1996||de OUveira-Costa et al. |1997||Leitc 


1 et al. 


1997 


1, which occurs at these frequencies and is thought 



to arise from spinning dust grains (e.g. Erickson 1957 



[Draine & Lazarian 1998) , it is currently unclear whether 
this component contributes significantly to globally in- 
tegrated measurements ( Murphy et al.|[201 0). Thus, 
higher frequency radio observations may be particu- 
larly powerful for precisely measuring the star forma- 
tion history of the Universe unbiased by dust. 

For example, surveys at > 10 GHz start to probe 
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factor of ^ 2 deeper than a corresponding 1.4 GHz 
survey to detect the same population of star forming 
systems. For an intrinsic magnetic field of 10 /iG, this 
small discrepancy between point source sensitivity re- 
quirements for 10 GHz and 1.4 GHz surveys occurs at 
a z > 2. Given that high-frequency observations are 
almost as sensitive to high-z star-forming galaxies as 
low-frequency observations, and that higher frequency 
observations provide a much more robust measure of 
star formation activity, surveys at > 10 GHz will be 
invaluable for accurately tracing the star formation his- 
tory of the Universe. 



Figure 9 The fraction of radio emission of SFG due 
to thermal processes as a function of redshift. Galaxy 
magnetic fields of 10 (solid lines), 50 (dashed lines), 
and 100 /iG (dotted lines) are shown. The radio emis- 
sion from observations at 10 GHz is dominated by free- 
free processes beyond a redshift of z > 2 even for 
magnetic field strengths of ^100 /iG, making it an 
ideal measure for the current SFR of high-z galaxies 
( |Murphy,2009) . 



the rest-frame > 30 GHz emission from star-forming 
galaxies by z > 2 for which >50 % of the emission 
is thermal (free-free). This is illustrated in Figure |9] 
where the thermal fraction of the observed 10 GHz and 
1.4 GHz emission is shown against redshift for star- 
forming galaxies having intrinsic magnetic field strengths 
of 10, 50, and 100 [iG while still obeying the FIR- 
radio correlation at z = 0. In addition to measuring 
higher rest-frame frequencies with increasing redshift, 
this calculation also accounts for the suppression of a 
galaxy's non-thermal emission due to rapid cooling of 
CR electrons from inverse Compton (IC) scattering off 
the cosmic microwave background (CMB), whose ra- 
diation field ene rgy densit y scales roughly as J/cmb ~ 
(1 + z)4 (Mur phyl|2009| . By z > 4, nearly 80% of 
observed 10 GHz radio continuum is due to free-free 
emission for a galaxy having a large (i.e., 50 /fG) mag- 
netic field. 

Additionally, we can compare the relative sensi- 
tivity requirements of surveys at 10 GHz and 1.4 GHz 
to detect cosmologically important galaxies at early 
epochs. 



This is shown in Figure 10 where we show the ex- 
pected 10 GHz and 1.4 GHz flux densities from galax- 
ies having a range in IR (8—1000 [ivs\) luminosities and 
an intrinsic magnetic field of 50 /iG. This Figure essen- 
tially shows the selection function of radio surveys to 
star forming galaxy populations as a function of red- 
shift. While higher frequency surveys need to be much 
more sensitive to detect the same luminosity class of 
galaxies at 1.4 GHz, this becomes less of a problem 
at higher redshifts. For example, by z > 4, the sen- 
sitivity requirement of a 10 GHz survey is less than a 
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Figure 10 The expected 10 GHz (solid line) and 
1.4 GHz (dotted line) flux densities for galaxies of dif- 
ferent IR luminosities assuming an intrinsic magnetic 
field strength of 50 /iG. The depths of possible future 
surveys taken by the VLA, MeerKAT (MIGHTEE), 
and SKA are shown. Since the non-thermal emission 
from star-forming galaxies should be suppressed by in- 
creased IC scattering of CR electrons off of the CMB, 
the discrepancy between the point source sensitivity re- 
quirements of surveys at 1.4 GHz and 10 GHz falls be- 
low a factor of ^ 2 by z > 4. 



Figure 10 also shows the expected depths from po- 
tential continuum surveys using the VLA, MeerKAT, 
and SKA. The tentative Tier 4 and 5 steps of the MIGH- 
TEE survey plan to map 0.25 and 0.01 square degrees 
down to an RMS of 1 and 0.2 /iJy at 12 GHz, respec- 
tively. Although the precise area on the sky has not 
yet been defined, obvious targets are the Chandra Deep 
Field South and COSMOS fields. While the depth of a 
megasecond exposure is shown for the VLA, the area 
of such a survey would be restricted by the small pri- 
mary beam at 10 GHz (6'pb ~ 4 — h arcmin). How- 
ever, VLA surveys reaching /iJy depths for entire fields 
such as GOODS-N at 10 GHz have already been pro- 
posed, allowing synergy with future ALMA observa- 
tions. Clearly, the SKA will revolutionize any such 
high frequency continuum surveys using either the VLA 
or MeerKAT by being >2 orders of magnitude deeper 
for the same integration times. 
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3.3 Evolution of AGNs 

3.3.1 The AGN component in Deep Radio Fields 

Radio source counts at /iJy levels are dominated by 
star-forming galaxies, but their contribution decreases 



AGN Luminosity Function 



rapidly with increasing flux density (e.g. Seymour et 



aL]|2008 1, and AGNs contribute significantly at radio 
fluxes below 1 mJy (e.g. [G ruppioni et al.|1999 



gakakis et aL]| 19991 IMagliocchetti et al.||2000 



eor- 



doni et al.|200Tp .fonso et al. 2006 ; No rris et al.|2006| 



Pran- 



Mignano et al. 2008 ; Smolcic et al. |2008| [Padovani et 
al.|2009 1. Recent estimates of their relative contribu- 
tions are summarised in Figure |5] 



Models of the sub-mJy radio population (e.g. Wilm in 



et al.|2008 1 therefore include three main components: 



• Star-forming galaxies, 

• The extrapolation to low flux densities of the clas- 
sical radio-loud AGN population (radio galaxies 
and radio-QSO), 

• A radio-quiet (or low-luminosity) AGN compo- 
nent. 

The unexpected presence of large numbers of AGN- 
type sources at sub-mJy levels implies a large popula- 
tion of low/intermediate power AGNs, which may have 
important implications for the black-hole-accretion his- 
tory of the Universe. It will be particularly important 
to determine the accretion mode (see § |3.3.2| i of these 
low-power AGNs. 

Radio and optical studies of the Australia Telescope 
ESQ Slice Project (ATESP) sub-mJy sample (|Pran 



|doni et al.|2006||Mignano et al. 2008 ) have shown that 
at fluxes 0.4 < S' < 1 mJy a prominent class of low- 
luminosity (P < 10'^^ W/Hz) AGN associated with 
early-type galaxies appears. These sources have typ- 
ically flat or inverted radio spectra (a > —0.5) and 
compact sizes {d < 10 — 30 kpc). This source class 
produces a flattening of the average spectral index at 
flux densities around 0.5-1 mJy. A similar flattening 
has been found also in deep low-frequency (350 MHz) 
surveys ( |Owen et al.|2009[|Pra"ndoni et al.|201 l|l. Ra- 
dio spectral index studies (e.g. Prandoni et al.||20T0 i 
suggest that this class of objects is presumably asso- 
ciated with core-dominated self-absorbed FRl-like ra- 
dio galaxies, perhaps representing a low luminosity ex- 
tension of the classical radio-loud radio galaxy pop- 
ulation, which dominates at flux densities S >> 1 
mJy. On the other hand, deeper radio fields have shown 
that flat/inverted radio spectra become less frequent at 
lower flux densities (100 < S < 400/xJy) ( |Owen et 
|al.^2009; Prandoni et al. 2009), where steep-spectrum 
star-forming galaxies and radio-quiet AGN are increas- 
ingly important. 

Our current understanding of the AGN component 
in deep radio fields can be summarised as follows: 
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Figure 1 1 AGN luminosity function at different red- 
shifts expected from the combination of ASKAP and 
Meerkat deep surveys (based on models described by 
[Prandoni et aT] ( [20d7| ). 



• above flux densities S ^ 400 /iJy there is no 
clear evidence for a radio-quiet AGN compo- 
nent; 

• sub-mJy radio-loud AGN are mostly low-power, 
compact, (self-absorbed) jet-dominated systems, 
and seem to be a low power counterpart of the 
classical bright radio galaxy population; 

• at lower flux densities radio-quiet AGN become 
increasingly important and at S* ^ 100 /iJy they 
account for about 45% of the overall AGN com- 
ponent; 

• radio-quiet AGN show radio/optical/lR proper- 
ties consistent with radio follow-up of optically 
selected radio-quiet AGN (Kukula et al. 1998|, 
compact radio sizes, steep radio spectra, lumi- 
nosities of P < 10^4 W Hz-\ and Seyfert-2 
optical spectra and/or mid-lR colours. 

Multiwavelength analyses of larger and deeper ra- 
dio samples are needed to confirm these results, and to 
obtain quantitative constraints for evolutionary mod- 
elling of faint AGN. We can expect this field to move 
forward significantly in the next years, thanks to the 
combination of wide-field and all-sky surveys (sam- 
pling large local volumes) and very deep fields (sam- 
pling low powers at high redshifts) planned with the fa- 
cilities described in this paper An example of what can 



be obtained is illustrated in Figure 1 1 where we show 
the AGN luminosity function at different redshifts that 
can be derived from the combination of EMU and MIGH- 
TEE. 
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3.3.2 Accretion Modes 

AGN activity, and associated black-hole growth, oc- 
curs in at least two different modes, each of which may 
have an associated feedback effect upon the AGN host 
galaxy: 

(i) a fast accretion mode (often referred to as cold- 
mode accretion or quasar-mode accretion) asso- 
ciated with quasars (Silk & Rees 1998 1; this 



radiatively-efficient accretion mode may be im- 
portant in curtailing star formation at high red- 
shifts and setting up the tight relationship be- 
tween black hole and bulge masses observed in 



the nearby Universe (e.g.lMagorrian et al. 1998 1. 



(ii) a radiatively-inefficient slow accretion mode (of- 
ten referred to as hot-mode accretion or radio- 
mode accretion) ( |Croton et al.||2006] [Bower et| 
|al. 2006), the observational manifestation of 
which is low-luminosity radio sources; this mode 
is thought to be responsible for maintaining el- 
liptical galaxies at lower redshifts as "old, red 
and dead" (e.g. Best et al. 2006) and for pre- 
venting strong coohng flows in galaxy clusters 
(e.g. |Fabianetal.|2003) . 



One of the most fundamental issues in understand- 
ing the role of AGN in galaxy formation is the need 
to measure accurately the cosmic evolution of quasar 
activity and the accretion history of the Universe, and 
to compare this with the build-up of the stellar popula- 
tions of galaxies. Do black holes and their host galax- 
ies grow coevally, or does one precede the other? What 
is the primary mode of black hole growth? 

Much of the growth of black holes is believed to 
occur in an obscured phase, and these "Type-2" AGN 
are difficult to identify. Even the deepest current X- 
ray observations do not detect the most heavily ab- 
sorbed sources which Cosmic X-ray Background syn- 
thesis models predict exist in abundance (e.g. GiUi 
et al.pOO? I, implying that much high-z quasar activ- 



ity has yet to be detected directly. These "radio-quiet" 
quasars are not radio-silent, and their radio luminosity 
distribution peaks at about Li.4ghz ~ lO^'^WHz"^ 
(e.g. Cirasuolo et al.|[2003| l. Deep radio surveys with 
the SKA pathfinders therefore offer an alternative route 
to identifying these distant AGN, in a manner unbi- 
ased by dust and gas absorption at other wavelengths. 
Combined with multiwavelength datasets to separate 
source populations, the relative contribution of radio- 
quiet AGN to the faint radio source population can 
be determined. This will enable investigation of the 
dependence of the fraction of obscured AGN on lu- 
minosity and cosmic epoch, and thus the history of 
radiatively-efficient accretion in the Universe to be de- 
termined. 

The deep and wide radio surveys will also allow 
study of the role of low-luminosity radio sources in 



galaxy evolution. The low-luminosity radio- AGN pop- 
ulation is dominated by a population of sources in which 
there is little evidence for radiative emission from an 
accretion disk, and the bulk of the accretion power is 
channelled into the expanding radio jets (e.g. Merlonil 



f 



,& Heinz 2007 ; Hardcastle et al. 2007; Best et al. 2012 
These jets pump energy into their environments, in- 
flating cavities and bubbles in the surrounding inter- 
galactic and intra-cluster medium, from which estimates 
of the mechanical energy associated with the jets can 
be made (e.g. |Cavagnolo et al.||2010[ ). Studies of 
radio- AGN in the nearby Universe have shown that the 
fraction of galaxies that host radio-loud AGN (with 
ii.4GHz > lO^'^ WHz^^) is a strong function of stel- 
lar mass (e.g. Best et al. 2005 | l, and that the time- 
averaged energetic output associated with recurrent ra- 
dio source activity may indeed be sufficient to control 
the rate of growth of massive galaxies. However, sen- 
sitivity limits of current large-area radio surveys mean 
that these low luminosity sources are only observed in 
the nearby Universe. How does the relation between 
galaxy mass and radio-AGN fraction (the radio source 
duty cycle) evolve with redshift out to the peak epoch 
of galaxy formation? Out to what redshift does radio- 
AGN-heating continue to balance cooling? What is the 
differential evolution of the radiatively efficient and in- 
efficient accretion modes (c.f. [Best et al.|2012[ ), and is 
evidence for "down-sizing" in the radio source pop- 
ulation (e.g. Ri gby et al.||201l| ) confirmed? What 
drives these processes? These are all questions for the 
next generation radio arrays, coupled with high-quality 
multiwavelength datasets. Large-area surveys are re- 
quired in order to study a sufficient volume to include 
the full range of galaxy environments, given the very 
important role that large-scale environment can play in 
the evolution of galaxies. 

3.3.3 Ultra-steep spectrum radio sources below the 
mjy barrier 

Another exciting prospect for the upcoming generation 
of radio continuum surveys is that of finding complete 
populations of high-redshift radio galaxies (HzRGs). 
These are among the most luminous galaxies and seem 
to be associated with the most massive systems (e.g. 



van Breugel et al. 1999; Jarvis et al. 2001a; Willott et 
al..2003, Rocca-Volmerange et al. 2004, De Breuck et 



al.|2005i [Seymour et al.||2007| l often going through a 
phase of violent star formation at at the 1000 Mq yr~^ 
level and showing large gas and dust reservoirs (e.g. 
I Dunlop et alTj| 1994; Ivison 1995; Hughes et al. 1997) 
'ivison et al.'1998 ; Papadopoulos et al. 2000; Archibald] 
et al. 2001; Klamer et al. 2005; Reuland et al. 2003| 
[2004, 2007; but see RawHngs et al. 2004). HzRGs are 
considered to be the progenitors of the brightest cluster 
ellipticals, and have been used as beacons to identify 
overdensities in the distant universe, i.e., proto-cluster 
environments at z ~ 2.5 (e.g. Stevens et al.|2003„Vei>] 
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|emans et al.|2007| . Identifying and tracing the evolu- 
tion of HzRGs offers a unique path to study galaxy and 
large-scale structure formation and evolution from the 
earliest epochs, and extensive studies will finally be- 
come possible over the next few years. 

Until recently, searches for HzRGs have been lim- 
ited to wide-area (and thus relatively shallow) surveys. 
This is appropriate as the luminosities of HzRGs are 
large enough to make them detectable to very high- 
redshifts (beyond z ~ 5) at flux densities of tens or 
even hundreds of mJy. For example, TN J0924-2201, 
the most distant known radio galaxy at z = 5.2, ex- 



hibits a 1 .4 GHz flux of over 70 mJy ( De Breuck et al. 
2000 van Breugel et al.|[l999 1. But understanding a 



population, its origin and evolution and how it fits in 
the overall galaxy population, requires more than just 
the detection and study of the most extreme objects. 
The selection of more complete samples, over a wider 
redshift range, requires deeper surveys over wide ar- 
eas. However, the large population of low-z (z ^ 1—2) 
SF galaxies at sub-mJy levels presents a considerable 
difficulty for studies of HzRGs, as efficient ways of 
distinguishing between AGN and SF are then neces- 
sary. 

One of the most successful tracers of HzRGs re- 
lies on the relation between the radio spectral index 
and redshift (e.g. |Tielens et al.|1979{ [Chambers e t al.| 



1996 1. Although an ultra-steep (radio) spectrum (USS; 
a ^ — 1) is not a necessary condition for a high red- 
shifts, and, in fact, the USS selection misses a possi- 



bly large fraction of HzRGs (e.g. Waddington et al. 
[79991 IJarvis eTa l. 2001b 2009; Schmid t et al.|2006| l, 
a higher fraction of high-redshift sources can be found 
among those with the steepest radio spectra. The USS 
criterion is so effective that most of the radio galaxies 
known at z > 3.5 have been found using it (see |Mi-| 
ley & De Breuck||2008 for an overview of the selec 



tion of HzRGs). Surprisingly, there is still no satisfac- 
tory explanation for this spectral-index-redshift cor- 
relation, as neither (a) a combination of an increased 
spectral curvature with redshift and the redshifting of a 
concave radio spectrum (e.g., Krolik & Chen 19911) to 
lower radio frequencies nor (b) radio jets expanding in 
dense environments, more frequently observed in dis- 
tant proto-clusters, which would appear with steeper 



spectral indices ( 


Klamer et al.|2006||Bryant et al.|2009t 


Bomancini et al. 


2010|l seem to explain the observa- 



tions. 

With the next generation of radio continuum sur- 
veys, the search for USS HzRGs will be extended to 
the whole sky and to microJansky flux levels. Although 
such collections have been assembled before (e.g. Bondi 
et al.|2007[[Owen et al. 2009; Afonso et al. 2009; Ibai- 



et al.| 2009) only very recently has it been possible 
to characterise them. Using deep optical and infrared 
imaging of the Lockman Hole, [Afonso et al. ( 201 1\ 
have started a detailed analysis of the /iJy USS radio 
source population, something that wiU be possible over 



a significant fraction of the sky using radio surveys 
such as EMU, WODAN and LOFAR. 

This first detailed analysis of the USS /iJy popu- 
lation used the deep multiwavelength coverage of the 
Lockman Hole. The sample selection was made us- 
ing ultra deep VLA 1 .4 GHz (reaching 6 /iJy rms) and 
GMRT 610 MHz data (15 ^Jy rms), and resulted in 58 
sources with a spectral index a|fo° < —1.3. The 
deep IR coverage at 3.6 /im and 4.5 /im provided by 
the Spitzer Extragalactic Representative Volume Sur- 
vey (SERVS: [Mauduit et al.|2012 ) resulted in an ex- 
tremely high identification rate, above the 80% level. 
Mid-infrared colours are compatible with a significant 
AGN presence among the radio-faint USS population. 
Spectroscopic redshifts for 14 sources and photometric 
redshifts for a further 19 sources were used for the red- 
shift distribution of these sources, ranging from z ^ 
0.1 to z ^ 2.8 and peaking at z ^ 0.6 (see Figure 12 1. 
Twenty-five sources have no redshift estimate, includ- 
ing the faintest sources at infrared wavelengths, which 
indicates that higher redshifts are likely in this popula- 
tion. In this respect they may be similar to the Infrared- 



Faint Radio Sources (IFRS), first identified by Norris 
et al.] ([2006), and which appear to be very high redshift 
dusty radio galaxies (IG am & Alexander^2008,|Huynh 
et al.|2010|pi ddelberg ~t al.|2008b | ' 



j2011b| [Norris et 
|Zinn etal.[20TT] r 



al.[2007[|2()TTat|Cameron et al.|201 1 



A comparison with the SKADS Simulated Skies 



models (Figure 13 1 indicates that FRls and RQ AGNs 
may constitute the bulk of the USS-population at /iJy 
radio flux densities (FRlls are thought to dominate above 
the mJy level). According to the models, the USS 
technique will be as efficient for the selection of very 
high-z radio sources at fxJy fluxes as when applied at 
much higher flux density levels. This raises exciting 
prospects for the next generation of wide-field radio 
surveys, as it increases the potential for detailed stud- 
ies of sources at very high redshifts and understanding 
of the evolution of the AGN activity in the Universe. 

These prospects are strengthened by the recent work 
of [Ker et aL] ( [201 2| who have shown that redshift and 
spectral index are correlated in spectroscopically com- 
plete samples, and argue that this is due both to the 
dense intergalactic medium and the early evolutionary 



stage of the AGN, although they, together with Zinn et 
al. ( 2012a| l, also show that the redshift-spectral -index 



correlation is weaker than other potential redshift indi- 
cators. 



3.3.4 The luminosity-dependent evolution of the 
radio luminosity function 



Section 3.3.2 described how AGN radio jets interact 
with the surrounding intergalactic and intracluster medium 
to prevent both large-scale cluster cooling flows and 



the continued growth of massive elUpticals (e.g. Fabian 
'et al.'2006'; Best et al.|2006l [20071 [Croton et aLpOOgf 
|Bower et al.„2006| ). Determining the evolution of the 



18 



Publications of the Astronomical Society of Australia 



15 



I I I — I — I — I — I — I — I — I I I I — I — I — I — I — I — I — I — I I I — I — I — I — r 



Figure 12 Redshift distribution for radio-faint USS 
sources in the Lockman Hole, from [Afonso et al.| 
( |201 l| l. Filled histogram denotes sources with a spec- 
troscopic redshift determination, while the open region 
refers to photometric redshift estimates. A further 25 
USS sources (43% of the full sample) exist but with no 
redshift estimate, mostly at fainter 3.6 fj,m fluxes and 
likely to be found at higher redshifts. 



Survey Model: 0.6 deg^. S„o>100 /nJy 




Figure 13 Predictions from the SKADS Simulated 
Skies models for the redshift distributions of radio 
source populations, irrespective of their radio spectral 
indices, for a radio survey reaching a detection sensi- 
tivity of 100 /.tJy at 610 MHz over 0.6 square degree, 
similar to the Lockman Hole radio survey considered 
in [Afonso et al.| p011| l. The observed redshift distri- 
bution for USS sources in that work is also displayed. 



radio luminosity function (RLF) is therefore important 
for understanding the timescales on which they impose 
these effects. 

A key early study of radio-loud AGN by [Dunlop 



|& Peacock] ( |1990| l found increases of two to three or- 
ders of magnitude at a redshift of two, compared with 
the local Universe, in the comoving number density of 
both flat and steep spectrum radio-loud sources. They 
also saw indications of a high redshift (z ~ 2.5) de- 
cline in density, but their work, and that of subsequent 
studies in this area (e.g. [Shaver et al.||1996[ [Jarvis et 



[al.[[2001b[ [Waddington et al.[|2001[ nias been limited 
by the inability to probe the depth and volume neces- 
sary to probe the high-redshift behaviour This has mo- 
tivated the development of the Combined EIS-NVSS 
Surve y of Radio Sources (CENSORS: jTSest et al7~[ 



2003 p , which has been designed to maximise the infor- 



mation for high-redshift, steep spectrum sources close 
to the break in the RLF CENSORS is a 1.4 GHz se- 
lected sample, and contains 135 sources complete to 
a flux density of 7.2 mJy. It is currently 78% spectro- 



scopically complete ( Ker et al. 2012 1, with the remain- 
ing redshifts estimated via either the K — z or I ~ z 
magnitude-redshift relations. 

The CENSORS sample, together with additional 
radio data, source counts and local RLF, is used to in- 
vestigate the evolution of the steep spectrum luminos- 
ity function via a new grid-based modelling method in 
which no assumptions are made about the high-redshift 
behaviour ( Rigby et al.|2 01 1). Instead, the RLF varia- 
tions are determined by allowing the space densities at 
various points on a grid of radio luminosities and red- 
shifts to each be free parameters, and then simply find- 
ing the best-fitting values to this many-dimensional 
problem. The modelling finds conclusive evidence, at 
> 3(7 significance, for a luminosity-dependent high- 
redshift decline in space density for the steep-spectrum 
sources. At low radio powers (Pi. 4ghz = 10^^~^^W/Hz) 
the space densities peak at z ^ 1, but move to higher 

redshift for the higher powers (z ^ 3 for Pi.4ghz > 

lO^^W/Hz). This is illustrated in Figure 14 and is 
similar to the "cosmic downsizing" seen for other AGN 



populations (e.g. De Zotti et al. 2010,.Wall[2008][Hasr 



|er et al.|2005[ [Richards et al.|2005] i7 These results are 
robust to the estimated redshift errors and to variation 
in the radio spectral index with redshift. 

The evolution of low-luminosity radio AGN is less 
well understood than that of radio-loud AGN, with some 
studies finding no evidence for any evolution of the 
RLF for low-luminosity radio AGN (e.g. [Clewley[ 
& Jarvis|2004 1 and others finding that low -luminosity 
AGN do evolve with redshift, although more slowly 
(less than a factor of 10 from z = to z = 1.2) 



than their high-luminosity counterparts ( Smolcic et al. 
McAlpine & Jarvis [[20n] i. [Best et al [ ( [2012| i 



2009 



suggest that the luminosity dependence of the evolu- 
tion of the AGN RLF may be attributed to the varying 
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Figure 14 The dependence of the redshift of the peak 
in space density (zpcak) on radio power for the best- 
fitting steep spectrum grid, showing that the numbers 
of the most powerful radio sources peak at the high- 
est redshifts ( ,Rigby et al.||201 1 1. The shaded region 
and error bars give two different determinations of the 
spread in the measurements. 



fractions of hot and cold mode sources with redshift. 



Recent results ([M auch & Sadler|[20 07 ; Sad ler et al.| 
[20071 IPadovani et al.|2011||Mao et al.|2012| ) fail to re- 
solve this discrepancy, with ambiguity added by factors 
such as evolution, cosmic variance, differences in clas- 
sification, differences in terminology, and the difficulty 
of distinguishing low-luminosity AGN (LLAGN) from 
SF galaxies, or the bulge from the disk of a galaxy. 
Next-generation radio surveys using SKA pathfinders 
can resolve this discrepancy provided we develop reli- 
able techniques for distinguishing the AGN component 
from the SF component of a galaxy. 

3.3.5 Environment 



Studies of galaxies in groups and clusters by Best et al. 
( |2007] l; [Wake et al.| ( |2008j have shown that the en- 
vironment of a radio galaxy has a significant effect on 
the radio emission. [Sabater et al.| (2008 ) and |S abater] 
et al.| ( |201 2|) built on this result by comparing (a) 



the radio AGN activity in a sample of isolated galax 
ies (AMIGA sample; Verde s-Montenegro et al.|2005" 



with that of a matched sample of galaxies in clusters 
( [Miller & O wen 2001 ; Reddy & Yun 2004) and, (b) 
the optical AGN activity for the AMIGA sample with 
that of a matched sample of galaxies in compact groups 
( [Martinez et al.|2 0T0). The results of all these studies 
show that radio AGN activity is strongly influenced by 
the environment while this influence is less clear for 
optical AGN activity ( [Miller et al.|2003| l. 

3.3.6 Black Hole Spin 

A further area of study from radio observations is the 
study of black hole physics. For a given accretion rate, 
as traced by the X-ray to mid-infrared SED, AGN show 



a huge range of radio luminosities. This suggests that 
some extra parameter, independent of the accretion rate, 
is controlling the jet power of the AGN. The spin of 
the black hole is a very attractive candidate, for both 



observational and theoretical reasons (e.g. Wilson & 



Colbert] 1995 i. Spin is effectively a hidden parameter. 



which will barely show itself in the SED of the AGN. 
In addition, spinning black holes can convert rotational 
energy into jet power, meaning that jets can be powered 
with nominal efficiencies close to or exceeding unity. 
Observations of the most powerful radio AGN suggests 
the jet power is comparable or larger than the rate of 
accretion of energy (e.g. Fernandes et al.|2()TT I. 

Radio observations allow us to estimate the jet power, 
and combining these with observations at X-ray ener- 
gies or optical wavelengths can allow an estimate of 



the accretion rate and spin. For example, Martmez- 



Sansigre & Rawlings ( 201 l| l consider the jet powers of 
supermassive black holes (SMBHs) with m,>10* M©. 
The authors use a suite of semi-analytic models and 
GR (general relativistic) MHD (magneto-hydrodynamic) 
simulations that predict the jet efficiencies as a func- 
tion of spin, and find that these efficiencies can approx- 
imately explain the radio-loudness of optically-bright 
quasars. 

Martinez-Sansigre & RawUngs ( 2011[ l also find 



that they can explain the local RLE of high- and low- 
excitation galaxies independently (HEGs and LEGs, 
respectively), and that the HEGs and LEGs have differ- 
ent best-fitting spin distributions. The LEGs, modelled 
as ADAFs (advection-dominated accretion flows) or 
low-accretion rate SMBHs, have a bimodal spin distri- 
bution, with the SMBHs having typically low (a ^ 0) 
or high (a ~ 1) spins. The probability density of the 
high-spin peak is significant, typically ^33% of the to- 
tal probability density. On the other hand, the HEGs, 
modelled as SMBHs accreting at a high fraction of 
their Eddington-rate, are found to have typically low 
spins (a ^ only), with a small or non-existent com- 
ponent of high spin. The radio LFs of HEGs and LEGs 
given by these spin distributions can be extrapolated to 
high redshift, where they correctly reproduce the total 
radio LF of all AGN at z = 1. Due to the strong evo- 
lution of the high-Eddington rate SMBHs (as reflected 
by the evolution of the X-ray LF), this modelling also 
predicts that the z — 1 radio sources above a radio 
luminosity density --10^'^ WHz^^ sr^^ (at 151 MHz) 
should show exclusively HEG spectra, while the frac- 
tion of LEGs should increase rapidly at lower radio lu- 
minosity densities. 

Since the space density of HEGs increases rapidly, 
the typical spin of SMBHs will change with redshift. 
At z ~ 0, most of the SMBHs will have the bimodal 
distribution inferred for the LEGs, with a high fraction 
of SMBHs having a high spin (Figure 15 top). At z ^ 
1 the space density of HEG radio sources is larger than 
that of LEG sources, so that the SMBHs will typically 
have the distribution of spins of the HEGs, dominated 
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Figure 15 Inferred spin distributions for SMBHs at 
z = and z = 1, from |Martmez-Sansigre & Rawl- 
ings (2011 1. The vertical dashed lines mark the mean 
value at each redshift. (Bottom) The spin distribution 
at z = 1 is dominated by the low-spin population, with 
a minor component of high-spin SMBHs. (Top) At 
z — a larger fraction of SMBHs possess a high spin, 
so that the mean spin is higher. 



by a low-spin component (Figure[T5] bottom). 

The evolution of the spin from SMBHs (with 
m,>10^ M0) can hence be described as evolving from 
predominantly low spins in high-Eddington rate sys- 
tems at z 1, to a bimodal spin distribution at z ^ 0, 
for mostly low-Eddington rate systems. At the epoch 
of highest accretion, the typical spin is low and at later 
epochs the fraction of SMBHs with high spins increases 
(Figure [T6]l. This can be understood in terms of chaotic 
accretion (e.g. King et al. 20081 spinning the ac- 



creting SMBHs down, whereas major mergers spin the 
SMBHs up (e.g. |Rezzolla et al.|2008 ). At eariier epochs, 
there is a plentiful supply of cold gas, so that accretion 
spins SMBHs down. As this gas disappears, however, 
the braking mechanism that spins SMBHs down is re- 
moved, while some SMBHs will still experience ma- 
jor mergers. Hence a population of high-spin SMBHs 
gradually appears. 

3.4 Galaxy Clusters and Large-Scale Struc- 
ture 

Clusters of galaxies are the most massive systems in 
the present Universe. They are not isolated regions. 
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Figure 16 (T op) Mean spin of SMBHs with m, 



> 



10^ Mq, fromlMartmez-Sansigre & Rawlings (201 1 



There is a gradual increase of the mean spin between 
z ^ 1 and z = 0, which is most likely due to a fraction 
of SMBHs undergoing major mergers and subsequent 
chaotic accretion. (Bottom) Cosmological simulation 
of the mean spin from Fanidakis et al. ( 201 l[ l, which 



includes chaotic accretion and mergers, and which is in 
excellent agreement with the inference based on radio 
observations (top). 



but evolving regions at the intersections of filaments 
and sheets in the large-scale structure. 

Tens of thousands of clusters are currently known. 



but only a few are known at z > 1 ( Wilson et al. 2008 



|Kodama et al.|2007 1, with the highest cluster redshift 
at z = 2.07 ( |Gobat et al.pOTT] ). 

Radio observations of galaxy clusters are power- 
ful probes of the physics of the intra-cluster-medium 
(ICM), but they also provide a complementary way to 
trace the evolution of large-scale structure of the Uni- 
verse up to very large distances. These studies also 
have the potential to explore the way in which galaxy 
evolution (star formation, nuclear activity) depends on 
its environment. 

Radio observations have shown the existence of dif- 
fuse large scale synchrotron radiation from the ICM, 
implying that non-thermal components, magnetic fields 
and relativistic particles, are mixed with the hot ICM. 
The origin of these components is still unclear and a 
subject of hvely debate. Potentially these components 
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contribute to the energy of the ICM and drive com- 
plex physical processes that may significantly alter our 



present (simplified) view of the ICM ( Schekochihin et 



al..2QQ5 i Su bramanian et al.|2006[|Brunetti & Lazarian 
201 ID 



Radio emission from clusters of galaxies is gen- 
erated both by the constituent galaxies and by diffuse 
emission in the ICM. 

The radio emission from constituent galaxies often 
extends well beyond the galaxy optical boundaries, and 
interacts with the ICM. This interaction is observed 



in tailed radio galaxies (see p.4.3 i, and radio sources 
filling X-ray cavities at the centre of cool-core clus- 
ters (McNamara & Nulsen|2007i|Feretti & Giovannini] 



|2008i ). 



Diffuse sources are typically grouped into three classes 
( |Feretti & Giovannini(1996; Kempner et al. 2004j : 

• Halos, which are steep-spectrum Mpc-scale dif- 
fuse sources in the central regions of X-ray lu- 
minous clusters (e.g. 'Ferrari et al. 2008 Venturi 
[20n]|Ferettier al. 2012 ), 

• Relics, which are located in the peripheral re- 
gions of both merging and relaxed clusters and 
caused by cluster-cluster shocks, 

• Mini-halos, which are hosted in relaxed cool- 
core clusters, are centrally located, and usually 
surround a powerful radio galaxy, such as 3C84 
in the Perseus cluster. 

Even amongst these classifications, there is ambi- 
guity, with the term "relic" being used to describe both 
old galaxies and shock accelerated regions of diffuse 
electrons. Furthermore, some amorphous radio objects, 
such as the diffuse "Phoenixes" ( Kempn er et al.|20()4) l 
and "roundish relics" ( |Feretti et al.|2012 1, cannot eas- 
ily be so classified. 

There are two ways in which these classifications 
could be improved. First, we could build in more physics 
to link observations to thermal and non-thermal phys- 
ical processes in galaxy clusters, but must ensure that 
current observational/theoretical expectations do not bias 
our classification. Second, we could use a comparison 
of radio and X-ray surface brightness distributions. 

SKA pathfinder surveys will discover hundreds of 
diffuse radio halos in regions of the sky where there 
are no known clusters ( |Cassano et al.|2012)l an d even 
more tailed radio galaxies ( Mao et al.||2010a Norris 
] |201 1 1. Radio emission in clusters is therefore likely 
to become an important tool both for studying clusters 
themselves, and for detecting large numbers of clus- 
ters to study cosmology and trace large-scale structure 
formation. 

3.4.1 Radio Halos 

About 40 giant radio halos are known so far in the red- 
shift range < z < 0.55. The GMRT Radio Halo 



Survey ( [Venturi et al. |2008] l, in combination with the 
NVSS and WENSS surveys, allows a statistical explo- 
ration of radio halos in galaxy clusters in the redshift 
range < z < 0.4 ( [Cassano et al.|2008| l. The cluster 
radio properties were found to be bimodal: only a frac- 
tion of massive clusters host halos, while no Mpc-scale 
radio emission is detected in the majority of clusters at 
current sensitivity levels. The radio luminosity of these 
"radio quiet" clusters is less than one tenth of that of 
the observed radio halos (Brunetti et al. 2007 2009). 
X-ray follow up of these observations shows that radio 
halos are found only in merging clusters while "radio 
quiet" clusters are systematically more relaxed systems 
( [Cassano et al.|2010a) . 

The radio properties of halos correlate with the ther- 
mal properties of the parent clusters and the radio power 
of halos increases with X-ray luminosity, suggesting a 
direct connection between the thermal and non-thermal 
cluster properties (Feretti et al."1999' 'Bacchi et aLJ 
2003; Cassano et al. 2006, Brunetti et al. 2007, 2009) 
Feretti et al?f20T2| . 



According to a prominent scenario for the origin of 
giant radio halos, a fraction of the gravitational binding- 
energy of Dark Matter halos released during cluster 
mergers can be channelled into turbulence. This turbu- 
lence may in turn accelerate relativistic particles ( [Brun( ;tti 
et al.|200T Petrosian|2001 1 explaining the connection 



observed between radio halos and cluster mergers. 

The source of the synchrotron electrons in these 
halos is unclear, as acceleration of electrons from the 
thermal pool to relativistic energies by MHD turbu- 



lence in the ICM cannot provide sufficient energy ( Pet 



[rosian|200 8l). Instead, there must be a pre-existing pop- 
ulation of relativistic (or supra-thermal) seed particles 
that are then reaccelerated by turbulent acceleration 
during mergers. One possibility is that the seeds are 
secondary electrons generated by collisions between 
cosmic rays and thermal protons ( Brunetti et al.|2005 



Brunetti & Lazarian 201 1) . These hybrid models pre- 
dict that Mpc-scale synchrotron emission from secondary 
electrons ("pure hadronic halos") must also be present 
in relaxed clusters, with a radio luminosity 10-20 times 



smaller than that of merging clusters ( Brunetti & Lazar 
[ian[[201l] l. This emission should be detectable with 
EMUAVODAN. Figure [17] shows the 1.4 GHz Radio 
Luminosity Function (RLF) of radio halos in the lo- 
cal Universe (z = 0.1): at higher radio luminosities, 
turbulent cluster halos dominate, but at lower lumi- 
nosities the RLF is dominated by "hadronic" halos in 
relaxed systems. Additional sources of seed supra- 
thermal electrons include shocks, AGN, galaxies, and 
magnetic reconnection, so relaxed clusters may have 
an even higher radio luminosity. If EMU fails to de- 
tect radio halos in relaxed clusters, this will provide a 
strong constraint on secondary particles and on the en- 
ergy content of cosmic rays in the ICM. 

An extrapolation of the radio power of halos at 1 .4 
GHz versus the cluster X-ray luminosity indicates that 
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clusters with < 10"'^erg/s should host halos of 
power Li.4 < 10^^ W/Hz. With a typical size of 1 Mpc 
their surface brightness could be too low to be detected 
at 1 .4 GHz. However less powerful radio halos should 
be also smaller in size (Ferett i et al.||2012) l implying 
that they should be easily detected by the EMU sur- 
vey. A comparison between 1.4 GHz data with low fre- 
quency data (LOFAR) will be crucial to to test turbu- 
lence properties. A large number of "hadronic halos" 



may be detected at 1.4 GHz (see Figure 17 1, emphasis- 
ing the importance of comparing the EMU and LOFAR 
surveys. 

If turbulence is responsible for the origin of ha- 
los, then the radio-emitting electrons are accelerated 
up to energies rrieC^jm < several GeV, where the syn- 
chrotron and inverse Compton losses quench the accel- 
eration by MHD turbulence. Radio halos will therefore 
steepen at a frequency cx j^b, making it difficult 



to detect halos at higher frequencies ( [Brunetti et al. 
120081). 



Only the few most energetic mergers in the Uni- 
verse can generate giant radio halos with i^s > 1 GHz 
( Cassano et al.|2005] l. Most mergers instead consist of 
major mergers between less massive systems, or mi- 
nor mergers in massive systems ( Cassano et al.|2006 



2010b[ l, which will therefore have steeper spectra (smaller 
i/g)- Recent Monte Carlo calculations based on tur- 
bulent acceleration in galaxy clusters predict that LO- 
FAR, reaching rms= O.lmJy/beam at 120 MHz over 
the northern hemisphere, will detect more than 350 
radio halos at redshift z < 0.6, increasing the num- 



ber of known halos by a factor of nearly ten ( EnBlin 
& Rottgerin"gl|2002[ [Cassano et aL|[2006j . More than 
half of these halos will be "ultra-steep spectrum" ha- 
los, with < 600 MHz ( [Cassano et al. '2010b), pro- 
viding a clear test of the importance of turbulence in 
particle-acceleration in merging clusters. 

3.4.2 Radio Relics 

Relics (occasionally termed radio phoenixes) are irreg- 
ular, elongated diffuse sources which are generally lo- 
cated at the periphery of the cluster, with the Unear 
size of 400 - 1500 kpc ( [Giovannini & Feretti [[2002l l. 
Most have a steep spectrum (a < —1), and are lin- 
early polarised at a level of ~ 10 - 30%. They show an 
asymmetric transversal profile and spectral index dis- 
tribution. Current models suggest they are caused by 
diffusive shock acceleration of electrons in old radio- 
emitting regions, triggered by a cluster merger event. 
Currently, ^ 34 clusters are known to contain at least 
one relic source. Eleven of these also have a radio halo, 
and ten have a double relic. 

The existence of relic sources suggest cluster-wide 
magnetic fields of about 0.1 — 1 /iG and relativistic 
electrons of ^ GeV energies in large peripheral regions 
where the density of the hot ICM is low and galaxies 
are scarce. 
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Figure 17 Luminosity functions at 1.4 GHz of giant 
radio halos in the redshift range 0.1 < z < 0.2. The 
red dashed line marks the contribution from "hadronic" 
halos in relaxed clusters, while the blue dashed line 
marks the contribution from halos in merging "turbu- 
lent") clusters. The black solid line gives the total lu- 
minosity functions ( Cassano et al.[[20T2 i.The arrows 
show the EMU and NVSS sensitivities at that redshift. 



Few radio relics have been well studied and few 
have deep high resolution radio maps, so that their size 
and flux density may be under-estimated. Most radio 
spectra have been obtained with a limited frequency 
range and very few spectral index images have been 
published. Because most relics are located at the clus- 
ter periphery, they are often outside, or at the edge, of 
available X-ray images, making it difficult to compare 
radio and X-ray data. However, a few relics have been 
found and studied in det ail (e.g.[Jo hnston-H ollitt[2003t 
[van Weeren et al.|2010|[Brown et al.,201 la| . 

A few relics show a curved shape radio spectrum 
and a diffuse, filamentary morphology. They often have 
been found near the central cD galaxy suggesting a 
possible connection with a past activity of this galaxy. 



but in A 1664 (see Figure [T8]l and in A548b ( [Feretti 
[2006[ ) the relic is far from the cluster centre. 



Figure 19 compares the radio and X-ray power of 
radio halos and relics, showing that both the relic and 
halo power are correlated with the parent cluster X-ray 
luminosity, although this may be influenced by selec- 
tion effects. The relics have a larger dispersion than 
the halos, confirming that shocks are more efficient 
than turbulence in accelerating relativistic electrons or 
amplifying magnetic fields. This correlation confirms 
the link between relics and cluster mergers and shows 
that the SPARCS surveys will enable the study of large 
scale magnetic fields in low density regions at the clus- 
ter periphery. Such correlations have recently been dis- 
cussed by Feretti et al. ( 2012[ |. 
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Figure 18 Radio image (contours) overlaid on X-ray 
images (colour) of three representative clusters: (top) 
the radio relic in the cluster A 1664 ( |Govoni et al . 
(centre) the giant radio halo in A2163 



2001 1 



(Fer- 



etti et al.||2001|l, (b ottom) the double relic in A3376 
( |Bagchi et al.|2006| l. 



3.4.3 Tailed radio galaxies 

Tailed radio galaxies include wide-angle tails (WATs) 
and narrow-angle tails (NATs). They are common in 
large clusters ( Blanton et aL][2000l [200T] [2003] l, and 
appear to represent radio-loud AGN in which the jets 
are distorted by the ICM (Mao et al. "2010a|l. |Mao et 
al. ( [201 Ibj ) have found 6 tailed galaxies in the seven 
square degrees of the Australia Telescope Large Area 
Survey (ATLAS), suggesting that 40,000 new clus- 
ters will be detected by EMU, using the tailed galaxies 
as tracers. With higher spatial resolution, more tailed 

( pon] i. 



galaxies become visible, and Dehghan et al. 



using high-resolution images of the same ATLAS fields, 
find 12 tailed galaxies in 4 deg^, implying ^10^ tailed 
galaxies in EMU. Importantly, such galaxies can be 
detected out to high redshifts of z 0.5 and beyond 
( [Wing & Blanton|201 l[[Mao et al.|2010a| , providing a 
powerful diagnostic for finding clusters. 
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Figure 19 Halo (blue) and relic (red) radio power at 

1 .4 GHz in units of W/Hz versus cluster X-ray lumi- 
nosity. 

3.5 The Magnetic Sky 

Magnetic fields are very important at all scales in the 
Universe, exerting their influence at very small scales 
in the interstellar medium (ISM), at large scales in galax- 
ies and clusters, and on the vast scale of cosmic fila- 
ments. Yet we still have limited understanding of the 
magnetic contribution to the energy balance in these 
structures, and their dynamical influence and contri- 
bution to the evolution of different classes of object. 
Our best tools for measuring cosmic magnetic fields 
are radio observations of synchrotron radiation and its 
linear polarisation, which are sensitive to total and or- 
dered magnetic fields, respectively, and Faraday rota- 
tion, which is sensitive to magnetic fields along the line 
of sight. 

One of the areas which will benefit most from SKA 
pathfinder surveys is the study of radio continuum po- 
larisation. The prospects for wide-field deep imag- 
ing, with broad wavelength coverage, and (in PAF- 
equipped telescopes) with excellent control over the 
primary beam response, combine to make the suite of 
pathfinders excellent for the study of polarisation. The 
two major polarisation surveys which are currently plan- 
ned with SKA pathfinders are (i) the LOFAR Key Sci- 
ence Project on Cosmic Magnetism (see 
(ii) ASKAP-POSSUM project (see 



2.8 1, and 



2.4 1. They are 



highly complementary, spanning more than a decade in 
observing frequency, and will probe widely differing 
regimes of magnetised plasmas. Moreover, together 
they can observe both hemispheres. There are a num- 
ber of other surveys aimed at recovering high density 
grids of rotation measures (RM grids) towards differ- 
ent astrophysical regions, notably the APERTIF-BEO- 
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WULF and APERTIF-FRIGG surveys (§2.2i, the e 



MERLIN Super-CLASS survey (Battye et al. 2011 



Jansson et al.1[2009l |Nota & Katgert|[20T0l |Mao et al\ 



2.5 I and the RM grid component of the MeerKAT 



MIGHTEE Survey (§|Z9i. 

The primary tool for analysing the polarisation data 
from the pathfinders is the "RM Synthesis" technique 
(Burn 1966 Brentjens & de Bruyn 2005| l. The essence 
of this technique is to perform a Fourier transform of 
the observed complex polarisation vector from the fre- 
quency (strictly, A^) domain to the "Faraday depth" 
domain. The Faraday depth refers to the line-of-sight 
integral of the thermal electron density and magnetic 
field strength. In simple situations it is more com- 
monly known as the Faraday rotation measure (RM), 
while this simple relation breaks down for synchrotron 
emission and Faraday rotation in the same medium, 
various emitting and rotating layers along the line of 
sight, or more than one polarised feature with different 
RMs within an angular resolution element. The output 
of the RM Synthesis technique is a reconstruction of 
the Faraday dispersion function, which describes the 
complex polarisation vector as a function of the Fara- 
day depth ( [Brentjens & de Bruyn|2005||Heald|2009] l. 
Other techniques for identifying magnetic field struc- 
tures from RM synthesis observations continue to be 
developed (e.g. Beck et al.|2012 i. 

The technique allows the detection of very faint po- 
larised emission even when the RM is very high (i.e. 
corresponding to the polarisation angle rotating by more 
than a turn across the observing band) and avoids mr 
ambiguities. Better precision in the Faraday depth do- 
main is obtained by observing with wide bandwidth, 
and access to high values of |RM| is provided by high 
spectral resolution - both of which are included in the 
specifications of all of the pathfinders. 

High quality polarisation observations address di- 
verse science goals, ranging from local to cosmologi- 
cal scales. Here, we summarise a subset of these, and 
discuss the impact of the pathfinder surveys on these 
fields. 



3.5.1 The large-scale structure of the Milky Way 
magnetic field 

Knowledge of the strength and orientation of the large- 
scale magnetic field in the Milky Way is crucial for 
studies of the origin and evolution of Galactic mag- 
netism, but is also valuable for extragalactic topics such 
as ultra-high energy cosmic rays, cosmic microwave 
background polarisation, and magnetism in external gal- 
axies, galaxy clusters and intergalactic space. 

Models constrained by observational RM data from 
pulsars, extragalactic sources and diffuse Galactic syn- 
chrotron emission, attempt to determine reversals in 
the spiral magnetic pattern of the Milky Way, and to 
measure the large-scale configurations of the field in 
the Milky Way disk and halo (e.g. |Beck|201 1[ iNouFl 
Isos et al.||2008l|Van Eck et al.||2011[|Sun et al.i2008t 



2010b| l, all of which are key to understanding the evo- 
lution of galactic magnetic fields. However, contro- 
versy still exists about the existence and location of 
reversals, and the overall 3D configuration of the mag- 
netic field in the Milky Way's disk and halo. 

Major constraints to these large-scale magnetic field 
models will come from the Rotation Measure Grid of 
RMs from compact extragalactic sources, most of which 
probe Galactic Faraday rotation. The SKA pathfinders 
ASKAP, LOFAR, APERTIF, and MeerKAT have com- 
plementary programs to observe RMs of extragalactic 
point sources to construct an RM grid with an aver- 
age density of 100 RMs per square degree over the en- 
tire sky. These efforts will provide a polarised source 
density two orders of magnitude higher than the recent 
NVSS RM catalog ( [Taylor et al.T2009) , which has al- 
ready proven useful for a wide range of analyses (e.g. 
'McClure- Griffiths et al.|20T0l|Schnitzeler|20101|Stil et| 
,al..2011j . 

3.5.2 The turbulent, magnetised interstellar medium 

Magnetised turbulence in the interstellar medium (ISM) 
plays a major role in many physical processes such as 
star formation, total ISM pressure equilibrium, cosmic 
ray scattering, gas mixing, and radio wave propagation 
(Scalo & Elmegreen 2004) . 

At high Galactic latitudes, in the Galactic thick disk 
and halo, RM variations are small, and can be observed 
at long wavelengths. Spectropolarimetric observations 
of the Galactic foreground emission at mid-latitudes 
around 350 MHz, made using the Westerbork Synthe- 
sis Radio Telescope (WSRT) show a wealth of small- 
scale polarisation at angular resolutions down to a few 
arcmi n ([Haverkom et al.||2003a|b| [Schnitzeler et al.| 
[20071 |2009|l. At even lower frequencies of 150MHz, 
^Bernardi et al.| ( |2009| l observed abundant structure in 
diffuse polarisation at intermediate latitudes. There- 
fore, LOFAR RM synthesis observations at mid- and 
high-latitude are expected to discover a wealth of small- 
scale RM structures in the magneto-ionised gas. This 
will be used to characterise the power spectrum of mag- 
netic field fluctuations at high latitudes, i.e. in the Galac- 
tic gaseous halo and disk-halo connection. 

In addition, calculation of the spatial gradient of 
diffuse polarisation can allow direct visualisation of 
turbulent structures in the ISM, from which we can in- 
fer fundamental properties such as Mach number and 
Alfven speed ( [Gaensler et al.|[20lT] [Burkhart et al.[ 
[2012] ). 

Closer to the Galactic plane, the larger RM values 
will cause partial or complete depolarisation of low- 
frequency emission, but are well-suited to observations 
with ASKAP, APERTIF and MeerKAT. The RM grid 
allows statistical analysis of interstellar magnetised tur- 
bulence as a function of position ( Haverkom et al.[ 
[2008| ), while RM synthesis of the diffuse polarised emis- 
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sion allows disentangling of synchrotron emitting and 
Faraday rotating layers (Brentjens '201 1 1. It will also 



be possible to map out magnetic field strength and struc- 
ture in discrete structures such as supernova remnants, 
planetary nebulae, super-bubbles, Hll regions. Galactic 



regions. A concise review of the properties of magnetic 
fields in galaxies is provided by Beck ( 2009| l. 



chimneys, and high-velocity clouds (see, e.g.. Beck & 
|Gaensler|2004l l. 

3.5.3 Faraday Screens 

At 1.4 GHz the diffuse polarised sky typically bears Ut- 
tle resemblance to the total intensity sky (e.g. Lan- 



decker et al.'2010). A widely accepted interpretation 



for this characteristic is that the detected polarised fea- 
tures are the signature of Faraday rotation rather than 
structure in the Galactic synchrotron emission. The ex- 



ceptions are supernova remnants (§ 3.8.3 i and pulsar 
wind nebulae (§ 



3.8.41 



Observations have shown that the synchrotron emis- 
sion of external galaxies is very smooth and the WMAP 
results have confirmed that for our Galaxy too. How- 
ever, at lower radio frequencies, surveys show a lot of 
small-scale structure caused by Faraday rotation. A 
typical radio telescope operating in the 1.4 GHz wave- 
length range proves to be more sensitive to the Faraday 
rotation of an ionised and magnetised medium than to 
its bremsstrahlung signature. 

The Faraday rotation structures we find in linear 
polarisation images are called "Faraday Screens". These 
Faraday Screens are still not well understood. Ques- 
tions about these features we have to answer include: 

• Are those features created by variation of the 
magnetic field or the free electron distribution 
(or a combination of both)? 

• Can we correlate (or anti-correlate) these fea- 
tures with other tracers of the ISM? 

• How can we explain features that do not corre- 
late with anything? 

• Can Rotation Measure Synthesis help to reveal 
the characteristics of these features? 

• Can we relate these features with structures seen 
in the rotation measure distribution determined 
from extra-galactic point sources? 

A large wide-field survey such as POSSUM will help 
to better understand these features and reveal their true 
nature 

3.5.4 Nearby galaxies 

Spiral galaxies have typical magnetic field strengths of 
order 10 /xG. They typically show ordered field pat- 
terns consistent with the spiral pattern in the plane, a 
higher level of magnetic turbulence in the arms them- 
selves, and an X-shaped configuration in the high-latitude 



The WSRT-SINGS Survey has investigated the global 
magnetic field properties across a small sample of nearby 
galaxies with diffuse polarised emission (Brau n et al.| 
[2010) . This work showed that azimuthal patterns in 
polarised intensity and RM ( Heald et al. 2009 [ l are con- 
sistent with a magnetic field configuration consisting 
of an axisymmetric spiral in the disk and a quadrupo- 
lar poloidal pattern away from the plane, as predicted 
by the dynamo theory of magnetic field amplification 
(e.g. I Widrow 2002 ). A critical component of the model 
is that turbulent magnetic fields in the star-forming disk 
depolarise synchrotron radiation from the far side of 



the disk (see Bum 1966 for a description of this Fara- 
day dispersion effect). In face-on galaxies this effect, 
together with trailing spiral arms, results in the ob- 
served patterns. At higher and lower radio frequencies, 
the effects of Faraday dispersion will be reduced and 
enhanced, respectively, giving rise to different patterns 
in the polarised synchrotron radiation. 

The polarisation surveys planned with ASKAP and 
LOFAR will thus be complementary in that they will 
trace ordered magnetic fields in significantly different 
volumes within nearby galaxies. Although depolarisa- 
tion by Faraday dispersion is much stronger at LOFAR 
frequencies, this is mitigated by the fact that cosmic- 
ray electrons (CRs) that have diffused far from their ac- 
celeration sites (i.e. regions of massive star formation) 
will have lower energies because of synchrotron losses. 
Thus synchrotron radiation from this aged CR popula- 
tion will peak at LOFAR frequencies, and should trace 
ordered fields in the far outer regions of galaxies. At 
the other end of the spectrum, high frequency observa- 
tions with MeerKAT and the VLA can "see through" 
the star forming disk, and trace the fields in the star 
forming regions directly. The suite of SKA pathfinders 
can thereby provide complementary polarimetric stud- 
ies leading to an "onion-peel" polarimetric view of dif- 
fuse synchrotron emission in nearby galaxies. 

3.5.5 Magnetic field evolution in galaxies 

Beyond the nearby Universe, SKA pathfinder telescopes 
have insufficient angular resolution to map the syn- 
chrotron radiation across disks of galaxies. However, 
|Stil et al.] (2009) have shown that polarised radiation 
with a preferred, and wavelength-independent, orienta- 
tion will be detected in unresolved galaxies. Thus sta- 
tistical studies of the magnetic fields in distant galaxies 
will be possible. These same sources are also expected 
to be excellent background probes in the RM grid. 

Within the picture of the turbulent dynamo, the am- 
plitude of the regular magnetic fields in galaxies is ex- 
pected to decrease at higher redshift. The variation of 
galactic magnetic field strength with epoch may be re- 
flected in the redshift dependence of the FIR/radio cor- 
relation (see |Murphy ,2009) . 
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3.5.6 Galaxy clusters 

The ICM of galaxy clusters consists not only of hot 
thermal gas (T ^ 10^^^ K) emitting in the X-ray do- 
main, but also of magnetic fields and relativistic par- 
ticles that may give rise to synchrotron radio emis- 
sion. The most direct evidence for the presence of 
magnetic fields are the so-called radio halos and ra- 



dio relics (see § 3.4 1. Other evidence for cluster mag- 
netic fields comes from the Faraday rotation of radio 
sources that lie within or behind the cluster along our 
Une of sight (e.g. [Clarke et al.|2001[ [Johnston-Hollitt 
|& Ekers|2004l ). 

Magnetic fields are thought to play an important 
role in the development of large-scale structure in the 
Universe, but their origin and evolution are still poorly 
understood. Cosmological simulations indicate that the 
magnetic field strength resulting from the adiabatic com- 
pression of the gas falling into the cluster potential 
well cannot account for more than 1% of the magnetic 
field strength that is inferred from observations, so that 



other amplification mechanisms are required (Dolag 



et al. 2008 1. Cosmological simulations also predict 



that processes related to the cluster formation, such as 
merger events and shear flows, can in principle amplify 
the magnetic field up to the observed values. Hence, 
measurements of the ICM magnetic field, such as the 
power spectrum and radial profile, are necessary to test 
these predictions, and to measure the amplification of 
the magnetic field during cluster formation. Magnetic 
fields may also affect thermal conduction in the ICM, 
and so measuring the magnetic field strength is crucial 
to understand the origin of the relativistic particles re- 
sponsible for the radio halo and relic emission. 



3.5.7 Faraday Rotation Measures and ICM mag- 
netic fields 

A powerful method to constrain the magnetic field prop- 
erties in the ICM is to study the Faraday rotation of 
sources that are within or behind the cluster along our 
line of sight (e.g. [Krause et al. 12009 ). Recently Bonafe ie 
et al.[ ( [2010 1 have constrained the magnetic field in 
the Coma cluster by observing 7 sources at different 
projected distances from the cluster centre. The re- 
sulting RM images have been compared to the mock 



RM images obtained with the FARADAY code (Mur- 



gia et al. 2004 1 for different magnetic field configu- 
rations. The magnetic field that best reproduces the 
Coma observations has a central value Bq — 4.7 /iG, 
and declines with t he gas density profil e n according 
to B{r) = n(r)°-^. Murgia etal.' ('2004|l, Govoni et al. 



2006| l, and Vogt & Enfilin (2005) have achieved sim 



ilar results on smaller samples. These studies require 
deep multi-frequency observations of several sources 
located at different projected distances from the clus- 
ter centre, so only a small number of clusters have been 
studied. The sensitivity of current radio telescopes is 



insufficient to study a large number of galaxy clusters 
with many RM probes per cluster, but the next genera- 
tion of radio telescopes will make this possible. 

3.5.8 Fractional polarisation of radio sources and 
magnetic field in the ICM 

When synchrotron emission from a cluster or back- 
ground source crosses the ICM, regions with similar 
intrinsic direction of the polarisation plane, go- 
ing through different paths, will be subject to differ- 
ential Faraday rotation. If the magnetic field in the 
foreground screen is tangled on scales smaller than the 
observing beam, radiation with similar 'i'int but oppo- 
site orientation will be averaged out, and the observed 
degree of polarisation will be reduced (beam depolari- 
sation). In the central region of a cluster, B and n are 
higher, resulting in a higher RM, and lower fractional 
polarisation Fp. Sources at larger radii experience a 
lower RM, and so suffer from less depolarisation. 
Bonafede et al. ( 20II| l selected a sample of mas 



sive galaxy clusters and used the NVSS data (Condon 
et al. 1998[) to analyse the polarisation of radio-sources 
as a function of the projected distance from the clus- 
ter centre. They find that, statistically, the fractional 
polarisation decreases with the cluster projected dis- 
tance. By comparing this trend with that predicted by 
magnetic field models, they estimate the magnetic field 
in the cluster centre to be ~ 5 ^G. The advantage of 
this approach is that it does not need multi-frequency 
observations. This effect can be investigated with ra- 
dio surveys at a single frequency in full polarisation 
mode. Although such studies do not provide detailed 
information on magnetic fields in specific clusters, they 
allow us to understand the average properties of mag- 
netic fields in the ICM, and potentially even on larger 
scales. 

3.5.9 The origin of cosmic magnetism 

Galaxies and the ICM can only account for about one 
third of the baryon density in the local Universe ex- 
pected from a concordance cosmology. The majority 
of the missing matter is likely to reside in a warm-hot 
intergalactic medium (WHIM), which in turn is ex- 
pected to reside in the cosmic web of the large-scale 
structure. Detection of the magnetic field found in this 
cosmic web, or placing stringent upper limits on it, will 
provide powerful observational constraints on the ori- 
gin of cosmic magnetism. 

Various mechanisms have been suggested for the 
origin of magnetic fields in our Universe. One pos- 
sibility is that the fields are truly primordial, and that 
a seed field formed prior to recombination in the very 



early Universe ( Banerjee & Jedamzik,2003 ). Alterna- 
tively, fields could have been produced via the Weibel 
instability, a small-scale plasma instability formed at 
structure formation shocks ( |Medvedev et al.[[2004| i. 
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Both of these mechanisms occurred early in cosmic 
history, prior to galaxy formation, and are therefore re- 
ferred to as early-type mechanisms. In contrast, possi- 
ble late-type mechanisms are processes where the field 
was injected into the WHIM via the action of jets from 
AGN and other outflows such as galactic winds (e.g. 



Kronberg|2004l l. 



The SKA pathfinder polarisation surveys will to- 
gether generate a rich dataset from which a wide, high 
density, broad bandwidth RM Grid may be built, and 
weak magnetic fields in intergalactic filaments may be 
detected for the first time. A lower limit for the inter- 
galactic magnetic field has been determined (10^^^ G; 
IDolag et al.|20TT] i. 

In the filamentary structures of the cosmic web, 
field strengths are expected to be far more sensitive to 
the origin of the magnetic fields than in larger gravi- 
tationally bound structures where turbulent amplifica- 
tion will have led to saturation of the field. Specif- 
ically the strength and penetration of magnetic fields 
in filaments is expected to be much lower if the field 
arises from a "late" astrophysical source (AGN, galac- 
tic winds etc), compared to that arising from a primor- 
dial one, and simulations have demonstrated that the 
strength of these turbulently amplified early-type seed 



fields range typically between 0.1 and 0.01 jiG (Ryu 
|et al.|[2008| . Consequently observations of the fila- 
mentary cosmic web are invaluable for inferring the 
origin of cosmic magnetism. Although the predicted 
field strengths are extremely low, such RMs are ob- 
servationally accessible with the SKA pathfinder tele- 
scopes. 

In addition, the rotation of emission from back- 
ground radio sources through the large-scale structure 
filaments within e.g. super-clusters of galaxies can also 
be used to infer the origin of cosmic magnetism. A 
statistical analysis of these RMs to measure the power- 



spectrum of the magnetic field of the cosmic web |Ko- 
latt 199811, using cross-correlations with other large- 



scale structure indicators, will also provide stronger 
constraints than individual RMs. These measurements 
will also probe the distribution and density of thermal 
gas at the interface between clusters and the cosmic 
web, which has important consequences for a number 
of related areas in astrophysics and will be highly com- 
plementary to X-ray and Sunyaev-Zel'dovich studies 
of these structures. However, because of the low tem- 
peratures of these regions, RM measurements may pro- 
vide the only way to probe them. 

3.5.10 High redshift rotation measures 



Studies of high-redshift objects (e.g. [Kronberg et al. 
|2008[ who examined the redshift dependence of the 
spread of RMs) suggest that the amplification of mag- 
netic fields occurred more rapidly than the expectations 
of mean-field dynamo theory The conclusions of this 
type of work are so far limited by small samples, to- 



gether with uncertainties in the removal of the Milky 
Way foreground RM structure. SKA pathfinders will 
provide RM grids that will enable this kind of study 
to be done with significantly lower uncertainty. A bet- 
ter comparison with the predictions of dynamo theory 
(see e.g. |Arshakian et al.(2009) ) wiU thus become pos- 
sible. 



3.6 Cosmology 

The SKA pathfinder surveys (here referred to as SPARCS 
surveys) will not only be an important preparation for 
SKA surveys, but will also be able to deliver important 
cosmological results themselves. Earlier surveys were 
able to provide valuable constraints on cosmology and 
cosmic evolution (e.g. Blake & Wall| |2002[ |Overzjer] 
|et al.|2003j ). Here we show that SPARCS surveys will 
be able to constrain cosmological models and parame- 
ters via several probes: the Integrated Sachs-Wolfe ef- 
fect (Sachs & Wolfe|1967| l, the lensing effect known as 
cosmic magnification, and source count correlations. 

The source count correlation function (or if mea- 
sured in Fourier space, the power spectrum) is a well- 
known statistic measuring the excess probability that 
sources are clustered together on different scales. In 
the radio surveys considered here, we assume that in- 
dividual redshifts will not be available, so the power 
spectrum will be measured on the 2-D celestial sphere, 
by measuring the spatial auto-correlation function of 
the radio sources in the SPARCS surveys. 

Cosmic magnification is a gravitational lensing ef- 
fect, where the number density of background sources 
is affected by foreground structures in two competing 
ways. First, lensing dilates the background patch of 
sky, reducing the background number density. Second, 
the lensing expands the size of background sources 
while maintaining their surface brightness, resulting 
in an increase of apparent luminosity of sources, in- 
creasing the number density above any given detection 
threshold. Either of these effects can dominate depend- 
ing on the luminosity function, so there may be either 
a correlation or anti-correlation of background number 
density with foreground structures. The effect is sensi- 
tive to the geometry of the Universe and the growth of 
structures in the Universe. In our case, we correlate 
background SPARCS radio sources with foreground 
sources from shallower optical surveys, such as SDSS 
( [Abazajian et al.|2009|l f or northern surveys, Skymap 



per ([Keller et al.||2007|l or DES ( [The Dark Energy 
Survey Collaboration|2005 1 for southern surveys. 



The ISW effect is a gravitational redshift, due to 
the evolution of the gravitational potential while pho- 
tons pass through under- or over-densities in their path 
from the last scattering surface to us. In an Einstein-de 
Sitter universe, the blueshift of a photon falling into a 
well is cancelled by the redshift as it climbs out. But 
in a cosmology where there are evolving gravitational 
potentials, photons will experience a net blue- or red- 
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shift, leading to a net change in photon temperature, 
which accumulates along the photon path. This trans- 
lates into CMB temperature anisotropics proportional 
to the variation of the gravitational potentials: 



ST 



*(z) + $(z) 



dz . 



Here we have introduced the gravitational potentials 4' 
and <!>, which are the time-part and space-part pertur- 
bations of the FRW (Friedmann-Robertson-Walker) 
background in the conformal Newtonian gauge, ds^ — 
-a^{{l + 2^)dT^ - (1 - 2^)dx^). In GR, in the ab- 
sence of anisotropic stress, we expect and $ to be 
identical. 

The ISW effect can contribute significantly to the 
CMB temperature fluctuations on large angular scales, 
but it enhances only the low £ multipoles, and is smaller 
than other CMB anisotropics. For this reason, in order 
to detect the effect, we have to cross-correlate CMB 
temperature maps with tracers of the potential wells 
( Crittenden & Turok|1996j ) such as galaxies; if the po- 
tential evolves, then we should observe a correlation 
between CMB temperature anisotropics and the galaxy 
distribution. For the tests considered here, we cross- 
coiTclate the SPARCS radio sources with the CMB data 
from WMAP or Planck. 

3.6.1 Calculating probe predictions 

Here we consider the measurements of the three probes 
discussed above (source count correlations, cosmic mag- 
nification and ISW), and calculate the resulting con- 
straints obta ined with the LOFAR, WODAN and EMU 
surveys (see |Raccanelli et al.|20I I| for details about the 
surveys and the assumed distributions and bias). 

To predict the redshift distribution for these sur- 
veys, we use the latest empirical simulations devel- 
oped for the Square Kilometre Array continuum sur- 
vey ( jWiknan et al.|[2008] ). This simulation provides 
different prescriptions for the redshift evolution of the 
various populations which dominate the radio source 
counts at five different radio frequencies: 150, 610, 
1400, 4860 and 18000 MHz. Catalogues are generated 
from the S'^ databas^ using the radio flux-density lim- 
its of different surveys, and we perform our analysis by 
extracting a catalogue from the database with cuts at 
the 10 a detection levels of the forthcoming surveys. 

The simulation provides us with a source cata- 
logue in which sources are identified by type (e.g. SB, 
SFG, FRI, FRH, RQQ). The bias is computed sepa- 
rately for each galaxy population using the formalism 



of Mo & White ( 1996 ), in which each population is as- 
signed a dark matter halo mass chosen to reflect the ob- 
served large-scale clustering. With this framework, the 
increasing bias b{z) with redshift would lead to exces- 
sively strong clustering at high redshift, so the bias for 



each population is held constant above a certain cut-off 
redshift (see [Wilman et al.|2008| for details). 

We calculate the angular power spectra for each 
of the probes, i.e. source-source, radio background - 
optical foreground, and source-CMB. Each of these is 



a projection of the density power spectrum; see Rac 



[canelli et al.| ( |20I l| l for details. For instance, the ISW 
cross-power spectrum is given by 



Cf = («L«F™> = 47ry' fA'{k)W,\k)Wl{k), 

(2) 

where and are the galaxy and CMB window 
functions respectively, and A^(fc) is the logarithmic 
matter power spectrum today. 

The galaxy window function, which enters all three 
power spectra, is written as: 



dN 
dz 



biz)D{z)je[kr,{z)]dz, (3) 



where {dN / dz)dz is the mean number of sources per 
steradian with redshift z within dz, brighter than the 
flux limit, h{z) is the bias factor relating the source to 
the mass over-density, D{z) is the linear growth factor 
of mass fluctuations, ji{x) is the spherical Bessel func- 
tion of order t, and rj{z) is the conformal look-back 
time. Similar window functions exist for the CMB and 
foreground elements for the different power spectra. 

We can write corresponding correlation functions 
as a function of the angular separation 6: 



2e + i 



CiLi{cos9), 



(4) 



where Lg are the Legendre polynomials of order £. 

3.6.2 Predictions for the pathfinders 

We have calculated correlation functions or power spec- 
tra for the ISW, cosmic magnification and source counts 
for the forthcoming radio surveys. For example, in 



Figure 20 we show our result for the CMB-LSS cross- 
correlation using the predicted redshift distributions and 
bias for the LOFAR Tier 1 survey, where we assume 
the standard ACDMh-GR model and the shaded areas 
are cosmic-variance errors, calculated via: 



\ 



c. 



r<99r<TT 



{2e- 



(5) 



l)/sky 

where /sky is the sky coverage of the survey; Cf^ and 
Cf^ are the auto-correlations density-density and temp- 
erature-temperature, and contain shot noise and detec- 
tor noise; the errors are then converted to real space 
as: 
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Figure 20 Predicted correlations of the CMB with the 
LOFAR Tier 1 survey data. As a comparison we 
also show error bars from current NVSS measure- 
ments. Note the improvement we obtain with one of 
the pathfinders over current constraints on all scales. 



Now we wish to predict the combined cosmologi- 
cal measurements that the SKA pathfinders should al- 
low. In order to do this, we perform a joint analysis in- 
cluding CMB-nSNIa priors, IS W, cosmic magnification 
and the auto-correlation of radio sources. We test de- 
viations from the standard ACDMh-GR model for two 
cases: we first parameterise modifications of the theory 
of gravity; and then alternatively, we assume GR to be 
the correct model to describe gravity, and consider a 
dynamical dark energy model. 

3.6.3 Modified Gravity Constraints 

We want to test deviations from GR using the {77, /i} 
parameters, that are defined via: 



fc2^ ^ -inGa^nia, k)pA, 



(7) 
(8) 



where A is the gauge-invariant comoving density con- 
trast: 

A^d + 3^v; (9) 

in General Relativity T]{a, k) = /^(a, k) = 1, while 
in alternative models they deviate from this value, and 



can also be functions of time and scale (e.g. Zhao et 
|al.|2010| l. 

We model the time evolution of /i and 77 as (see 
Zhao et al. 2010 for details): 
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Figure 21 Predicted constraints on parameters using 
EMUh-WODAN. The outer dotted ellipse shows the 
current constraints, and the innermost grey ellipse 
shows the constraints available using EMUh-WODAN 
with the three probes described. Bottom panel: dark 
energy parameters, showing current and predicted 
68%CL constraints. Top panel: modified gravity pa- 
rameter constraints. 



In Figure|2T](top) we show predictions for how well we 
can constrain the modified gravity parameters using the 
EMU survey. Also displayed are the 68% confidence 
level for current surveys (SNe, WMAP and SDSS; 
[Raccanelli et al.|[2011| l. We use the current best-fit 
model as our fiducial model; note that if this current 
best fit remains the best fit, the standai'd GRh-LCDM 
model would be excluded with forthcoming radio sur- 
veys at high significance. In any case, the surveys im- 
pressively reduce the permitted parameter space. 

3.6.4 Constraints on Dark Energy evolution 

To study how well the SKA pathfinders will constrain 
the evolution of dark energy we use the {wq, Wq} pa- 
rameterisation ( |Luider|2003) of the dark energy equa- 
tion of state: 



w{a) — wq + Wa(l — a)- 



(12) 
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In Figure |2T] (bottom panel) we show predictions for 
the constraints on the parameters {wq, Wa}, using the 
EMU survey. Again, the constraints are substantially 
improved over the currently permitted region of param- 
eter space. 

3.6.5 Bias 

To obtain constraints on cosmological parameters from 
continuum surveys, the redshift distribution and bias 
of radio sources needs to be known. Here we suggest 
a method for measuring the bias of radio sources at 
high redshift. We cross-match FIRST radio sources 
([Becker et al. 1995) with SDSS-DR7 data ( |Abaza-| 
jian et al. 2009 ) and use photometric redshifts in the 
DR7 catalog to obtain an estimate of the redshift dis- 
tribution of the optically identified sources. To probe 
the high-z radio continuum population, we assume that 
the z-distribution of all FIRST sources is given by the 
SKADS simulation d Wilman et al.|2008| and we sub- 
tract off the distribution of the optically identified sources 
to estimate the z-distribution of the unidentified sources. 
By measuring the angular correlation function of the 
unidentified sources and comparing this with predic- 
tions for dark matter clustering, we will then be able 
to estimate the bias at z ~ 0.7. Some initial results 
using this technique are presented in Passmoor et al (in 
preparation). 

If SKA pathfinder surveys are to be used to mea- 
sure cosmological parameters, it will be important to 
continue and extend this work to obtain a better esti- 
mate of bias for different classes of radio source, and 
also to understand the sensitivity of the derived cos- 
mological parameters to uncertainty in the value of the 
bias parameter. 



3.6.6 The effect of redshifts on cosmological tests 

The cosmology discussion to this point makes the con- 
servative assumption that no redshifts are available for 
individual radio sources, although we do know the red- 
shift distribution of the different classes of source, and 

are 



An alternative approach has been suggested by |Nor-| 
^ ( 201 lb| l, using radio data alone. Sources 



ris et al 



with detected polarisation in the POSSUM survey (see 



the akeady impressive results shown in Figure 21 



calculated using this assumption. While less than about 
1% of radio sources in SPARCS surveys are likely to 
have spectroscopic redshifts at the time of the data re- 
leases, photometric and statistical redshift information 
is likely to be available for ^50% of sources, as de- 
scribed in § |4.7| 

[Camera et al.| ( [2012| l have explored the impact 
of photometric redshifts on the cosmological measure- 
ments from EMU and WODAN. Using only the galaxy 
clustering angular power spectrum and dark energy con- 
straints, they find that even approximate photometric 
redshifts obtained from SDSS and SkyMapper yield a 
significant improvement to the constraints provided by 
the radio data alone. 



2.4 1 are virtually all AGN, with an expected median 
redshift of < z >^1.8, while the remaining unpo- 
larised EMU sources will be a mixture of AGN and 
SFG, with an expected median redshift of < z >^1.1. 
The unpolarised sources can therefore be treated as a 
separate population from the polarised sources for the 
purposes of the cosmological tests described above, 
giving two radio source samples at different redshift 
ranges, thus potentially yielding a measurement of the 
time-dependent component of the dark energy and mod- 
ified gravity parameters. However, modelling is re- 
quired to establish whether this technique yields a sig- 
nificant improvement over other techniques. 

3.6.7 Weak Lensing 

Weak gravitational lensing by large-scale dark matter 
in the Universe distorts the images of faint background 
galaxies. Measurements of this "cosmic shear" effect 
as a function of redshift can provide powerful con- 
straints on dark energy and its evolution over cosmic 
time (e.g. [Albrecht et al.|2006[ [Peacock et al.|2006| l. 
The measurements require a large number of well-resolved 
galaxies to reduce the "shape noise" associated with 
the intrinsic variety of galaxy shapes. Most cosmic 
shear studies to date have therefore been at optical wave- 
lengths. Radio measurements have typically been less 
effective because of smaller number densities and poorer 
spatial resolution. [Chang et al.| ( [2004| l, however, made 
a statistical detection of cosmic shear in the VLA FIRST 



survey, and Patel et al. (20101 used VLA and MER 



LIN data to explore radio weak lensing techniques. 

Measuring weak lensing at radio wavelengths po- 
tentially offers a number of advantages over optical 
surveys in terms of minimising instrumental and astro- 
physical systematic effects. For example, radio tele- 
scopes have highly stable and well-understood point 
spread functions which makes them particularly suited 
for measuring the small distortions caused by weak 
lensing. In addition, a future weak lensing survey con- 
ducted with the SKA could yield redshifts for a signif- 
icant fraction of the lensed galaxies through the detec- 



tion of their Hi emission lines (e.g. Blake et al. 2007 ). 
Uncertainties and biases associated with photometric 
redshift errors would consequently be greatly reduced 
with an SKA lensing survey. 

A major astrophysical cause of systematic error in 
weak lensing measurements is the intrinsic alignment 
in galaxy shapes, caused by tidal torquing during the 
galaxy formation process. Cleanly separating these from 
the extrinsic alignments caused by weak lensing will 
be crucial for interpreting precision weak lensing mea- 
surements in the future. Brown & Battye ( 201 la| l 



showed that polarisation, which is routinely measured 
at radio wavelengths can, in principle, discriminate be- 
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tween the lensing distortions and intrinsic alignments. 

|Blake et al.| ( |2007| l have explored the potential 
of the SKA to perform a precision weak lensing sur- 
vey. The potential advantages described above would 
make such a survey highly complementary to future 
large-scale optical lensing surveys. The SKA pathfind- 
ers offer excellent opportunities to develop the field 
of radio weak lensing and to demonstrate the above 
techniques and advantages on real data. For exam- 
ple, large-area surveys with arcsec resolution with LO- 
FAR can test weak lensing techniques on large radio 
datasets. Moreover, the e-MERLlN telescope offers 
an exciting opportunity to perform deep sub-arcsec ra- 
dio imaging over a large field of view which will be 
an ideal test-bed for investigating radio lensing tech- 
niques. This potential has recently been recognised 
with the award of ^ 800 hours of e-MERLIN legacy 
survey time for the Supercluster Assisted Shear Survey 
(SuperCLASS), which will image a 1.75 deg^ region 
of sky containing a supercluster to a depth of 4 /iJy/beam 
with 0.2 arcsec resolution, including polarisation infor- 
mation. This study is expected to detect ^ 1-2 galaxies 
per sq. arcmin, which will be used to map the dark 
matter distribution in the vicinity of the supercluster 
( |Brown & Battye|201 lb| l. 



3.7 Variability and Transients 

Transient and variability studies are key science drivers 
for SKA pathfinder instruments, and potentially for the 
SKA as well. Here we briefly discuss some of the 
key issues that need to be considered and addressed on 
the road to the SKA, potentially via pathfinder instru- 
ments. We also highlight some of the outcomes that 
can potentially be achieved through future surveys. 

Blind and dedicated radio transient surveys can po- 
tentially 

• discover new classes of energetic objects 
( 'unknown-unknowns ' ), 

• monitor known transient and variable radio sources 
( ' known-knowns ' ), 

• search for previously undetected transient and 
variable sources ('known-unknowns'), 

• explore the intervening medium via extreme scat- 
tering events, and 

• search for a SETl signal. 

Table[T]summarises our current knowledge of vari- 
able radio sources (see also [Condon et al.| ( |1979| l 
and references therein), showing observations at a va- 
riety of cadences, frequencies and Galactic latitudes. 
A complex variety of variability is reported; in some 
cases short time-scale variations at low frequencies are 
attributed to interstellar scintillation, while long-term 



variations at GHz frequencies are in some cases at- 
tributed to AGN activity. A systematic experiment to 
map the surface density of variable radio sources as a 
function of flux, frequency and cadence would be an 
invaluable step on the road to the SKA. 

Source variability (including bright calibrators, see 
[Gaensler & Hunstead|2000] l may affect the calibration 
and clean strategies required for future surveys and will 
ultimately impact on the final dynamic range and im- 
age fidelity (e.g. see [Stewart et al.pOlT] ). Therefore 
characterising the variability and abundance of known 
radio sources is a necessity for planning future sur- 
veys and calibration strategies. Many of the compact 
nJy sources to be studied by the SKA will be affected 
by scintillation (Dennison & Condon 1981 Rickett 
'1986 ). Furthermore, at these sensitivities highly vari- 
able sources in the range 1/iJy ImJy will yield a dy- 
namic range of 1 000: 1 (and greater). CarilU et al. ( 2003 | l 
report a surface density of highly variable sources of 
up to 18 deg^-^ (see Table [T]) at a detection threshold 
of ^0.1 mJy: such a high abundance could seriously 
impact SKA surveys. 

Possible transient sources range from extrasolar plan- 
ets at frequencies below 50 MHz ("OrieBmeier et al. 
|2007) to brown dwarf flares at 8 GHz (Berger||2002| l. 
Early blind searches for transient sources suggest that 
a bright (> O.lJy), frequent, population of GHz tran- 
sients does not exist (see e.g. Croft et al.|2010 Bower 



et al.|20101 [Bower & Saul|201 Ij [Croft et al.|201 lt|Bell 



et al.|20lT I. The benchmark survey of (Bower et al. 



2007]) which predicted a substantial number of radio 



transients in the sub-mJy range has recently been re- 



vised by Frail et al. (2012l, lowering the predicted 
surface density of sources. An all-sky survey could 
however still potentially detect several thousand new 
transients, though a much lower yield would also be 
consistent with the current uncertainties (see e.g. [Ban-| 
nister et al.|20n]purphy et al.|20l"2) . 

So far many surveys have been conducted in the 
GHz regime, at low frequencies, synchrotron sources 
may be optically thick with lower peak fluxes. Coher- 
ent processes therefore might dominate transient de- 
tections below 1 GHz, and rate estimates should dif- 
ferentiate between the two processes. Exploring the 
transient rates at low frequencies is an especially im- 
portant goal for pathfinders such as LOFAR and the 
LWDA ( )Fender et al..2008,,Lazio et al.|2010) and ul- 
timately SKA-Low. 

A goal for the transient key science projects is to 
sample a given field on a logarithmic range of time- 
scales, perhaps commensally with a project revisiting 
a given field at a variety of time-scales to build up the 
necessary sensitivity. A further goal of the variability 
studies will be to follow up on objects (such as pulsars 
or radio stars) which may be detected in the large con- 
tinuum surveys, and may be manifested in the image 



as having spoke-like features, as discussed in S 3.8.4 



When the EMU and WODAN surveys start, transient 
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and variability studies through projects such as VAST 
( [Chatterjee eTar '2010'; Murphy et aL] [2012) can be 
performed commensally via comparison with previous 
catalogues such as SUMSS (|Bock et al.||1999|l and 



NVSS ( [Condon etal.|1998 V Through commensal and 
targeted observations a more complete census of the 
abundance of transients and variables can be achieved 
prior to SKA operations. 

3.8 Galactic science 

The proposed wide-field surveys such as EMU and POS- 
SUM include the Galactic plane, giving the opportu- 
nity to create a sensitive atlas of discrete Galactic radio 
continuum sources. Most of these sources represent 
the interaction of stars at various stages of evolution 
with their environment. The high sensitivity of the sur- 
veys allows access to all stages of Hll regions from the 
optically thick hyper-compact kind presumably related 
to the earliest stages of a "turned-on" star to large Hll 
regions related to massive mature or giant stars pro- 
ducing strong stellar winds. Due to the high resolu- 
tion of these surveys we can detect the youngest, most 
compact supernova remnants and pulsar wind nebu- 
lae. We can detect and study planetary nebulae giv- 
ing us insight in the late stages of the evolution of 
low mass stars. And we can detect radio stars and 
pulsars. High sensitivity polarisation surveys such as 
POSSUM will also open the door to the identification 
of currently unexplained polarisation features known 
as Faraday screens. 

At 1 .4 GHz the main mechanisms for radio emis- 
sion from discrete Galactic objects are thermal free- 
free emission from Hll regions and planetary nebulae 
and non-thermal synchrotron emission from supernova 
remnants and pulsar wind nebulae. Resolved Galac- 
tic sources can be separated by their spectral and lin- 
ear polarisation properties. Young HII regions such 
as compact, ultra-compact, and hyper-compact sources 
are optically thick at frequencies around 1.4 GHz, and 
thus easy to separate from not only Galactic but also 
extragalactic sources even if they are not resolved. Large 
evolved Hll regions have flat optically thin spectra com- 
pared to supernova remnants that may show similar 
structures but have steeper spectra, and are typically 
linearly polarised. Pulsar wind nebulae also have flat 
radio spectra, but are linearly polarised. Comparisons 
with surveys at other wavelength such as infrared or X- 
ray will further help to discriminate individual sources 
of interest for follow-up studies. 



3.8.1 HII regions 

The high sensitivity of the planned wide-field 1 .4 GHz 
surveys such as EMU and WODAN allows access to 
all stages of the evolution of Hll regions, even though 
free-free emission near 1 .4 GHz will be optically thick 



for the very dense ultra- and hyper-compact HII re- 
gions (UCHII and HCHII, respectively), which are too 
small to be resolved by a lOarcsec beam, unless they 
are very nearby. UCHII and HCHII represent a very 



young phase in the development of an HII region ( Chur ;h- 
iwell,2002j . 




10 
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Figure 22 Predicted peak radio flux density of unre- 
solved UCHIIs and HCHIIs at 1 .4 GHz as a function of 
distance from the sun. UCHIIs and HCHIIs with dif- 
ferent diameters are displayed. We assumed a constant 
electron temperature of 8000 K, a homogenous distri- 
bution of material, and high optical depth (r 3> 1.0). 

Expected flux densities at 1.4 GHz for the most com- 



pact HII regions are displayed in Figure 22 The sep- 
aration between ultra-compact and hyper-compact HII 
regions is at a diameter of about 0.02 pc. Hence, EMU 
will be able to detect all ultra-compact HII regions and 
also many hyper-compact objects. There is still con- 
fusion about what a hyper-compact HII region actually 
is. Are they just extremely compact ultra-compact HII 
regions or a new class of objects? There is evidence 
that HCHIIs have somewhat different spectral charac- 
teristics than the UCHIIs, that might be explained by a 
density gradient within the ionised gas ( Kurtz [2005 | l. 
So far only a small number of HCHIIs are known. A 
more complete census of these objects with a sensi- 
tive wide field survey is required to get a more com- 
plete picture of the evolutionary process of HII regions. 
Follow-up observations at higher radio frequencies or 
a comparison with high radio frequency surveys is then 
required to determine physical characteristics which 
will help to unravel the true nature of these objects. 

We can also probe the magnetic field inside HII re- 
gions by either studying the rotation measure signature 
of polarised point sources "shining through" the ob- 
ject (e.g. Harvey-Smith et al. 201 l[ l or by studying 
overlapping or even related diffuse polarised emission 
(e.g. foster et al.|2006[ ). 
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Table 1 A selection of variable source statistics taken from the literature, p gives the snapshot rate of sources 
(deg~^). tchar gives the characteristic timescale on which the variability was sampled. ^ gives the fractional 
change in flux (please refer to individual publications for further details). 



_s_ 



Study 



Flux (mJy) 



p(deg ^) 



tr 



V (GHz) 



B annister etal J 20 1 1 1) 
CariUi et aTT poteF 



Frail et al. (20051 



Becker et al. 



Gaensler & Huiistead' ('2000') 
Condon & Backer (1975) 



Dennison et al. (T98T 



Simonetti ^Tordes (1990i 



Ryleet alj ([T97g 



I'aylor & Gregory (1983i 



>14 

>0.1 

>0.1 
>0.25 
> 2500 

- 1000-25000 

- 500-33000 
400-12000 
200-2000 

18-1200 



0.268 
<18 
1.6 
5.8 



days - years 
19days & 17mths 
^15 years 

^ 1 day 
days - years 

days 
5 to 10 years 
days 
months 
days-months 



0.843 
1.4 
4.8 
5 and 8.5 
0.843 
2.695 & 8.085 

0.318 
0.820 & 1.41 
2.7 to 15.4 
5 



>50% 
>50% 
>50% (d) 

>50% 
<20% (d) 
0.5% & 0.98% (a) 

8-100% 
4.1% & 3.5% (b) 
10%-50% (c) & (d) 
10%-400% (c) & (d) 



(a) 0.5% at 2.695 GHz and 0.98% at 8.085 GHz. Values are derived from the average (daily) fractional change in 
16 sources. 

(b) 4.7% at 0.820 GHz and 3.5% at 1.41 GHz. Values are derived from the average modulation index ( nns/mean) 
in 13 flat spectrum sources. 

(c) also see \Rickett\ ^1986^ 

(d) Observations in the direction of Galactic plane. 



3.8.2 Planetary nebulae 



3.8.3 Supernova Remnants 



Only 274 supernova remnants ( SNRs) have been dis 



covered so far in our Galaxy (Green 2009f 



Planetary nebulae (PNe) are the most abundant com- 
pact Galactic sources in the NRAO VLA Sky Survey 
(NVSS, [Condonetal.|1998| l. [Condon & Kaplan] ( [1998) 1 
identified 680 NVSS radio sources brighter than 2.5 mJy 
with planetary nebulae. Assuming a detection limit of 
50/iJy, EMU will be able to detect unresolved PNe 
up to seven times farther away than NVSS, potentially 
leading to a significant increase of the number of PNe 
detected at radio frequencies. PNe can be very useful 
tools for measuring extinction and probing star-formation 
rates of stars not massive enough to produce super- 
novae or HII regions ( Condon & Kaplan" 1998', 'Con- 
[don et al.|1999| l. Identification of these PNe, however, 
will be difficult. 

The ionised low density circumstellar material of 
nearby planetary nebulae, released during the AGB wind 
phase of the central stellar object, is illuminated by 
background linear polarised emission through Faraday 
rotation. Free thermal electrons embedded in a mag- 
netised medium rotate background linear polarisation 
emission to create a rotation measure signature observ- 
able at 1.4 GHz. A typical radio telescope operating 
in this wavelength range proves to be more sensitive 
to the Faraday rotation of a plasma region than to its 
bremsstrahlung signature. Linear polarisation studies 
of nearby planetary nebulae can reveal important infor- 
mation about the mass loss history of old AGB stars or 
white dwarfs that cannot be obtained easily any other 
way (e.g. [Ransom et ar|2008,,2010| ). 



How- 
ever, the estimated population of SNRs is much higher 
(maybe 500-1000; He lfand et al.|2006 [). This discrep- 
ancy is mainly caused by a strong bias towards bright 
extended objects. These are the mature SNRs expand- 
ing into a medium to high density environment. There 
are two groups of missing SNRs. The compact, pre- 
sumably very young, objects were missed by previous 
wide-field surveys, because they are either confused 
by other nearby objects or mistaken for extragalactic 
sources if they are not resolved. The second group 
of missing SNRs are low surface brightness SNRs. A 
population of more than 1000 SNRs for our Galaxy 
can also be extrapolated using the Canadian Galactic 
Plane Survey (CGPS, Taylor et al. "2003) and its su- 
pernova remnant catalogue ( Kothes et al.^2006) . The 
CGPS is a radio continuum and HI survey of the Outer 
Galaxy. Since the distance to the edge of the Galaxy, 
and hence to every object in the survey, is very short, 
the ^ 1 arcmin resolution at 1420 MHz radio contin- 
uum (< Ipc) is sufficient to resolve any SNR - even 
those which are very young. 

EMU and POSSUM will have a very similar spatial 
resolution than the CGPS but much higher sensitivity. 
The 10 arcsec beam translates to about the same spa- 
tial resolution at the most distant edge of the Galaxy 
through the inner part of the Milky Way (w 20kpc) 
than the CGPS beam does towards the Outer Galaxy. 
Hence, we would expect to detect a similar number 
density of SNRs. This indicates that EMU in con- 
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junction with POSSUM could discover more than 500 
new SNRs. Resolved SNRs can be easily distinguished 
from extragalactic sources and Galactic Hll regions by 
an analysis of structure and spectral and polarisation 
signatures. 

Radio polarisation observations of SNRs can also 
be a powerful tool to probe the magnetic field into which 
these objects are expanding. SNRs act as magnifying 
glasses of ambient magnetic fields, since they freeze 
and compress it into the expanding shell of the swept 
up material. Radio polarisation and rotation measure 
studies can reconstruct the 3-dimensional ambient mag- 
netic field, as demonstrated in the studies by[U yaniker, 
Kothes, & Brunt] (|2002ll; |Kothes & Brown| (,2009j); 



Harvey-Smith et al.| (2010|l. Radio polarisation stud 



ies of SNRs with EMU and POSSUM can probe the 
large-scale magnetic field of the Inner Galaxy, detect 
potential field reversals, and probe the transition from 
the Galactic plane to the halo. 



3.8.4 Pulsars and their wind nebulae 

Pulsars generate magnetised relativistic particle winds, 
inflating an expanding bubble called a pulsar wind neb- 
ula (PWN), which is confined by the expanding super- 
nova ejecta. Electrons and positrons are accelerated 
at the termination shock some 0.1 pc distant from the 
pulsar and interact with the magnetic field to produce 
synchrotron emission across the entire electromagnetic 
spectrum. At GHz frequencies the flat spectrum PWNe 
stand out from steep spectrum SNRs and can be distin- 
guished from Hll regions by their linear polarisation 
signal. 

The high resolution and sensitivity of the proposed 
wide field 1.4 GHz surveys will enable us to built a 
more complete census of Galactic pulsar wind nebu- 
lae down to unprecedented levels. 

The number of pulsars at high Galactic latitudes is 
poorly constrained since most surveys concentrate on 
a narrow band along the Galactic plane. Out of the 
1433 pulsars with measured 1400 MHz fluxes listed in 
the ATNF Pulsai- Catalogue ( [Manchester et al.|2005"] l, 
1410 have integrated fluxes above 50 /iJy. Only 310 of 
those have Galactic Latitudes |6| > 5°. Sensitive high 
resolution radio continuum surveys will enable us to 
get a more complete picture of the pulsar distribution 
in the Milky Way Galaxy and their rotation measures 
will help to constrain the configuration and strength of 
the large-scale Galactic magnetic field. 

Pulsars should be distinguishable from other point- 
like radio sources because of their steep radio contin- 
uum spectra, their polarisation, and artefacts caused by 
short term variability. These observational signatures 
will be complemented by the variability and transient 



3.8.5 Radio stars 

Radio stars emit a tiny fraction of their total luminos- 
ity in the radio band. For example, the quiet Sun has 
a radio luminosity of 10^^ W, which is only ^ 10^ ^•^ 
of its bolometric luminosity. Nevertheless, in many 
cases, radio observations of stars and stellar systems 
have revealed astrophysical phenomena, not detectable 
by other means, that play a fundamental role in our 
understanding of stellar evolution and of physical pro- 
cesses that operates in stellar atmospheres. 

Generally, the brightest stellar radio emission ap- 
pears to be associated with magnetically-induced phe- 
nomena, such as stellar flares, related to the presence 
of a strong and/or variable stellar magnetic field (high 
brightness temperature) or with enhanced mass-loss 
(large emitting surface). 

Among the brightest radio sources are active stars 
and binary systems, including flare stars, RS CVns, 
and Algol binary systems, characterised by strong mag- 
netic activity which drives high energy processes in 
their atmospheres. Their radio flux density is highly 
variable and it usually shows two different regimes: 
quiescent periods, during which a basal flux density of 
few mJy is observed, and active periods, characterised 
by a continuous strong flaring which can last for sev- 
eral days ( Umana et al. 1 1 995) l . 

Both quiescent and flaring radio emission show spec- 
tral and polarisation characteristics consistent with non- 
thermal radio emission, probably driven by the mag- 
netic activity manifested at other wavelengths. In this 
scenario, the radio flux is due to gyrosynchrotron emis- 
sion ( |Guedel|[2009i , caused by the interaction of the 
stellar magnetic field with mildly relativistic particles. 
A similar mechanism causes the radio emission from 
pre-main sequence (PMS) stars and X-ray binaries. Non- 
thermal radio emission is also seen from shocks of col- 
liding winds in massive binaries. 

Stellar radio flares can occur also as narrow band, 
rapid, intense and highly polarised (up to 100%) radio 
bursts, that are particularly common at low frequency 
(< 1.5 GHz). They are believed to be the result of co- 
herent emission mechanisms, requiring a strong mag- 
netic field (which may be variable) and a source of en- 
ergetic particles. Coherent burst emission has been ob- 
served in RS CVns, flare stars. Brown dwarfs (BD), 
and chemically peculiar stars (CPs) ( |Slee et al.|2008t 
Berger|2002 Trigiho et al.|2000 i. Coherent emission 
has been detected in just a few tens of stars, because of 
the limited sensitivity of the available instruments. 
Stellar radio emission often exhibits significant cir- 



surveys discussed in p. 7 



cular polarisation ( [Trigilio etal. 2011 Ravi et al. 2010 1. While 
some extragalactic sources do exhibit circular polar- 
isation (e.g. |Rayner et al.||2000| l, the degree of po- 
larisation is small (typically ^ 0.01%), compared to 
the levels up to 100% seen in radio stars. It is there- 
fore important, for the purposes of distinguishing stars 
from extragalactic objects, that surveys record circular 
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Figure 23 Typical radio spectrum of several classes of 
radio emitting stars. Fluxes have been derived from the 
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assuming an appropriate distance for each type of radio 
star (10 pc for flares stars, including late M-L, 100 pc 
for active binary systems, 1 kpc for supergiants, OB 
and WR, 500 pc for CP stars). The bandwidth and sen- 
sitivity of EMU have been also indicated. 



as well as linear polarisation. 

Thermal emission (bremsstrahlung emission) is ex- 
pected from winds associated with WR and OB stars, 
shells surrounding PNe, and Novae and jets from sym- 
biotic stars and class PMS stars (Gu edel 2002). In all 
these cases, radio continuum emission is an important 
diagnostic to understand the underlying physics. For 
example, radio continuum observations are the most 
precise way yet to determine the mass-loss of stellar 
winds. This is particularly important when other di- 
agnostics cannot be used, such as in the case of dust 
enshrouded objects ( Umana et al.|2005) l. 

Since no sensitive survey for radio stars has yet 
been conducted, we have little statistical information 
on stellar radio emission is available. Forecasts of the 
number of sources detectable with SPARCS surveys 
rely on radio observations of small samples of stellar 
sources, usually selected on the basis of observed pecu- 
liarities at other wavelengths. Assuming typical radio 
luminosities for each type of radio star, we conclude 
that stellar winds and non-thermal radio emission from 
many active binaries, flares stars, PMS will be eas- 
ily detected with an rms sensitivity of 10 /iJy/beam 
(EMU), while SKA (1 ^Jy/beam rms) will be able to 
detect a quiescent Sun at 10 pc (see Figure [23]l. 

Until now, targeted observations of well known ra- 
dio stars have constituted the best approach to investi- 
gate their radio emission mechanisms (e.g. flare devel- 
opment, spectral and polarisation evolution, emission 
mechanism, etc.). An all-sky survey would signifi- 
cantly enlarge the known stellar radio emitting database, 
free from selection effects. Results from such a survey 



would provide new insights into the physics of active 
stellar systems and plasma processes. 

Two areas of stellar radio emission will particu- 
larly benefit from SKA and SPARCS surveys. One 
is the study of gyrosynchrotron stellar flares, where 
multi-epoch, multi-frequency observations will enable 
the detection of serendipitous flaring activity, allowing 
the derivation of a typical behaviour (occurrence rate, 
variation amplitude, etc.) from a statistical study of 
a larger source population. Moreover, detailed stud- 
ies of a large number of stellar coronae will be pos- 
sible, allowing us to understand the nature of energy 
release in the upper atmospheres of stars of different 
mass and age and to investigate the correlation between 
radio and X-ray emission and thus on the Neupert ef- 
fect (Guedel 2009). The expected sensitivity will also 
test the presence of super-flares on solar-type stars, and 
search for correlations with an orbiting hot Jupiter ( |Ma^ 
|hara et al.|20T2l ). 

The other important area is related to the study of 
coherent emission. The EMU plus WODAN surveys 
will offer the best opportunity yet to determine how 
common coherent radio emission is from stellar and 
sub-stellar systems. The detection of coherent emis- 
sion from a large fraction of surveyed stars will have 
immense implications for our understanding of both 
stellar magnetic activity and the dynamo mechanism 
generating magnetic fields in fully convective stars and 
brown dwarfs (Halhnan et al..2008,,Ravi et al.|20l"T] l. 
Coherent emission observed in active and BD stars 
shares several characteristics with that observed in CP 



stars ( Trigiho et al. 20 00), since both require a large- 
scale magnetosphere, and are similar to the low fre- 
quency coherent radio emission observed from the mag- 
netised planets in our solar system ( [Trigilio et al. 201 l| l. 
The well-known topology of CPs magnetic fields, in- 
dependently derived from optical observations, makes 
this kind of object an excellent laboratory for stellar 
magnetosphere studies. Matching of predicted emis- 
sion with that observed yields parameters of the radio 
source, such as the surface magnetic field, the num- 
ber density of the emitting electrons, and the energy 
spectrum of the electrons. In turn, this provides clues 
to the acceleration process, the size of the magneto- 
sphere (Alfven radius) and the inclination of the rota- 
tional axis. 

If coherent emission is present in many radio active 
stars, with the same characteristics, it will constitute 
an excellent diagnostic for star magnetospheres, and 
a powerful probe of magnetic field topology. In one 
CP star (CU Vir) the coherent emission is stable on a 
time-scale of years, and has been used to time the rota- 
tion of the star, revealing a likely change of its period 
( [Trigilio et al.|2008[|Ravi et al.|2010| . The discovery 
of other similar radio lighthouses will enable high pre- 
cision studies of the rotation period, and thus angular 
momentum evolution, in different classes of star. 
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3.9 Unexpected discoveries 

Experience has shown that many great discoveries in 
astronomy have been made, not by testing a hypoth- 
esis, but by observing the sky in an innovative way. 
The necessary conditions for this to take place are (a) a 
telescope observing an unexplored part of the observa- 
tional phase space (frequency, resolution, time-domain, 
area of sky, etc), (b) an intelligent observer who under- 
stands the instrument sufficiently well to distinguish 
between artefact and discovery, (c) a prepared and en- 
thusiastic mind ready to accommodate and interpret a 
new discovery. 

SKA pathfinder surveys will easily satisfy (a), if 
only in terms of the numbers of objects surveyed. How- 
ever, their petabytes of data, and arms-length access, 
may prevent an observer from satisfying (b) and (c). 
Although we may hope that someone will eventually 
stumble across any unexpected phenomena in the data, 
the impenetrable size of the database implies dark nooks 
that may never be fully explored. If we rely on serendip- 
ity, discoveries may remain undiscovered, for ever. 

An alternative is to harness data-mining techniques 
to help the intelligent observer search for the unex- 
pected. For example, the WTF (Widefield ouTlier Finder) 
(e.g. Norris et al. 201 lb) project will systematically 
mine the EMU database, searching for the unexpected 
by discarding objects that already fit known classes of 
object, using a variety of approaches including: 

• decision tree analysis 

• cluster analysis 

• k-furthest neighbour (kFN) by analogy with the 
kNN technique ( |Zinn et al.|2012b[ ) 



• Bayesian analysis 

Identified objects/regions will be either processing 
artefacts (which are themselves important for data qual- 
ity control), statistical outliers of known classes of ob- 
ject, or, in a few cases, genuinely new classes of object. 



4 Technical Challenges 
4.1 Survey design & quality 

The SKA pathfinder radio continuum surveys may be 
characterised by the following parameters: 

• point-source detection limit, 

• angular resolution, 

• brightness sensitivity, 

• dynamic range, 

• frequency coverage, 

• polarisation, 

• accuracy: 

- flux densities, 

- positions, 

- polarisation, 

- image fidelity, 

- uniformity in sensitivity. 

All of these need to be addressed for these sur- 
veys. Ironically, the more ambitious a survey, the more 
susceptible it will be to small systematic errors. The 
ambition of the surveys therefore needs to be matched 
by a corresponding level of effort to understand subtle 
sources of error and uncertainty. For example, many 
radio-astronomers have been content to tolerate 10% 
uncertainty in their measurements of flux density, while 
we know that our instruments are, in principle, capable 
of measuring flux densities to 1%. Reasons for this 
failure include deconvolution errors, errors in primary 
beam shape and pointing, failing to correct for miss- 
ing extended flux, resolution bias, bandwidth smear- 
ing, time smearing, CLEAN bias, Eddington bias, etc. 
To overcome these effects (and probably find new ef- 
fects which we have not yet met) will doubtless require 
a significant amount of effort to understand the instru- 
ment and its errors, but that effort is justified by the 
scientific return from these enormous surveys. 

It will also be important to present the data to users 
in ways that are clear and unambiguous, with uniform 
coverage and easily understood survey parameters. For 
example, most surveys produce catalogues chosen at 
a ^ 5(7 level of confidence. For some purposes, a 
less reliable catalogue at a 3 tr level (for greater num- 
bers of sources) may be useful, while for other studies 
a threshold of 8 tr or 10 cr level (for greater reliabil- 
ity) may be more appropriate (e.g. |Raccanelli et al.| 
201 H assumes 10 a). While an 8 cr subset can always 
be chosen from a 5 cr sample, it may cause less con- 
fusion amongst users for the survey to generate these 
different levels of catalogue. 
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It is also important to remember that these surveys 
will be used for purposes that haven't yet been thought 
of. For example, one of the most exciting results from 
the NVSS was the measurement of the RM sky (and 



served by all existing and new radio telescopes. Where 
possible, they have been chosen to overlap with a field 
that is well-studied at other wavelengths, to maximise 
the science to be obtained from these observations. The 



hence the Galaxy's magnetic field) (Taylor et al. 2009 1. 
That was not planned at all during the survey, but was 
subsequently possible because the survey had excellent 
continuum calibration. 

Management of large surveys is likely to be an is- 
sue, with previous large radio surveys having had rela- 
tively few people with the responsibility for conducting 
the survey in contrast to the emerging model in which 
large teams are responsible. 

4.2 The SPARCS Reference Fields 



The uniformity required from the SPARCS surveys can- 
not be achieved by simply ensuring that the calibrators 
are on the same flux scale, although this is an essen- 
tial first step. For example, the fluxes of strong sources 
in flie first ATLAS data release ( |Norris et a"L||2006[ 
Middelberg et al.|[2008a| l agreed with those measured 
by the VLA, while weaker ATCA sources were sys- 
tematically weaker than weak VLA sources, an effect 
which turned out to be caused by incorrect assump- 
tions about bandwidth smearing in mosaiced ATCA 
data (Hales et al .|2012a I. Other subtle flux-dependent 
effects doubtless await us, and consistency between 
our calibrator sources will not be sufficient to detect 
or fix them. There have even been suggestions of dec- 
lination-dependent calibration effects (Parra et al. 2010[ l. 

During the commissioning of the SPARCS tele- 
scopes, particularly those using the relatively untried 
PAF technologies, it will be important to repeatedly 
survey well-studied fields containing a grid of sources, 
pushing the sensitivity and dynamic range until all the 
effects are well-understood and can be corrected for. 

The measured flux densities of sources in radio sur- 
veys are subject to a large number of subtle corrections 
and bias effects, and even expert radio-astronomers do 
not always agree on these. For example, there is a 
continuing debate on the measured scatter in the radio 
source counts at low flux densities (see Fig. [3]l, with 
some proponents arguing that it is caused by cosmic 
variance, and others arguing that it is due to processing 
and calibration differences. In some cases, reprocess- 



ing or reobserving (e.g. Condon et al. 2012 1 has shown 
that some variations are caused by different processing 
paths, although cosmic variance may be sufficient in 
other cases. The problem may be even worse for EMU 
and WODAN, as both rely on new phased array feed 
(PAF) technology. To extract the maximum science 
output from these major projects, relying on new PAF 
technology that has yet to be sufficiently battle tested, 
requires careful characterisation of individual survey 
systematics. 

To achieve this, we propose three reference fields at 
declination - +30°, 0°, and -30°), which can be ob- 



two Southern fields are in the CDFS-ATLAS field (Nor 



ris et al. 2006) and the COSMOS field ( |Schinnerer et 
al.||2007| . A Northern field (15:30, H-29:00) has been 



chosen within the overlap range of ATCA and Wester- 
bork, within the SLOAN survey area, and at least 2.5° 
away from any strong (>0.5 Jy) sources. It is planned 
to observe all three fields, as far and deeply as possi- 
ble, with all existing survey telescopes (VLA, ATCA, 
Westerbork, LOFAR, GMRT, etc.) as well as the new 
SPARCS surveys. 

4.3 Calibration 

Calibration and imaging of SKA pathfinder surveys 
is likely to present significantly new challenges com- 
pared to earlier surveys. For example: 

• Direction-dependent calibration is likely to be 
important not only for low frequencies, such as 
LOFAR, but also for intermediate frequencies, 
such as EMU and probably even for VLA and 
MeerKAT. 

• Because all data are acquired in a spectral-line 
mode, flagging the data forRFI affects the weight- 
ing of the data. 

• The time scale of gain variations, such as those 
resulting from pointing errors, is important. For 
example, wind gusts have a particularly dramatic 
effect on pointing errors, in part because the er- 
rors scale as the square of the wind velocity. The 
result is quite rapid pointing variations, which 
may be difficult to track and correct. 

• Experience from multi-epoch observations with 
the WSRT has illustrated the importance of mon- 
itoring the calibration information over time. In 
particular, there can be time variability of the 
calibration information (e.g., a variable source in 
the field of view of the calibrator), which could 
affect the calibration or result in apparent vari- 
ables. 

• The scientific drivers for high dynamic range imag 
ing need to be identified and kept clear. There 
may not be a need for uniform and high dynamic 
range across the entire sky. 

4.4 Imaging 

At the excellent sensitivities offered by SKA pathfinder 
telescopes, imaging and calibration algorithms need to 
deal simultaneously with issues of 
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• wide-bandwidth (e,g, primary beam corrections 
and polarisation leakage terms change with fre- 
quency), 

• wide-field imaging, 

• advanced multi-scale deconvolution techniques 
for imaging extended emission, 

• directionally-dependent effects such as instrumen- 
tal and ionospheric errors, 

• the need for high dynamic range and high polar- 
isation purity. 

Furthermore, the inherent data volumes are large 
(typically petabytes) so that computational efficiency 
and the number of data traversals are important param- 
eters to consider in designing post-processing strate- 
gies. 

4.4.1 Existing algorithms 

Two distinct approaches are being pursued in the com- 
munity to correct for wide-field and direction-dependent 
effects, which can be classified as (a) projection algo- 
rithms and (b) faceting algorithms (Bhatnagar et al. 
|2009l ). 

Projection algorithms are based on physical mod- 
elling of the various direction dependent (DD) terms 
of the measurement equation and incorporating them 
as part of the forward and reverse transforms. These 
algorithms therefore do not require assumptions about 
the source brightness distribution. This allows them to 
take advantage of the FFT algorithm for Fourier trans- 
forms, as well as integrate well with advanced decon- 
volution techniques for scale-sensitive wide-band imag- 
ing. 

Faceting algorithms on the other hand are based 
on data partitioning in such a way that standard di- 
rection independent techniques can be applied to the 



partitions ( Nijboer & Noo rdam 2007, Mitc hell et al. 
2008 Smimov||201 1 1. This approach is phenomeno- 
logical by design requiring no physical understanding 
or modelling of the DD term in the measurement equa- 
tion. They can however suffer from non-optimal use 
of the available signal-to-noise ratio (due to data par- 
titioning), curse of dimensionality (too many degrees 
of freedom) and computational load. To alleviate some 
of these problems, in addition to the problem of de- 
grees of freedom, assumptions about the structure of 
the brightness distribution typically have to be made. 

The W-Projection algorithm ( jComwell et al.|2008) l 
incorporates the effects of non co-planar baselines and 
corrects for its effects during imaging. The A-Projection 
algorithm ( Bhatnagar et al.|[2008[ l similarly accounts 
for the time, frequency and polarisation dependence 
of the antenna primary beams modelling the antenna 



aperture illumination patterns (or the measured illumi- 
nation patterns). The combined AW-Projection algo- 
rithm therefore can be used at low frequencies where 
the effects of the W-Term as well as that of antenna 
primary beams limit the imaging performance. An- 
tenna pointing errors are thought to limit the imaging 
performance at high frequencies and in mosaic imag- 
ing at all frequencies. Using A-Projection for forward 
and reverse transforms, the Pointing SelfCal algorithm 
(Bhatnag ar et al.||2004 l solves for the time dependent 
antenna pointing errors. The pointing solutions are 
then incorporated as part of the A-Projection to correct 
for the effects of antenna pointing errors. 

The Mul ti-Scale Clean (jCornwell et al.|2008 | and 
Asp-Clean ( Bhatnagar & CornweLl||2004 1 algorithms 
are advanced image deconvolution techniques for imag- 
ing fields with extended emission. The MS-MFS algo- 
rithm ( |Rau & Cornwell||2010| ) for wide-band multi- 
frequency synthesis not only deals with the deconvo- 
lution of extended emission, but also accounts for spa- 
tially resolved spectral index variation across the field 
of view. The combination of AW-Projection and MS- 
MFS algorithms can therefore be used to account si- 
multaneously for all the dominant wide-field wide-band 
effects which are expected to limit the imaging perfor- 
mance with SKA-pathfinder telescopes. 

The algorithmic design of the combined AW-Proj- 
ection and MS-MFS algorithm is such that it falls in the 
category of embarrassingly parallel algorithm. To miti- 
gate the fundamentally higher computing load, work is 
in progress to deploy these combined advanced algo- 
rithms on High Performance Computing (HPC) plat- 
forms, typically consisting of a cluster of computers 
with multi-core CPUs and high bandwidth intercon- 
nect (Bhatnagar et al.'2009'). 

Since these algorithms also need to iterate over a 
large volume of data, computing clusters connected to 
a parallel file system (e.g. the Lustre file system) and 
using multi-threaded I/O techniques are currently be- 
ing used. Tests with up to a few Terabytes of data 
show close to linear scaling in computing with num- 
ber of computing nodes, although this is subject to the 
achievable compute-to-I/0 ratio. Work is in progress 
to address the problem of improving this ratio while 
keeping the resulting memory footprint in reasonable 
limits. For SKA-sized data sets, iterating over the en- 
tire data set might not be feasible. Work is in progress 
to devise computing schemes and algorithms that may 
eliminate or at least reduce the number of iterations in- 
volving the entire data. 

4.4.2 Development of new deconvolution algorithms 

In the last few years, the theory of Compressed Sens- 



ing (CS: Donoho 2006 Candes & Romberg & Tao 



2006 1 has been developed by the signal/image process- 
ing community, and may be applicable to radioastro- 
nomical image deconvolution. 
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The theory of CS mostly addresses random sam- 
pling, sparse signals, and coarse-grained dictionaries. 
In radio interferometry, sampling is usually structured, 
astronomical images are not sparse, and dictionaries 
must be fine-grained to adequately represent the im- 
ages. Theoretical CS results therefore do not usually 
improve estimates of astrophysical sources from in- 
terferometric data. However, a key ingredient of the 
CS framework, sparse representations ( Mallat|2()08 



gies may be computationally and qualitatively more ef- 
ficient for large dictionaries (Sta rck et al. [2010 1, the 



has demonstrated significant improvements in radio- 
astronomical deconvolution, and has been used as a 
tool for deconvolution of extended and diffuse sources 
in the image plane (|Wiaux et al. 2009a| |Li et al. 



|2011a[ [Dabbech, Mary, & Ferrai-i,2012| , as well as in 
Faraday depth ( |Li et al. |201 lb| l. 

Sparse representation theory ("Fornass ier|2010| l may 
be seen as a generalization of ideas exploited in CLEAN 



(e.g.|Hogbom||1974[rWaBcer & Schwarz|[T988l |Con> 
[well et al.|2008| l. CLEAN performs an iterative decon- 
volution by using a dictionary of shifted point spread 
functions (PSFs) to construct a sparse dictionary of 
point sources. Sparse representation theory aims to 
construct a dictionary of more complex geometrical 
features than point sources, addressing three fundamen- 
tal issues of astronomical image restoration. First, the 
sparsity of the dictionary reduces the indeterminacy 
and instability caused by the zeros of the transfer func- 
tion of an interferometer. Second, the theory allows 
flexible and sophisticated models that can cope with 
complex astrophysical sources. Third, it offers compu- 
tationally efficient optimisation techniques to solve the 
resulting deconvolution problem. 

Sparse representations can take either of two ap- 
proaches: synthesis and analysis ( |Elad et al.|2005| l. 

In the synthesis approach, the unknown intensity 
distribution x (of size {N, 1), say) is assumed to be 
sparsely synthesisable by a few atoms of a given full 
rank dictionary S of size {N, L). Hence, we write x 
as X = S7, where 7 (the synthesis coefficients vector) 
is sparse. In contrast, the analysis approach assumes 
that X is not correlated with some atoms of an analysis 
dictionary A of size {N, L): A'^x is sparse. 

Since real images can often be approximated by 
a linear combination of a few elementary geometrical 
features, the synthesis approach is more intuitive and 
has been the focus of more research. Its design sim- 
plicity has also made it popular in image processing 
applications such as compression, de-noising, and in- 
painting. However, because a synthesis-based decon- 
volved image is restricted to a subspace of the synthe- 
sis dictionary, the resulting astrophysical images may 
be too rough to be realistic. On the other hand, since 
the signal is not built from a small number of atoms, 
the analysis approach may be more robust to "false de- 
tections" . However, the number of unknowns in the 
synthesis approach (the number of atoms in the dic- 
tionary) may be much larger than in the analysis ap- 
proach. Thus, while analysis-based optimisation strate- 



question of which approach is best remains open ( Gri- 
[bonval & Schnass|2009| and references therein). 

Regardless of the approach chosen, a sparsity -based 
deconvolution method still requires a dictionary, cho- 
sen from a class of images ( Mallat|2008 1. Astronom- 
ical wavelet dictionaries are widely used, but some- 
times fail to represent asymmetric structures adequately. 
In such cases, it may be necessary to use other trans- 
forms designed for specific classes of objects, such as 
curvelets, which are optimum for curved, elongated 
patterns such as planetary rings or galaxy arms, and 
shapelets, which are often optimum for galaxy mor- 
phologies. Modelling complex images may require sev- 
eral dictionaries to be concatenated into a larger dic- 
tionary (Chen & Donoho & Saunders] ( |1998| l; |GrP] 



bonval & Nielsen ( 2003 1), although computational is- 



sues limit the size of dictionaries, especially for radio 
synthesis images containing hundreds of thousands of 
Fourier samples. 

Several deconvolution methods that combine these 
approaches have recently been suggested (e.g.|Suksmono| 



'2009; Wiaux et al."2009a b ; Vannier et al. 2010 ; Wengej 
,et al. 2010, McEwen & Wiaux 2011)), including suc- 
cessful simulations of SKA pathfinder observations. O 



et al. (2011al) present a classical synthesis approach 



with an lUWT (Isotropic Undecimated Wavelet Trans- 



form) synthesis dictionary, Carillo, McEwen, & Wiaux 



(2012| used an analysis approach using a concatena- 



tion of wavelet bases, and Dabbech, Mary, & Ferrari 
(2012) defined a hybrid analysis-by-synthesis approach 



While the theoretical results of CS have not yet 
generated tools which can be used routinely for synthe- 
sis imaging, they have brought new perspectives from 
the domain of sparse representations: improved spar- 
sity models, and improved optimization algorithms. The 
increased research effort on fundamental models for 
sparse representations (analysis / synthesis / hybrid mod- 
els) leads to improved models and thus to improved 
reconstruction. Research in sparse representations is 
active and growing, and promises valuable future de- 
velopments in radio astronomical deconvolution algo- 
rithms. 

4.5 Source Extraction and Measurement 

There are many approaches to source extraction, each 
with different strengths and weaknesses, so it is impor- 
tant to mount a joint effort to determine how an opti- 
mum source finder would work, and what algorithms 
to implement. In such a joint effort, a reference data 
set should be subjected to the various algorithms, for 
optimum comparability. 

Building a source finder is likely to be subject to 
compromise, and current source finders fall well short 
of what may be theoretically achievable. For example, 
no source finder currently accounts for variations in the 
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point spread function across the image. Development 
of an "optimum" source finder may well require col- 
laboration with computer scientists with no astronomy 
background to develop signal extraction algorithms. It 
is likely that no single algorithm will be optimum for 
all scenarios. 

4.5.1 Compact Source Extraction 

Two main areas of recent focus have been investiga- 
tions into the impact of background and noise estima- 
tion on the detection of compact sources, and a rig- 
orous comparison of existing source finding software 
tools. The first of these focuses on source detection, 
where the background and noise levels in an image are 
characterised ( [Huynh et al.|2012 l 



sation process. 

These issues have triggered the development of new 



The second focuses on source characterisation, where 
islands of pixels are described or fit in some well de- 
fined manner ( [Hancock et al.|[20T2l ). The most com- 
mon way to describe a source in radio astronomy is as 
a group of Gaussian components. 

The process of background characterisation has been 
explored bylHuynh et al.J ( |2012p , who tested the steps 
of background and noise measurement, and of the choice 
of a threshold high enough to reject false sources yet 
not so high that the catalogues are significantly incom- 
plete. This analysis explored the results from testing 
the background and thresholding algorithms as imple- 
mented in the SExtractor, Duchamp/Selavy, and SFIND 
tools on simulated data. The result of this analysis 
shows clearly that it will be crucial to first develop and 
then implement an automated algorithm for establish- 
ing the appropriate scale size on which to estimate the 
"local" background and rms noise values for each im- 
age, and that this may not be common to all images. 
This analysis also suggests that the False Discovery 
Rate approach ( Hopkins et al.|[2002 i to thresholding 
may be the most robust, optimising for both complete- 
ness and reliability. 

The source characterisation stage has been found 
to have its own distinct issues ( [Hancock et al.|20"T2l l. 
For a selection of commonly used source finding tools 
(SFind, Sextractor, IMSAD, and Selavy) there are sig- 
nificant differences in the way that islands of sources 
are divided into components and subsequently fit. This 
is exacerbated with examples of sources which are blen- 
ded or adjacent, especially when there is a large dif- 
ference in the relative flux densities (even though the 
fainter component(s) may in their own right be very 
high signal-to-noise objects). Such behaviour has par- 
ticular implications for variable objects. If one such 
source were variable or transient, this behaviour could 
be missed if a source finding program isn't able to cor- 
rectly decompose the island. If a source finder's ability 
to decompose an island into components is related to 
the clipping limit, or background noise, then a source 
of constant flux could be interpreted as transient or 
variable due to irregularities in the source characteri- 



algorithms (e.g. Hancock etal. 2012 Hales etal. 2012a I 
that aim to determine the correct number of compo- 
nents within an island, and to generate an initial pa- 
rameter set for the fitting of multiple gaussian compo- 
nents. Hancock et al. ( 2012[ l use a Laplace transform, 
to identify regions of negative curvature in the pixel in- 
tensity profiles, combined with the map defined from 
the islands of pixels above a threshold, to automate the 
process of identifying the numbers of components in 
each island, and their initial position and peak flux den- 
sity estimates. 

Another clear requirement that has developed from 
these analyses is the need to quantify the false detec- 
tion rate of sources as a function of signal-to-noise 
level, and perhaps other image or telescope parame- 
ters, in order to assess the likelihood that any detected 
source is actually real. The impact of even a 0. 1 % false 
detection rate when catalogues of many tens or hun- 
dreds of millions of sources are being anticipated, is 
clearly large. 

These initial steps have demonstrated that while 
source-finding has been ubiquitous in radio astronomy 
for many decades, there is still scope for optimising 
and quantifying the process, and that careful attention 
will need to be paid to these details in the lead up 
to the new generation of radio telescopes. There will 
be many lessons to be learned in order to ensure that 
source measurement using the SKA, ultimately, can be 
made as robust as possible. 

4.5.2 Extended Source Extraction 

Many of the feature recognition algorithms widely ap- 
plied in fields such as optical character recognition, 
machine vision or medical imaging are absent from the 
repertoire of techniques routinely employed by astro- 
physicists. However, the increasing complexity and 
scale of astronomical data are driving an interest in 
more advanced feature extraction techniques for astro- 
nomical contexts. In particular, radio surveys on next 
generation telescopes will generate image data con- 
taining tens of millions of astronomical sources, many 
of which will be extended. 

Localised kernel transforms, such as the Wavelet 
transform, are often used for astronomical purposes, 
including diffuse source detection (e.g. detection of 
diffuse X-ray emission in the COSMOS field: Finoguerjov 
et al. 2006) . They are created either by modifying a 
conventional kernel transforms, or by using exotic ba- 
sis functions, and are supplemented by a localisation 
algorithm to allow the target to be located. The suc- 
cess of these methods relies on "matching" the basis 
function to the problem at hand. Given the variety of 
source morphologies in astronomical images, no single 
transform will adequately detect all sources. Localised 
kernel transforms work best in cases where a single 
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class of objects can be matched to a suitable set of ba- 
sis functions, and can also be used to separate diffuse 
from point source components using curvelet, ridgelet 



or a trous transforms (e.g. Fig. 24 1 
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Figure 24 Feature extraction in the HST image of 
Abell 370 from Starck et al."| (p003]l. The top left 



panel shows the original HST image, top right shows 
the co-added image from the ridgelet and the curvelet 
transforms which extract extended features, bottom left 
shows the reconstruction from the a trous algorithm 
which responds to point sources, and bottom right 
shows the addition of the results of all three transforms. 



Compressive processes rely on the fact that pix- 
els in an image are not independent, so it is possi- 
ble to encode the patterns in the image rather than the 
pixel-by-pixel information. This provides a more ef- 
ficient "compressed" representation of the information 
by providing a 'model' of the image. Object detection 
is then applied to this 'model' rather than to the real 
data, greatly reducing the complexity of the problem. 
Compression attempts to preserve the information in- 
herent in the original image but does not attempt to ab- 
stract geometric or semantic information. Compressive 
techniques therefore preserve large-scale noise features 
in addition to diffuse sources. The use of compres- 
sive sampling for deconvolution has been discussed in 



U.4.2 



Template matching is a process where an image is 
correlated with a target object of known form, and is 
similar to "matched filtering". This technique has been 
used in radio astronomy for pulsar timing prediction 
( |van Straten ||2006| [Oslowski et al.~||2011| l, transient 
source detection (Trott et al. ||2012 % and to remove 
radio point source contaminants from diffuse source 



backgrounds ( |Pindor et al. |201 1 \. Two disadvantages 
of template matching are (a) it is hard to control the 
required number of filters, resulting in a greater num- 
ber than for localised kernel transforms, and (b) each 
template must fit a range of orientation and scaling. 
Even when the computational complexity of template 
matching is reduced by combining it with a scale in- 
variant transform such as the Mellin transform, orien- 
tation must still be calculated making this technique 
very computationally expensive. On the other hand, it 
has the advantages that (a) the templates do not have 
to be precise, and cartoon-like models suffice, (b) the 
method is robust against artefacts, (c) it can be applied 
in either the image or Fourier domains, and (d) it re- 
turns the size, orientation, location type of detected 
sources. 

Hough transforms can find and characterise geo- 
metrical objects such as lines, circles and ellipses. In 
general, the Hough transform is equivalent to template 
matching for some geometric shape, although it inher- 
ently incorporates scaling and orientation constraints. 
Circle or elliptical Hough transforms are particularly 
relevant to astronomical source detection to determine 
questions such as the size and location of a circle/ellipse, 
and have the advantage that the geometrical abstrac- 
tion limits the number of filters required (e. g. different 
types of objects contain circular and arc-like features). 
The Circle Hough Transform (CHT) can find partial 
circular objects and is robust to noise. The CHT in- 
puts an image and outputs a three dimensional array in 
which two dimensions represent the possible locations 
for the centres of circles and the third dimension spans 
the set of possible circle radii. Peaks in this so-called 
"Hough space" therefore correspond to circles in the 
input image. In its simplest form the circle Hough 
transform is computationally challenging, with both 
computational effort and memory consumption scaling 
as 0{n'^) when transforming an n x n pixel image. Re- 
cently the transform has been recast via a convolution 
approach which reduces the complexity to 0{n^) log n 
( IHolUtt 12009] [20r2 ). Although this is not the fastest 
approach for astronomical source detection, it has been 
applied as a source detection method for a variety of 
extended radio sources with arc-like features includ- 
ing supernova remnants, tailed radio galaxies and radio 
reUcs (Hollitt & Johns ton-Holhtt ||2009| |20T2| . For 
example. Figure 25 shows the detection of a supernova 
remnant. Hough transforms return semantically useful 
information such as position, size, orientation and ec- 
centricity, but will produce somewhat less information 
than a full template matching process. For example, it 
cannot distinguish between a supernova remnant and a 
tailed radio galaxy. 

Segmentation methods assume nothing about a source 
but rather seek to group pixels together based on dif- 
ferences with the background. Three main classes of 
segmentation methods are used in astrophysics: region 
growing, edge detection, and waterfalHng. The sim- 
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Figure 25 Application of the circle Hough transform 
to detect radio supernova remnants in the Molonglo 
Galactic Plane Survey. Left to right: the original image 
of G315.4-2.3, the location and size of the remnant as 
found by the CHT, the original image with ten times 
the Gaussian noise and the response of the CHT in the 
noise case (yellow) compared to the original (black) 
which demonstrates the robustness of CHTs to noise. 
( IHolhtt & Johnston-Holhtt |2012| l. 



plest version of a segmentation method is to set a thresh 
old and assume everything above that threshold is a 
target. In practice adjoining pixels above the threshold 
are grouped together as a single target. This approach 
is usually too simplistic and requires an extension via 
one of the three subclasses above. Region growing, 
which includes floodfilling techniques, compares ad- 
jacent pixels at the edge of a target and groups them 
together based on some tolerance in pixel value (e.g. 
Ducham p and BLOBCAT Whi ting |2012[ [Hales et aL 
2012b] l). This technique alone is poor at dealing with 



disjointed areas, although it may be heuristic ally aug- 
mented to define nearby areas as one object. Another 
method, edge detection, uses thresholds to find the bor- 
ders between different "types" of pixel. This is rarely 
useful in astronomy because sources have "soft' edges 
. Waterfalling involves considering the image as a to- 
pographic surface and calculating the intensity gradi- 
ent, pixels with the lowest gradient denoting watershed 
lines which are used to bound regions. An advantage 
of segmentation methods is their speed, making them 
good for large area surveys. One of the downfalls of 
segmentation methods is that no semantic information 
is returned. An additional classifier step is required to 
be applied to the output of segmentation algorithms to 
determine semantic properties. 

While algorithms for the detection of point sources 
are well-established, automatic detection and charac- 
terisation of diffuse sources presents a significant chal- 
lenge to large radio surveys. While methods such as 
local kernel transforms, template matching and Hough 
transforms have been demonstrated to detect and char- 
acterise diffuse sources, they are too computationally 
expensive to apply to large images. The favoured ap- 
proach at present is therefore a combination of seg- 
mentation algorithms with heuristics to extract diffuse 
source catalogues. This combined process is not likely 
to be sufficient for future survey science goals and some 
combination of techniques, possibly applied on differ- 



ent resolution scales, will be required. 



4.6 Classification and Cross-Identification 

The large area surveys proposed by the SKA pathfinder 
projects will return many tens to hundreds of millions 
of sources, requiring automated techniques both to group 
components into sources and to identify the counter- 
parts at other wavelengths. 

Table[2]lists the current, or planned, large area mul- 
tiwavelength surveys that have significant overlap with 
next generation radio surveys. Based on analysis of 
current deep multiwavelength fields such as COSMOS 
and GOODS-N (|Scoville et al.|2007HSchinnerer et"aL] 



[2007lpomson et al.|2010||Giavahsco et al.|2004 l), we 
estimate that, assuming a limit of rms ~ 10 /xJy/beam 
at 1.4 GHz, ^ 50% will have counterparts in areas 
co-incident with deep optical surveys (i.e. DBS, Pan- 
starrs), while ~ 70% will have near-IR counterparts in 
the overlap regions with surveys such as VISTA VHS 
and VIKING. However, because most of the deeper 
surveys only cover small areas, about 50% of radio 
sources will have no counterparts at other wavelengths, 
with this fraction falling to at most 30% in the ^ 1000— 
2000 deg^ overlap with the deepest large area optical- 
near IR surveys. 

[Padovani] ( |201 ip has estimated cross-identification 
rates for the SKA down to nJy radio flux densities. Fig- 
ure 26 plots i?mag versus 1.4 GHz radio flux density 
and shows the expected i?inag for various classes of 
sources that will be found in the radio surveys. 



Figure 26 shows that, down to 50 /iJy, ~ 50% of 
SFGs should be detected by the SDSS (northern hemi- 
sphere) and SkyMapper (southern hemisphere), while 
most AGN wiU be fainter. PAN-STARSS and the LSST 
should detect most of the radio sources down to this 
flux density but on longer timescales (e.g. by ^2028 
in the case of the LSST, assuming a start of opera- 
tions in 2018). Surveys reaching 1 /iJy (e.g. those from 
the VLA and MIGHTEE) will obviously have much 
fainter counterparts and might require JWST and the 
Extremely Large Telescopes (ELTs) for the optical iden- 
tification of most of their sources, which however will 
be covering a relatively small field of view (up to a 
few arcmin^). X-ray surveys are intrinsically much 
less sensitive than optical/IR surveys, and will detect 
only a very small fraction of the AGNs detected by ra- 
dio surveys ( Padovani [201 1\ . 

The requirements for spectroscopic redshifts are ob- 
viously even more difficult. For example, long ex- 
posures (~ lOh) with 8/10 m telescopes can secure 
redshifts in the case of strong emission lines down to 
-Rmag ~ 26. Large-area spectroscopic surveys, how- 
ever, cannot reach these magnitudes. For example, cur- 
rent magnitude limits are r < 17.8 (galaxies) and i < 
19.1 (quasars) for the SDSS and bj < 20.85 for the 
2dF QSO redshift survey. 
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to be built into the process for selecting the most likely 
candidate. This can in principle even extend to includ- 
ing the knowledge that the core of an FRII is associated 
with a host galaxy while the lobes are not. This bal- 
ance between reliability of matches, and completeness 
as a function of flux, will be key in deciding which al- 
gorithm is appropriate for the scientific exploitation of 
radio surveys. 

Finally, it is likely that any automated algorithm 
will fail on complex radio sources, which may include 
those which are most scientifically interesting. For these 
cases, a "Citizen Science" approach is being developed 
to enlist the help of thousands of enthusiastic amateurs. 



in collaboration with the Galaxy Zoo group ( Lintott et 



al.|2008 1 to examine each of the millions of complex 
sources by eye. The "Radio Zoo" project expects to 
release its first prototype in mid-2012. 



Figure 26 i?mag versus the 1 .4 GHz radio flux density 
for faint radio sources. The diagonal dashed line indi- 
cates the maximum value for SFG and the approximate 
dividing line between radio-loud and radio-quiet AGN, 
with SFG and radio-quiet AGN expected to populate 
the top left part of the diagram. The typical R mag- 
nitudes of the three classes at 5*1 4 = 1 /iJy are also 
shown, with SFG split into starbursts and spirals. Fi- 
nally, the mean radio and i?mag values for sources from 
the VLA-CDFS sample, with error bars indicating the 
standard deviation, are also marked. The horizontal 
dot-dashed lines indicate the approximate point-source 
limits of planned and existing surveys. 



The source density of potential matches will proba- 
bly be > 1 per radio beam, so simply taking the nearest 
neighbour to the radio source in other catalogues is not 
a feasible approach. Furthermore, radio sources will 
typically have more than one component per source, 
so a simple nearest-neighbour algorithm will attempt 
to match sources to the lobes of an FRII galaxy as well 
as to the core. Clearly, more sophisticated matching 
algorithms must be considered 

A significant number of techniques exist to facili- 
tate the matching of catalogues. Techniques such as the 
p-statistic 



( Downes et al.|198 6 ) and the likelihood ratio 
( [Sutherland & Saunders] 1992[ ) make use of the number 



density of background (and true counterpart) sources 
to "weigh" potential counterparts within a given search 
radius, but cannot take into account the likelihood that 
several radio components may correspond to one op- 
tical component. More recently Budavari & Szalay 



(2008 1 outlined a Bayesian method to determine the 
most likely match. This method, while similar to the 
likelihood ratio, is distinct in that it does not consider 
the source density, and hence is less biased towards 
faint counterparts (although potentially less reliable for 
very high source densities). A Bayesian approach ( |Bu-| 
|davari|20TT) also allows a number of other constraints 



4.6.1 The Likelihood Ratio technique for radio 
source cross-identification 

The Likelihood Ratio (LR) technique (e.g. 'Richterj 



1975 Sutherland & Saunders 1992; Cilegi et al. 2003 



3] 



can be used to identify optical/near-infrared counter- 
parts with sources selected from low resolution survey 
data. Smith et al. (2011 1 discuss a recent applica- 



tion of this Bayesian technique to identify 2423 op- 
tical counterparts to 6621 sources from the Herschel- 
Astrophysical TeraHertz Large Area Survey (//-ATLAS; 
[Bales et al.|2010| l, selected at 250 /im. 

The LR technique relies on the brightness and as- 
trometric properties of the input catalogues in general, 
and each individual source in particular, to determine 
the ratio of the probability that two objects are related, 
to the probability that they are unrelated. After us- 
ing an empirical method to account for the fraction 
of sources which are not detected in the optical/near- 
infrared catalogue, we can determine the reliability of 
each association, with values ranging from (not re- 
lated) to 1 (a match). 

The likely positional uncertainties, coupled with 
the probable optical/near-infrared catalogues that will 
be available for identifying counterparts to the new cat- 
alogues that these surveys will generate, suggest that, 
apart from the problem of radio sources with multi- 
ple components, the problem will be very similar to 
that addressed using the LR technique in the Herschel- 
ATLAS. 



'McAlpine et al. 



( 2012 1 have begun detailed inves- 
tigations based on existing deep radio survey data and 
new near-infrared data from the VISTA VIDEO survey 
to quantify how effective the LR technique will be for 
these purposes, and how it will be affected by using 
data of different sensitivity and spatial resolution. 



4.6.2 Bayesian approaches 

Budavari & Szalay ( 2008| l propose Bayesian hypoth- 
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esis testing for cross-identification. Given a set of de- 
tections one can ask directly whether the measured di- 
rections are consistent with a common single object. 
This hypothesis is tested against its complement using 
the Bayes factor, which is the likelihood ratio 



B = 



(13) 



where the likelihoods are calculated as a sum over all 
possible model parameters, e.g., the true position of the 
source. Assuming the astrometric uncertainty is well 
approximated with a normal distribution, the deriva- 
tion is analytic. For example, in the case of two de- 
tections of point sources with high-precision circular 
uncertainties and an all-sky prior on the true position 



B 



cxp 



1^2 



2(a? 



(14) 



where tp is the angular separation of the detections and 
(Ti are the uncertainties; all in radians. When this di- 
mensionless quantity is larger than 1, the data support 
the association. The higher the value, the greater the 
evidence. The probability is analytically given in terms 
of B and the prior of the hypothesis, which is a func- 
tion of the source densities of the input data sets. The 
approach is applicable to different kinds of data, such 
as fluxes or polarisation. The corresponding Bayes 
factors are computed independently and simply mul- 
tiplied together. The method's firm statistical founda- 
tion along with its explicit dependence on the geometry 
and physical modelling provides a clean framework for 
catalog associations at all wavelengths and straightfor- 
ward extensions. Such examples include the associa- 
tion of stars with unknown pro per motions (jKerekes| 
|et al.|2010| or transient events ( |Budavari|2012| l. 

4.7 Measuring Redshifts 

To obtain the full value of a radio survey, it is essential 
to have multi wavelength identification, and redshifts. 
However, only about 1% of radio sources from these 



surveys will have spectroscopic redshifts (§ 4.6 1, and 
so it is necessary to rely on photometric redshifts, or 
even statistical redshifts. 

The most popular way to compute photometric red- 
shifts is via fitting of templates of galaxies to the mea- 
sured multi -band photometry of target galaxies. Clearly, 
the accuracy of the result depends not only on the li- 
brary of templates, but also on the number of pho- 
tometric points available and their associated errors. 
Other techniques, such as k-nearest neighbour (kNN) 
implementations. Support Vector Regression (SVR) 
models, Self-Organised maps (SOMs), Gaussian pro- 
cessing (GP), neural-networks (ANNz) and rainbow 
forest (RF) do not make any assumptions about the 
spectral energy distribution (SED) of the galaxies, or 
the type and amount of extinction. Instead, they build 



on a simple colour-redshift relation using a training 
sample of sources with reliable spectroscopic redshifts. 
These methods are very efficient and fast whenever 
the training sample is representative of the entire pop- 
ulation, and become less reliable outside the redshift 
range for which they were trained. 

About 20% of EMU sources will have photometric 
data from the 6 optical bands of the Skymapper survey, 
the 2 near-infrared bands of the Vista Hemisphere Sur- 
vey (VHS) and the 4 mid-infrared bands of WISE (see 
Table [2]l. Typically, with these available datasets, an 
error in redshift of ^ 0.03 can be reached for normal 
galaxies up to redshift z ^ 1.5 using either standard 
(SED template fitting) or empirical (machine learning) 
techniques. 

About 20% of identified sources will be classified 
as AGN, based on one or more of the following criteria: 
(a) radio morphology, (b) radio polarisation, (c) X-ray 
detection, (d) a radio spectral index very different from 
a ^ —0.7, (e) a radio/IR ratio very different from the 
canonical star-forming value. About 5% of the sources 
will be QSOs, detectable by the next X-ray all sky sur- 
vey, e ROSITA ( |Cappelluti et al.|2011|perloni et al. 



2012 1, and for these sources a special treatment must 
be adopted ( |Salvato et al.|2009|[20TT] i. 



About 80% of the EMU sources will have an opti- 
cal counterpart that is fainter than r = 22 magAB and 
for these sources the ancillary multiwavelength data 
set will be coarse and non homogenous. Finally, the 
training sample with spectroscopic redshift available at 
high redshift may not be as complete as at z < 1, ham- 
pering the empirical photometric redshift techniques. 

With the final goal of providing the most reliable 
photometric redshift for the larger number of sources 
possible, the EMU redshift working group is testing 
all the methods mentioned above, trying to understand 
the advantages and disadvantages of each of them, and 
thus applying the most suitable, or a combination, of 
the methods, depending on (a) the available depth of 
the multiwavelength data and (b) on the redshift range 
available for the training sample. We started experi- 
menting on COSMOS field |Scoville et al.| ( (2007] l rep- 
resentative of a pencil beam survey ( 2 deg'^) with a 
very deep and homogenous multiwavelength data set, 
but for which the training spectroscopic sample is lim- 
ited to low redshift radio sources, and with a strong 
component of X-ray detected sources. Other fields which 



we will be testing are the larger area ATLAS (Nor- 



ns et al. 2006 Middelberg et al. 2008a i which has 



deep but non homogenous ancillary data, and GAMA 
( 310 deg^) for which the training spectroscopic sam- 
ple is rich but representing only the nearby universe 
(r < 20 mag) and with limited availability of ancillary 
data. 

In addition to the optical, near, and mid-infrared 
data, we will make use of the radio properties of the 
sources. While for the empirical methods the radio flux 
will be used as an additional band, for the SED fitting 
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technique, we will use the radio properties as a prior 
for limiting the redshift range or the templates used for 
the possible solutions. 

The "statistical redshift" concept, introduced by |No} "- 
" ' recognises that a probabilistic estimate of 



ns 



(2011 



redshift for a large sample can be useful, even if many 
of them are incorrect, provided that the level of in- 
correctness and incompleteness can be accurately cali- 
brated. Even a non-detection can carry useful informa- 
tion, and radio data themselves can add significantly to 
the choice of SED template. For example, 

• high redshift radio galaxies can be identified from 
their strong radio emission coupled with a NIR 
non-detection ( jWillott et al.|2003t[Norris et al.| 
[2011a, ). 

• a steep radio spectral index increases the proba- 
bihty of a high redshift pe Breuck et al.|2002| l. 



• the angular size of a particular galaxy class can 
be loosely correlated with redshift (Wardle & 
|Miley| 1974) 1, 

• polarised radio sources are nearly always AGN 
( Hales et al. |2012b| ), while unpolarised sources 
are mainly star-forming galaxies. Therefore, for 
the EMU/WODAN surveys, polarised sources have 
a mean z ^ 1.8, while unpolarised sources have 

a mean z ^ 1.1. 



4.8 Data Issues 

4.8.1 Stacking and Data-Intensive Research 

Stacking is used to explore the properties of a class of 
objects which are below the detection limit of a survey. 
For example, the radio flux of RQQs can be measured 
by averaging the measured radio flux at the position of 
each RQQ in a radio survey, even though each individ- 
ual RQQ has a flux well below the detection limit. 

More generally, the process of stacking involves 
combining (typically by taking a censored mean or me- 
dian) the data at the position where such objects are ex- 
pected in the survey. The noise tends to cancel, while 
any low level of flux density in the sources adds, re- 
sulting in a detection threshold very much lower than 
that of the unstacked survey. Stacking at radio wave- 
lengths has been used very successfully (e.g.jBoyle et 



al.|2007[[mite|20071|Ivison et al.|2007| [20101 [Dunne 



et al.|2009t|Messias et al.|20T0l [Bourne et al.|201 l| l on 
high-resolution data for the purpose of studying faint 
populations which are below the detection threshold 
of the radio image, and has proven to be a powerful 
tool for studying star formation rates, AGN activity, 
radio-infrared correlation, and measuring the fraction 
of the extragalactic background contributed by vari- 
ous source populations. Jack-knifing techniques are 
currently being explored to provide more information 
about the flux density distribution of stacked sources 



than is available from a simple mean or median (Rees 
et al., in preparation). 

The unprecedented area-depth product of SPARCS 
surveys makes them very suitable for stacking. For 
example, stacking a sample of a million optically se- 
lected galaxies in the EMU data will result in a noise 
level of ^ lOnJy. Because of the wide area of the 
EMU, WODAN and LOFAR surveys, even rare classes 
of source can be stacked, and it should also be possi- 
ble to create stacked images of extended sources, such 
as clusters. However, the extent to which such deep 
stacking will be successful will depend on the extent 
to which imaging artefacts are cancelled by stacking. 
Stacking is also not without its hazards, and can be 
biased by a number of effects such as failure to ac- 
count for the point-spread function, particularly in the 
presence of confusion and in highly clustered fields 



(Bourne et al. 201 1[ Greve et al. 2010 




Penner et al.[ 


[ioi IfiChary et al.|2008| Bethermin et al 




20101. 



Other examples of data-intensive research include: 

• Identification of sources which do not fit into 
known categories of radio source, and so are likely 
to be artefacts or exciting new classes of source. 



as discussed in S 3.9 



• Extraction of low surface brightness radio emis- 
sion to detect the WHIM synchrotron emission 
from cosmic filaments or sheets, by cross-correlating 
with optical galaxies selected from other surveys. 

• Cross-correlation of low-redshift galaxies with 
high-redshift galaxies, or the CMB, to test cos- 
mology and fundamental physics as discussed in 
§[3l 



4.8.2 Survey versus Pointed data 

Since surveys compete for observing time with con- 
ventional "pointed" observations, it is useful to con- 
sider their role compared to projects doing detailed sci- 
ence on particular sources or fields Clearly both modes 
are valuable, and complementary. Surveys provide source 
lists for follow-up observations of particular sources, 
as well as finding and characterising calibrators for other 
radio observations. However, the huge volume of data 
from the new radio surveys will trigger a shift in the 
way astronomers "observe". The ingenuity currently 
applied to devising a key experiment, and crafting a 
successful proposal, will be applied to devising novel 
ways of mining the data. To some extent, this is already 
happening, and about three times as many citations are 
delivered by results from the HST/STScI archive than 
from the papers written by Pis. The surveys discussed 
here will also have an enormous legacy value, since 
they will be hard to improve on until the technology 
makes a further major step forward. 
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5 Conclusion 

The SKA Pathfinder Radio Continuum Survey work- 
ing group (SPARCS) has been estabUshed to coordi- 
nate the science and technical developments associ- 
ated with continuum surveys being planned or under- 
taken with the SKA pathfinder telescopes (APERTIF, 
ASKAP, VLA, e-MERLIN, e-EVN, LOFAR, Meerkat). 
This review paper was triggered by the first SPARCS 
workshop, and is an attempt not only to share the knowl- 
edge residing in the individual pathfinder projects, but 
also to estabUsh a baseUne of knowledge from which 
all projects can move on, and which can be used to 
optimise the development of the SKA. A primary out- 
come of the first SPARCS workshop was to rethink as- 
pects of the way we are constructing the SKA pathfinder 
telescopes and designing the surveys. For example, it 
became clear that the impact of these surveys on cos- 
mology and fundamental physics may dwarf the key 
science goals we have been working on for the last few 
years. That places demands on source-measurement 
accuracy and uniformity which in turn impacts on cal- 
ibration and imaging accuracy. 

Key demonstrable outcomes of the workshop in- 
clude this review paper, and the establishment of a se- 
ries of reference fields in both hemispheres which are 
now starting to be observed with major existing tele- 
scopes (VLA, ATCA, Westerbork) and will be observed 
with the new pathfinder telescopes (LOFAR, Aperitif, 
ASKAP, Meerkat) to ensure consistency and unifor- 
mity across the surveys. 

In this review paper, we have explored the current 
state of science to be tackled with the SPARCS sur- 
veys, and identified the most pressing challenges. We 
conclude that: 

• The SKA pathfinder telescopes and upgrades, and 
their planned continuum surveys, have comple- 
mentary capabilities which will lead to ground 
breaking science. Particularly exciting are (a) 
the abihty of EMU and WODAN to cover the en- 
tire sky with high sensitivity and resolution, de- 
tecting about 100 million galaxies, (b) the abil- 
ity of LOFAR to make complementary observa- 
tions of much of that region at low frequency,(c) 
the ability of POSSUM and BEOWULF to ob- 
tain detailed information about the magnetic sky 
over the same region, and (d) the ability of e- 
MERLIN, VLA, e-EVN, and MIGHTEE to study 
selected areas and objects with very high sensi- 
tivity and resolution. 

• Together, these surveys are likely to make ma- 
jor inroads on some of the pressing questions 
about the origin and evolution of galaxies, in- 
cluding distinguishing AGN and SF components 
of galaxies, tracing the evolution of the luminos- 
ity function for both, and exploring the causes 



and signatures of hot-mode and cold-mode ac- 
cretion. 

• Amongst the sources detected by the surveys will 
be very high redshift galaxies, possibly the high- 
est discovered, but identifying them, and mea- 
suring their redshift will be challenging. 

• The SPARCS surveys are likely to detect tens to 
hundreds of thousands of new clusters, and will 
also map the extended halo and reUc emission 
from a fraction of these. Comparison with X-ray 
observations is likely to yield new insights into 
the physics of clusters and the growth of struc- 
ture in the Universe. 

• The measurement of polarisation and rotation mea- 
sures is an important part of the surveys, and is 
likely to yield not only insight into the physics 
of the radio galaxies, but into the origin of cos- 
mic magnetism. There is a good chance that the 
intergalactic magnetic field will be detected and 
traced by these surveys. 

• Even without redshifts, several cosmological probes 
can use SPARCS survey data to measure cos- 
mological parameters and test models of grav- 
ity, with significantly more precision than cur- 
rent analyses. The combination of WODAN and 
EMU is likely to yield the best measurement yet 
of the Integrated Sachs- Wolfe (ISW) effect. The 
most powerfully constraining probe of dark en- 
ergy parameters turns out to be the source count 
correlations, while ISW provides the most strin- 
gent constraints on modified gravity theories. If 
photometric redshifts can be obtained for a sig- 
nificant subset of the radio sources, then these 
cosmological tests become even more powerful. 

• Variable and transient sources are currently poorly 
studied, and the SPARCS surveys are Ukely to 
change this field significantly. In addition to con- 
straining models of (e.g.) black hole accretion, 
these studies are also hkely to result in unex- 
pected discoveries. 

• While most SPARCS surveys are focussed on 
extragalactic science, they will also produce the 
best high-resolution data yet on the Galaxy, de- 
tecting essentially every SNR and UCHII region 
in the Galaxy, and vastly increasing the number 
of known radio stars. It is essential that this sci- 
entifically important aspect of the SPARCS sur- 
veys is not neglected in our drive to to answer 
the key questions of extragalactic astrophysics. 

• Because the SPARCS surveys are opening up a 
large area of unexplored parameter space, it is 
likely that they will make unexpected discover- 
ies. It is important not to leave this process to 
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chance, but to recognise it as a legitimate goal 
of large surveys, and to plan processes and soft- 
ware accordingly, to mine the deluge of data that 
these surveys will produce. 

To achieve these exciting scientific goals, many tech- 
nical challenges must be addressed by the survey in- 
struments. Given the limited resources of the global 
radio-astronomical community, it is essential that we 
pool our skills and knowledge. We do not have suf- 
ficient resources to enjoy the luxury of re-inventing 
wheels. We face significant challenges in calibration, 
imaging, source extraction and measurement, classifi- 
cation and cross-identification, redshift determination, 
stacking, and data-intensive research. As the SPARCS 
instruments extend the observational parameters we will 
face further unexpected challenges in calibration, imag- 
ing, and interpretation. If we are to realise the full 
scientific potential of these expensive instruments, it 
is essential that we devote enough resources and care- 
ful study to understanding the instrumental effects and 
how they will affect the data. The prime role of SPARCS 
is to facilitate the process of doing so by ensuring we 
share resources and expertise across the projects. 
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Table 2 Key multiwave length surveys with which EMUAV O DAN data will be cross-identified (restricted to surveys 
larger than 1000 deg^) adapted from Norris et al. (201 lb I. All magnitudes are in AB. The "detectable" column 
is the fraction of 1.4 GHz EMUAVODAN sources that are in principle detectable by the multiwavelength survey to 
its 5 cr limit. The "matched" column is the fraction of 1.4 GHz sources which are both detectable and in the area of 
sky covered by the multiwavelength survey. The sensitivity shown for the WISE survey is for the 3.5 /im band. 



Survey 


Area 


Wavelength 


Limiting 


Detectable 


Matched 


Data 


Name 


(deg2) 


Bands 


Mag. 






Release 








flux^ 


(%) 


(%) 


Date 


WISE^ 


40000 


3.4, 4.6, 12, 22 ^m 


SOfiJy 


23 


23 


2012 


Pan-Starrs^ 


30000 


g, r, i, z, y 


r < 24.0 


54 


41 


2020 


Wallaby3>b 


30000 


20 cm (HI) 


1.6 mJy'^ 


1 


1 


2013 


LSST^ 


20000 


u, g, r, i, z, y 


r < 27.5 


96 


48 


2020 


Skymapper^ 


20000 


u, V, g, r, i, z 


r < 22.6 


31 


16 


2015 


VHS^ 


20000 


Y, J, H, K 


K< 20.5 


49 


25 


2012 


SDSS^ 


12000 


u, g, r, i, z 


r < 22.2 


28 


8 


DR8 2011 


DES^ 


5000 


g, r, i, z, y 


r < 25 


71 


9 


2017 


VST-ATLAS^ 


4500 


u, g, r, i, z 


r < 22.3 


30 


4 


2012? 


Viking 1° 


1500 


Y, J, H, K 


K< 21.5 


68 


3 


2012 


Pan-Starrs Deep^ 


1200 


0.5 — 0.8, g, r, i, z, y 


g < 27.0 


57 


2 


2020 



^Denotes 5 a point source detection. However, in many cases, a priori positional information will enable 3 a data 
to be used, resulting in a higher detection rate. 

^Being an HI survey, WALLABY will measure redshifts for all detected galaxies out to z = 0.26. 
per 4km s~^ c hanne l achieved in 8 hours integration. 
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